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1.0  INTRODUCTION 


The  Articulated  Total  Body  Model  (ATBM)  computer  program  simulates  gross 
human  body  motion  in  response  to  internally  or  externally  applied  forces. 

Program  VIEW  was  written  to  provide  visual  computer-aided  data  analysis  fo.- 
users  of  the  ATBM  program.  It  plots  a  pictorial  representation  of  simula¬ 
tion-generated  body  position  data  at  any  constant  time  step.  These  plots, 
since  they  are  projected  views  as  seen  by  a  camera,  can  be  directly  com¬ 
pared  with  pictures  taken  from  tests  using  dummies  or  human  subjects.  The 
plots  from  the  VIEW  program  provide  an  efficient  method  for  preparing  and 
checking  initial  position  data  used  for  input  to  the  ATBM  program.  Still 
pictures  of  the  initial  position  of  the  pilot  and  cockpit  configuration  can 
be  overlaid  with  plots  from  program  VIEW  that  represent  body  position  data 
at  time  step  zero.  The  input  data  can  then  be  adjusted  until  the  desired 
positions  of  the  body  segment  contact  ellipsoids  are  obtained. 

Plots  from  program  VIEW  approximate  pictures  taken  with  a  camera  because  of 
the  projection  technique  used  by  the  program.  Three-dimensional  objects 
are  projected  through  a  point  onto  a  projection  plane.  This  is  similar  to 
3D  objects  being  projected  through  a  lens  onto  a  film  plane.  The  viewpoint 
used  by  program  VIEW  must  be  placed  at  the  corresponding  position  and 
orientation  of  the  camera  lens  to  reproduce  the  same  aspect.  Then,  by 
selecting  the  correct  scale  factor,  plots  from  the  program  VIEW  can  be 
overlaid  with  pictures  from  the  camera.  This  provides  a  convenient  means 
for  pictorally  comparing  simulated  and  photographic  experimental  data. 

The  basic  graphics  elements  used  in  program  VIEW  are  ellipsoids  and 
polygons.  These  basic  elements  correspond  to  the  ellipsoids  and  contact 
planes  used  in  the  ATBM  program.  The  VIEW  program  has  the  capability  to 
plot  convex  polygons  with  up  to  four  sides  which  allows  the  plotting  of 
realistic  figures.  The  graphic  elements  for  VIEW  are  defined  by  data 
stored  on  logical  unit  number  (LUN)  1  which  is  a  data  file  generated  by  the 
ATBM  simulation  program.  VIEW  reads  the  ATBM  input  file  to  obtain  linear 
position  and  orientation  data  for  the  elements  at  any  requested  constant 
time  step. 
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Each  element  Is  represented  by  a  set  of  contour  lines  that  are  projected  to 
form  an  Image  on  the  projection  plane.  Each  contour  line  consists  of  a  set 
of  connected  vectors  In  three-space  that  projects  to  a  set  of  connected 
vectors  on  the  projection  plane. 

The  first  step  In  generating  the  3D  ellipsoid  Images  Is  to  define  contour 
lines  for  ellipsoids.  This  is  done  by  setting  up  a  grid  mesh  In  the  local 
X-Y  plane  of  the  ellipsoid.  The  Z  values  for  the  contour  vectors  are  cal¬ 
culated  us'ng  a  constant  X  value  and  changing  Y  by  constant  Increments. 

This  results  In  ellipsoidal  contour  lines  that  are  concentric  about  the 
local  X  axis  of  the  ellipsoid.  Since  these  contour  lines  are  fixed  to  the 
local  coordinate  system  of  the  ellipsoid,  any  rotation  of  the  ellipsoid 
will  be  seen  as  rotation  of  the  projected  contour  lines.  Polygon  contour 
lines  are  defined  by  a  set  of  vectors  in  the  global  reference  frames,  each 
of  which  represents  a  side  of  the  polygon.  These  vectors  are  projected  in 
the  same  manner  as  ellipsoidal  vectors.  The  result  is  a  plot  of  the  pro¬ 
jected  contour  of  the  polygon. 

Hidden  line  routines  are  included  in  the  VIEW  program.  These  routines 
eliminate  sections  of  contour  lines  that  are  hidden  from  a  viewer  posi¬ 
tioned  at  the  viewpoint.  Some  lines  are  hidden  because  they  are  on  the 
surface  of  an  ellipsoid  that  faces  away  from  the  viewer.  Other  lines  are 
hidden  because  another  element  blocks  that  contour  line  from  the  viewer. 

The  routines  that  check  for  hidden  lines  are  formulated  to  handle  elements 
that  are  imbedded  in  other  elements.  For  example,  an  ellipsoid  can  be 
imbedded  in  another  ellipsoid.  The  contour  lines  will  be  drawn  only  up  to 
the  intersection  of  the  surfaces. 
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2.0  INSTALLING  THE  VIEW  PROGRAM 


The  VIEW  program  source  code  contains  1717  lines  of  FORTRAN  code, 
consisting  of  a  mainline  and  36  subroutines  and  functions.  The  subroutines 
and  functions  are  in  alphabetical  order.  VIEW  was  originally  developed  on 
a  CDC  6600  and  has  since  been  Implemented  on  a  Perkin-Elmer  3242  computer 
and  an  IBM  370  computer.  The  version  used  for  development  of  this  document 
was  the  Perkin-Elmer  version.  All  Calcomp  plotting  calls  are  standard 
except  the  call  to  the  Initialization  subroutine  PLOTS  which  will  be  dis¬ 
cussed  later.  All  other  subroutines  and  functions  are  either  contained 
within  the  program  or  are  standard  Intrinsic  FORTRAN  functions  (SQRT,  COS, 
etc. ). 


2.1  POSSIBLE  PROGRAM  MODIFICATIONS 

Although  the  VIEW  program  source  code  may  compile  and  execute  properly  on 
the  target  system,  slight  modifications  may  be  necessary  due  to  compiler 
differences.  They  include  the  following: 

a.  In  certain  FORTRAN  compilers,  seven  character  subroutine  names 
are  allowed.  VIEW  has  incorporated  this  feature.  If  the 
FORTRAN  compiler  on  the  target  system  does  not  allow  this, 
subroutines  BUILDIE,  CONVREC,  LSEGINT,  OVERLAP,  and  all  refer¬ 
ences  to  these  subroutines  must  be  trimmed  to  six  characters  or 
less. 

b.  This  version  of  VIEW  was  implemented  on  a  32  bit  machine.  The 
character  packing  per  word  factor  was  set  to  4  (A4).  If  the 
target  computer  is  other  than  a  32  bit  machine,  the  character 
packing  factor  will  have  to  be  changed  (16  bit  to  A2,  60  bit  to 
A7,  etc.).  This  is  done  by  altering  the  format  with  statement 
label  200  in  the  mainline  and  the  format  with  statement  label 
200  in  subroutine  TITLE.  Allow  for  8  bits  per  character. 
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c.  Hollerith  strings  In  this  version  of  VIEW  are  enclosed  by  single 

quotes  On  other  target  systems,  new  delimiter 

characters  may  be  needed  (e.g.,  *.  .  .*,  etc.). 

d.  In  READ  statements,  the  END  *  feature  was  used  to  check  for  end- 
of-flle  on  the  Input  files.  If  this  feature  is  not  In  the 
target  system  FORTRAN  compiler,  an  alternative  means  of  checking 
for  end-of-flle  will  have  to  be  used  (e.g.,  EOF  function  on  CDC 
systems ) . 

e.  The  call  to  the  Calcomp  subroutine  PLOTS  in  the  mainline  is  for 
a  Perkin-Elmer  system  [call  PLOTS (0,  0,  LUPLOT)  where  LllPLOT  is 
the  logical  unit  number  of  the  plotting  device].  Consult  your 
system  operator  for  the  protocol  used  on  the  target  system  call 
to  PLOTS. 

f.  The  subroutine  NFRAME  is  designed  for  a  Grinnell  color  graphics 
subsystem  on  a  Perkin-Elmer  3240.  If  this  option  is  chosen  for 
the  VIEW  program,  consult  your  system  operator  for  the  protocol 
needed  to  interface  to  your  graphics  system. 


3.0  VIEW  PROGRAM  MAINLINE  DESCRIPTION 


The  general  functions  of  the  VIEW  program  are  listed  in  the  numbered  blocks 
in  the  flow  chart  of  the  main  routine  (Section  5.1).  A  general  description 
of  each  block  will  be  given.  For  further  information  on  a  particular 
block,  see  the  subroutine  descriptions  later  in  this  report. 

Block  No.  1  calls  subroutine  INPUT  to  read  data  from  the  input  control  file 
and  the  ATBM  input  file.  The  input  control  file  defines  parameters  for 
selecting  what  data  are  to  be  plotted.  The  ATBM  input  file  contains  data 
generated  from  the  ATBM  program  on  unit  1  that  defines  object  sizes,  orien¬ 
tations,  and  positions. 

Block  No.  2  calls  subroutine  CONVREC  to  convert  rectangles  used  as  contact 
planes  in  the  ATBM  program  to  polygons  used  in  the  VIEW  program  and  to 
transfer  planar  coordinates  for  all  polygons  to  the  inertial  reference 
system.  See  the  discussion  of  CONVREC  (Section  4.4)  for  further  informa¬ 
tion. 

Block  Nos.  3,  4,  and  6  are  used  to  decrease  the  computation  time  of  the 
hidden  line  subroutines.  Without  these  three  blocks,  the  hidden  line  sub¬ 
routines  would  have  to  check  out  every  point  or  'vector  head'  against  all 
objects  to  determine  if  that  point  was  hidden  or  not  hidden.  The  subrou¬ 
tines  contained  within  these  three  blocks  project  all  objects  onto  the 
projection  plane  and  generate  perimeters  for  all  the  shadows  of  the 
projected  objects. 

In  Block  No.  3,  subroutine  PRJPLY  projects  the  ATBM  and  input  polygons  onto 
the  projected  plane.  This  is  the  form  needed  by  the  analysis  performed  by 
subroutine  BUILDIE  in  Block  No.  6. 

In  Block  No.  4,  subroutine  PRJELR  performs  the  projection  of  the 
ellipsoids.  Ellipsoids,  in  general,  project  as  complex,  nonelliptical 
shadows.  If  the  viewpoint  coordinate  system  points  directly  at  the  center 
of  the  ellipsoid,  the  ellipsoid  projects  to  an  elliptical  shadow.  The 
subroutine  assumes  that  this  is  the  case.  PRJELR  projects  all  ellipsoids 
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as  elliptical  shadows  circumscribed  with  rectangles  to  conform  with  the 
format  of  subroutine  BUILDIE. 

In  Block  No.  5,  subroutine  POL YD  Is  called  to  generate  the  directional 
cosine  matrices  for  all  polygons.  See  the  description  of  POLYD 
(Section  4.22)  for  further  Information. 

In  Block  No.  6,  subroutine  BUILDIE  takes  the  polygons  formed  by  the  third 
and  fourth  blocks  to  build  an  array  that  defines  object  overlap.  RUILDIE 
starts  with  the  first  object  and  records  the  object  numbers  of  all  polygons 
that,  overlap  with  the  first  object.  Then  the  second  object  is  selected  and 
so  forth  until  all  objects  and  their  overlapping  objects  are  stored  in  the 
array  IE. 

Blocks  No.  7  and  No.  8  call  subroutines  PSE  and  PLPLN  to  check  each  vector 
against  all  the  blocking  objects  recorded  In  the  array  IE  with  the  hidden 
line  subroutines.  If  the  vector  is  not  hidden,  it  is  projected  through  a 
point  onto  the  projection  plane.  The  resultant  vector  is  plotted.  Note 
that  while  projections  performed  by  PRJPLY  and  PRJELR  are  only  on  the 
perimeter  of  the  shadow  of  each  object,  In  PSE  and  PLPLN,  the  projections 
are  of  vectors  that  make  up  the  contours  of  the  objects.  The  program 
returns  to  Block  No.  1  to  get  the  next  data  set. 


4.0  VIEW  PROGRAM  SUBPROGRAM  DESCRIPTIONS 


No. 

Name 

Type 

1 

MAIN 

Main  Program 

2 

BUILDIE 

Subroutine 

3 

CLIP 

Subroutine 

4 

CONVREC 

Subroutine 

5 

CROSS 

Subroutine 

6 

DEI 

Real  Function 

7 

DOT 

Subroutine 

8 

DOTT 

Subroutine 

9 

DRCYPR 

Subroutine 

10 

ELIPSN 

Subroutine 

11 

EXTEND 

Subroutine 

12 

GENDCM 

Subroutine 

13 

HIDE 

Subroutine 

14 

HYDE 

Subroutine 

15 

INPUT 

Subroutine 

16 

LSEGINT 

Subroutine 

17 

MAT 

Subroutine 

18 

NFRAME 

Subroutine 

19 

OVERLAP 

Subroutine 

20 

PLPLN 

Subroutine 

21 

PNTPLT 

Subroutine 

22 

POLYD 

Subroutine 

23 

PREPLT 

Subroutine 

24 

PRJELR 

Subroutine 

25 

PRJPLY 

Subroutine 

26 

PSE 

Subroutine 

27 

ROT 

Subroutine 

28 

SOLVA 

Subroutine 

29 

SOLVR 

Subroutine- 

30 

TITLE 

Subroutine 

31 

TPOINT 

Subroutine 

32 

TRANS1 

Subroutine 

33 

XINTCP 

Real  Function 

! 

i 
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34 

XYZ 

Subroutine 

35 

YINTCP 

Real  Function 

36 

XYZ 

Subroutine 

37 

Z 

Subroutine 

4.1  MAIN  PROGRAM 

a.  Purpose 

Main  program  for  the  VIEW  plotting  package  providing  graphical 
representation  of  the  ATB  model  output.  Controls  the  Initial  1 
zatlon,  data  Input,  data  processing,  and  plotting  output  func¬ 
tions  of  the  program. 

b.  Subroutines  Required 

BUILDIE,  CONVREC,  DOT,  ELIPSN,  INPUT,  MAT,  NEWPEN,  NFRAME, 
NUMBER,  PLOT,  PLOTS,  PLPLN,  POLYD,  PRJELR,  PRJPLY,  PSE,  SYMBOL 
TITLE 

c.  Labeled  Common  Blocks  Used 

ATB,  DBUG,  ELLIPSE,  INTERS,  PLTT,  POLYGON,  VIEWP 

d.  Input  or  Argument  Parameters 
Input  cards  1.0,  3.0,  and  4.0 

e.  Optional  Output 

IDEBUG  (1)  =  1,  print  NIE  array 
IDEBUG  (2)  =2,  print  IE  array 
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f.  Procedure 


1.  Read  input  cards  1.0,  3.0,  and  4.0.  Call  subroutine  INPUT 
to  get  ATB  data. 

2.  Call  subroutine  CONVREC  to  convert  ATB  rectangles  to  VIEW 
polygon  and  to  transform  all  polygon  data  from  the  segment 
coordinate  system  to  the  inertial  coordinate  system. 

3.  Call  subroutine  PRJPLY  to  project  polygons  onto  the 
projected  plane. 

4.  Call  subroutine  PRJELR  to  circumscribe  ellipsoids  with 
rectangles  and  project  them  to  the  projected  plane. 

5.  Call  subroutine  POLYD  to  generate  directional  cosine 
matrices. 

6.  Call  subroutine  BUILDIE  to  define  object  overlap. 

7.  Call  subroutines  PSE  and  PLPLN  to  plot  ellipsoids  and 
polygons,  then  go  get  next  data  set. 

4.2  SUBROUTINE  BUILDIE 
a.  Purpose 

The  BUILDIE  subroutine  is  called  by  the  main  routine  to  build 
the  IE  and  NIE  arrays.  These  arrays  are  used  by  the  hidden  line 
routines  to  determine  if  any  objects  block  the  point  being 
plotted.  This  subroutine  uses  the  data  stored  in  the  CONVEC 
array.  See  Figure  4.1  for  a  pictorial  example  of  the  data  used 
by  BUILDIE  routine.  The  rectangles  represent  the  projected 
ellipsoids  of  the  ATB  program.  The  ellipses  are  circumscribed 
with  rectangles  because  algorithms  dealing  with  polygon  overlap 
require  less  memory  and  computation  time  than  with  algorithms 
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that  look  for  ellipse  overlap.  The  rectangles  shown  In 
Figure  4.1  are  sufficient  to  assure  that  no  object  overlap 
will  be  missed.  The  area  inside  the  rectangles  is  larger  than 
the  area  Inside  the  circumscribed  ellipse.  While  there  is  a 
chance  that  the  8UILDIE  routine  would  find  two  objects  to  over¬ 
lap  that  really  just  missed,  such  a  result  will  not  make  the 
plot  incorrect;  It  just  increases  the  computation  time  of  the 
hidden  line  routines. 

b.  Subroutines  Called 

OVERLAP 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  INTERS,  POLYGON,  REMOVE 

d.  Input  or  Argument  Parameters 
None 

e.  Optional  Output 
None 

f.  Procedure 

The  flow  of  BUILDIE  starts  with  Initializing  the  NIE  and  IE 
arrays.  If  the  number  of  segments  is  equal  to  zero,  the  routine 
goes  on  and  examines  the  polygons.  If  there  are  segments,  the 
first  one  Is  selected,  and  Immediately  Its  own  segment  number  is 
inserted  in  the  IE  array  to  indicate  that  it  can  hide  portions 
of  itself.  Next,  any  other  segment  that  Is  found  to  overlap 
with  the  first  segment  is  recorded  in  the  IE  array.  The  fact  of 
the  overlap  is  recorded  two  places,  once  In  the  portion  of  the 
array  for  segment  No.  1,  and  once  in  the  portion  of  the  array 
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Figure  4.1  Example  of  Data  Used  by  BUILDIE 
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that  represents  the  other  segment.  In  this  way  all  segments  are  checked 
against  each  other  until  all  segment -segment  overlaps  are  recorded  In  the 
IE  array.  After  the  segment -segment  overlaps  are  recorded,  they  are 
recorded  In  the  IE  array.  Segment  numbers  are  always  less  than  31,  and 
polygon  numbers  are  always  greater  than  30. 

4.3  SUBROUTINE  CLIP 

a.  Purpose 

The  purpose  of  this  subroutine  is  to  correctly  execute  the  first 
draw  that  crosses  a  boundary  (X  or  Y)  of  the  valid  plotting 
region.  What  is  meant  by  "correctly  execute"  is  that  the 
intended  draw  extends  to  the  boundary,  but  not  beyond  it.  All 
further  plotting  after  this  draw  is  "clipped,"  i.e.,  not  drawn, 
until  the  pen  returns  to  the  valid  plotting  region.  The  valid 
plotting  region  is  described  by  the  variables  XMIN,  XMAX,  YMIN, 
and  YMAX .  CLIP  also  takes  care  of  the  special  circumstance  that 
the  point  being  clipped  is  the  first  point  in  a  line  segment. 

b.  Subroutines  Called 
PLOT,  XINTCP,  YINTCP 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

X  —  X  Coordinate  of  Point  to  be  Plotted 

Y  —  Y  Coordinate  of  Point  to  be  Plotted 

XSAV  —  X  Coordinate  of  Last  Point  Plotted 

YSAV  —  Y  Coordinate  of  Last  Point  Plotted 

XMIN  —  X  Coordinate  of  Left  Side  of  Plotting  Region 
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XMAX  —  X  Coordinate  of  Right  Side  of  Plotting  Region 
YMIN  —  Y  Coordinate  of  Bottom  of  Plotting  Region 

YMAX  —  Y  Coordinate  of  Top  of  Plotting  Region 

IPEN  —  Calcomp  Pen  Control  Variable 

IPLOT  —  Flag  that  Tells  CLIP  if  Last  Pen  Move 

Was  Clipped  (IPL0T=(J)  or  Not  Clipped  (IPL0T=1) 

e.  Optional  Output 
None 

f.  Procedure 

The  first  function  of  CLIP  is  to  handle  the  special  circumstance 
that  the  point  being  clipped  is  the  origin  of  a  line.  The 
reason  this  merits  special  treatment  is  that  there  is  no  pre¬ 
vious  point  in  the  line  from  which  a  new  point  can  be  inter¬ 

polated.  The  subroutine  must  wait  until  the  next  call  to  do 
this.  So,  the  first  step  in  the  subroutine  is  to  check  the  last 
call  flag  and  see  if  it  is  set  (LCALL  =  1).  If  it  is  set,  it 
means  that  on  the  previous  call,  a  clip  was  performed  on  the 
origin  of  a  line  segment.  What  CLIP  does  is  determine  the  X  and 
Y  coordinates  to  move  the  pen.  The  X  coordinate  is  either  the  X 
intercept  if  the  Y  coordinate  was  off  the  plotter  or  the  X 
boundary  value  if  the  Y  was  on  the  plotter.  The  Y  coordinate  is 
either  the  Y  intercept  value  if  the  X  coordinate  was  off  the 
plotter  or  the  Y  boundary  value  if  the  X  coordinate  was  on  the 
plotter.  CLIP  moves  the  pen  (up,  IPEN  =  3)  to  the  point  denoted 
by  the  X  and  Y  coordinates,  resets  the  previous  X  and  Y  coordi¬ 
nate  variables  (XSAV  and  YSAV),  and  clears  the  last  call  flag 
(LCALL  =  0).  If  LCALL  was  not  set  at  the  beginning  of  the 
subroutine,  CLIP  skips  this  portion  of  the  code. 

The  second  step  of  this  first  function  is  to  check  if  the 
current  point  being  clipped  is  the  origin  of  a  line.  If  this  is 
true  (IPLOT  =  1  and  IPEN  =3),  the  coordinates  of  the  point  are 
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saved  In  XLSAV  and  YLSAV,  and  the  last  call  flag  Is  set  to  one 
(LCALL  ■  1).  If  this  second  step  was  executed,  the  subroutine 
exits  here.  If  not,  It  continues. 

The  second  function  of  CLIP  starts  here.  The  subroutine 
determines  If  the  clip  flag  has  been  set  (IPLOT  »  1).  This  is 
the  flag  that  enables  only  the  valid  portion  of  the  first  seg¬ 
ment  clipped  to  be  plotted.  IPLOT  Is  cleared  after  plotting  and 
Is  kept  clear  until  the  pen  returns  to  a  valid  X  or  Y  coordi¬ 
nate.  It  Is  then  reset  by  a  draw  In  the  valid  plotting 
region.  The  X  and  Y  coordinates  are  calculated  by  looking  at 
the  XOFF  and  YOFF  flags  and  using  either  XINTCP  and  YINTCP  func¬ 
tions  or  the  X  and  Y  boundaries.  A  line  Is  drawn  to  that 
point.  The  clip  flag  IPLOT  is  cleared  and  the  subroutine 
exits. 

4.4  SUBROUTINE  CONVREC 
a.  Purpose 

The  purpose  of  subroutine  CONVREC  is  as  follows.  The  ATB 
simulation  program  outputs  plane  information  in  the  form  of 
three  plane  equations  that  define  boundary  planes  that  bound  the 
rectangle.  The  coordinate  system  used  for  defining  these  planes 
can  be  any  one  of  the  following  three. 

1.  A  rectangle  defined  in  the  inertial  reference  frame. 

2.  A  rectangle  defined  in  the  vehicle  reference  frame. 

3.  A  rectangle  defined  in  the  segment  reference  frame. 

The  VIEW  plotting  package  works  with  polygons  in  an  Inertial 
reference  frame.  Therefore,  each  of  the  three  possibilities 
described  above  must  be  converted  to  the  format  used  Internally 
In  the  VIEW  plotting  package. 
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Subroutine  CONVREC  converts  rectangles  In  the  ATB  simulation 
format  to  polygons  In  the  VIEW  plotting  format.  It  then  trans¬ 
forms  the  vectors  to  the  vertices  of  all  the  polygons  from  the 
coordinate  system  of  the  segment  they  are  tied  to  the  Inertial 
reference  system. 

b.  Subroutines  Called 
DET,  DOT 

c.  Labeled  Common  Blocks  Used 

ATB,  CONECT,  DBUG,  ELLIPSE,  POLYGON 

d.  Input  or  Argument  Parameters 
None 

e.  Optional  Output 

I  DEBUG (4 )  =  1;  print  reference  segment  number  (ISG),  plane 
number,  and  plane  vectors  in  inertial  reference. 

f.  Procedure 

The  equations  that  represent  the  rectangle  in  the  ATB  output  are 
as  follows: 

1.  AOX  +  BOY  +  COZ  =  DO 

2.  A1X  +  B1Y  +  C1Z  =  D1 

3.  A2X  +  B2Y  +  C2Z  =  D2 

Solving  these  equations  simultaneously  gives  a  point  of 
intersection  of  all  three  planes  (see  Figure  4.2). 
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This  point  Is  the  point  PI  used  by  the  VIEW  program  for  defining 
polygons.  The  VIEW  program  defines  polygons  by  specifying  the 
position  of  all  corners,  r  Is  the  normal  to  plane  number  3,  and 
r3  Is  the  normal  to  the  plane  number  2.  The  length  of  r2  and 
?3  are  obtained  from  the  output  of  the  ATB  simulation  program. 
Each  corner  of  the  polygon  Is  found  by  adding  the  contour  vec¬ 
tors  of  the  rectangle  In  succession.  This  Is  done  for  all  poly¬ 
gons  read  from  the  ATB  input  file. 

REFERENCE  FRAME 


Figure  4.2  Point  of  Intersection  for  all  Three  Planes  in  CONVREC 
4.5  SUBROUTINE  CROSS 

a.  Purpose 

Computes  vector  cross  product  £  =  &  * 

b.  Subroutines  Required 
None 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  Parameters 

A,  B,  C:  Arrays,  each  consisting  of  three  elements,  that 
represent  vectors  where  the  cross  product  Is  defined  by 

£  =  A  *  B. 

e.  Optional  Output 
None 

f.  Procedure 

Computes  each  element  of  £  by 
ci  =  a 2b 3  -  a 3b 2 
c2  =  a  3b 1  -  a ib 3 

C3  =  a  ib 2  -  a2bi 
4.6  FUNCTION  DET 

a.  Purpose 

DET  finds  the  determinant  of  the  3x3  square  array  passed  to 
it. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  Parameters 

Values  Representing  Square,  3x3  Array 

e.  Optional  Output 
None 

f.  Procedure 

The  determinant  is  calculated  in  the  following  manner: 

DET  =  A11*(A22*A33-A23*A32)  -  A12*(A21*A33-A23*A31)  + 
A13*(A21*A32-A22*A31 ) 

.7  SUBROUTINE  DOT 

a.  DOT  performs  matrix  multiplication  C  *  If  A  and  &  are 

vectors,  C  is  the  dot  product  ft  • 

b.  Subroutines  Required 
None 

c.  Labeled  Common  Blocks  Used 
None 


All 

A12 

A13 

A21 

A22 
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A31 

A32 

A33 
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d.  Input  or  Argument  parameters 

A  -  Matrix  of  Size  (L,N) 

B  -  Matrix  of  Size  (L,M) 

C  -  Product  Matrix  of  Size  (N,M) 

N,M,L  -  Sizes  of  Matrices  A.B.C 

e.  Optional  Output 
None 

f.  Procedure 

Each  element  C(I,J)  of  the  product  matrix  C  is  computed  by: 
L 

C(I,J)  =  I  A(K.I)  *  B(K,J)  for  I  =  1,N  and  J  =  1,M 
K  =  1 

4.8  SUBROUTINE  DOTT 

a.  Purpose 

DOTT  performs  the  matrix  multiply  C_  =  AB1 

b.  Subroutine  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  Parameters 

A  -  Matrix  of  Size  (N,L) 

B  -  Matrix  of  Size  (M,L) 

C  -  Matrix  of  Size  (N,M) 

N,M,L  -  Sizes  of  Matrices  A#B,C 

e.  Optional  Output 
None 

f.  Procedure 

Each  element  C(I,J)  of  the  product  matrix  C  is  computed  as: 

L 

C(I,J)  =  Z  A(I ,K)  *  B(J,K)  for  I  =  1,N  and  J  =  1,M 
K  =  1 

4.9  SUBROUTINE  DRCYPR 

a.  Purpose 

Sets  up  direction  cosine  matrix  £  for  rotation  angles  A  given  in 
degrees  about  local  the  x,  y,  or  z  axis  of  the  segment  In  ques¬ 
tion  as  indicated  by  II,  12,  or  13. 

b.  Subroutines  Required 
MAT,  ROT 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  Parameters 


D  —  3x3  direction  cosine  matrix  to  be  computed. 

A  —  3  rotation  angles  given  In  degrees. 

ID  —  3  Integers  (II,  12,  and  13)  that  Indicate  axis  of 

rotation  for  each  of  the  three  angles  In  A  (1,  2,  or  3 
Indicates  x,  y,  or  z  axis,  respectl vely). 

e.  Optional  Output 
None 

f.  Procedure 

Computes  as  a  matrix  product 


H  = 


where  each  Dj  is  one  of  the  following 
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Note: 


(1)  For  the  normal  sequence  of  yaw  (f),  pitch  (e),  and  roll 
(4),  let  ID  a  3,  2,  1  to  obtain 

2  “  2$  2a  2*  “  £i  h  £  a 

(2)  For  the  reverse  sequence  as  required  by  Subroutine  INITIAL 
prior  to  Version  18  of  the  ATB  program,  let  ID  *  1,  2,  3 
to  obtain 

£  =  £4  £e  £4  =  £3  22  £i 

(3)  For  Euler  angles,  precession  (4),  nutations  (0),  and  spin 

(4),  let  ID  =  3,  1,  -3  to  obtain 

2  =  24  2e  £4  =  23  21  23 


4.10  SUBROUTINE  ELIPSN 

a.  Purpose 

Subroutine  ELIPSN  generates  the  XI  array  of  contour  vectors  for 
a  quarter  of  the  ellipsoid.  These  contours  are  used  by  sub¬ 
routine  PSE  to  generate  complete  contours  for  the  entire  ellip¬ 
soid. 

b.  Subroutines  Called 
SQRT 

c.  Labeled  Common  Blocks  Used 

ELLIPSE 
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d.  Input  of  Argument  Parameters 

INDEX  --  Number  of  Steps  Plus  One 
XI  --  Array  Containing  a  Quarter  of  the  Ellipsoid 

IN  --  Array  Containing  Number  of  Points  in  each 

Contour  Array 

e.  Optional  Output 
None 

f.  Procedure 

ELIPSN  is  called  once  for  each  ellipsoid.  For  each  step  for 
ellipsoid,  calculate  X  and  Y  coordinates.  Calculate  test 
variable  [TEST  =  (1-X2*SIMP1)  -  Y2*SIMP2].  If  result  is  nega¬ 
tive,  it  means  point  is  not  on  the  ellipsoid.  In  this  case,  Y 
coordinate  is  calculated  as 

/  1-X2  *  SIMP1 

V  - 5T  W2 - 

and  Z  is  set  to  zero.  The  results  are  stored  in  XI.  If  TEST 
was  positive,  Z  is  calculated  and  the  results  are  stored  in 
XI. 

4.11  SUBROUTINE  EXTEND 
a.  Purpose 

EXTEND  is  a  subroutine  used  in  conjunction  with  the  hidden  line 
routines.  This  routine  removes  large  gaps  that  could  exist 
around  boundaries  where  contour  lines  are. hidden  and  not  hid¬ 
den.  The  gaps  are  caused  by  dividing  up  contour  lines  into  a 
set  of  vectors.  Once  it  is  determined  that  a  particular  vector 
crosses  a  boundary,  only  a  portion  of  the  vector  must  be 
plotted.  If  the  entire  vector  were  left  unplotted,  a  gap  would 
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exist  around  these  boundaries.  The  size  of  the  gaps  would  then 
cover  the  range  from  very  small  to  the  size  of  the  vector.  The 
EXTEND  subroutine  Is  called  whenever  a  vector  passes  through  one 
of  these  boundaries.  EXTEND  finds  the  portion  of  the  vector  not 
hidden  and  returns  this  Information  to  the  plotting  subrou¬ 
tine. 

b.  Subroutines  Called 
HIDE,  HYDE 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  INTERS,  PLTT 

d.  Input  or  Argument  Parameters 

P  —  Contains  X,  Y,  and  Z  coordinates  of  the  two  end 
points  of  the  line. 

I  —  Points  to  unhidden  point  location  in  P  array. 

J  --  Points  to  hidden  point  location  In  P  array. 

e.  Optional  Output 
None 

f.  Procedure 

EXTEND  starts  by  finding  the  midpoint  for  a  line  between  the  X, 
Y,  and  Z  coordinates  found  In  array  P.  This  midpoint  is  passed 
to  either  subroutine  HYDE  (ellipsoids)  or  HIDE  (polygons)  to 
check  If  contour  line  Is  hidden.  If  It  Is  hidden,  the  midpoint 
coordinates  are  loaded  Into  row  J  of  array  P.  If  it  is  not 
hidden,  the  midpoint  Is  loaded  Into  row  I  of  array  P. 
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4.12  SUBROUTINE  GENDCM 


a.  Purpose 

The  name  of  GENDCM  stands  for  generate  direction  cosine 
matrix.  This  subroutine  creates  a  direction  cosine  matrix  from 
Information  that  is  readily  available  to  the  user  of  program 
VIEW.  GENDCM  requires  two  position  vectors  to  define  the  view¬ 
point  position  and  orientation.  Vector  V?  defines  the  position 
of  the  viewpoint  in  the  reference  frame  of  the  segment  to  which 
the  viewpoint  is  attached.  Vector  RA  defines  a  point  where  the 
viewpoint  coordinate  system  is  aimed.  The  aiming  of  the  view¬ 
point  coordinate  system  is  determined  by  the  direction  of  the 
positive  Z  axis.  The  viewpoint  always  looks  down  the  positive  Z 
axis. 

b.  Subroutines  Called 
SQRT 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

CAMERA  —  Position  vector  for  the  viewpoint  in  the 

reference  frame  of  the  segment  to  which  the 
viewpoint  is  attached  (X,  Y,  Z  coordinates). 

FOCUS  —  Position  vector  for  point  which  viewpoint  Z  axis 
is  aimed  (X,  Y,  Z  coordinates  in  the  reference 
frame  of  the  segment  to  which  the  viewpoint  is 
attached). 

D  --  Direction  cosine  matrix  for  viewpoint 

(transformed  from  the  inertial  coordinate  system 
to  the  viewpoint  coordinate  system). 
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e.  Optional  Output 


None 


f.  Procedure 


GENDCM  assumes  the  X  axis  of  the  viewpoint  coordinate  system 
remains  parallel  to  the  X-Y  plane  of  the  inertial  reference 

A  A 

frame.  Also,  X  axis  must  be  normal  to  the  Z  axis  In  the  X-Y 

A 

plane  of  the  Inertial  reference  frame.  Z  In  the  X-Y  plane  Is 
given  by 

Zl  _Z£ 

W  FT’ 

and 


A  A  A  A 

must  be  (Z2,  -  Zi,  0).  Now  both  the  X  and  Z  vectors  of  the 

A 

viewpoint  coordinate  system  are  determined.  The  Y  is  obtained 

A  A 

by  crossing  Z  into  X.  The  direction  cosine  matrix  has  the 
following  form. 


Z ( 2 ) /XNORM 

-Z(l) /XNORM 

0.0 

Z(1)*Z(3)/XN0RM 

Z(2)*Z(3) /XNORM 

-XNORM 

Z(l) 

1(2) 

Z(3) 

where  XNORM  is  |ft|,  the  length  of  it. 
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4.13  SUBROUTINE  HIDE 


a.  Purpose 

The  purpose  of  HIDE  is  to  determine  if  a  point  is  being  hidden 
by  another  polygon.  The  subroutine  corresponds  with  subroutine 
HYDE  (HYDE  checks  ellipsoids). 

b.  Subroutines  Called 

DOT,  MAT,  TPOINT,  TRANS1 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  POLYGON 

d.  Input  or  Argument  Parameters 

KK  --  Polygon  number  to  check  blocking. 

P3  —  Contains  position  vector  of  point  to  check. 

IFLAG  --  Flag  passed  back  to  caller  indicating  hidden 
(IFLAG=1)  or  not  hidden  (IFLAG=2). 

e.  Optional  Output 
None 

f.  Procedure 

The  subroutine  first  determines  if  the  projected  point  and  the 
projected  polygon  overlap.  If  the  projected  point  is  outside  of 
the  projected  polygon,  the  point  cannot  be  hidden  by  the  polygon 
being  examined.  If  the  projected  point  does  lie  within  the 
projected  polygon,  it  must  be  determined  which  Is  closer  to  the 
viewpoint.  If  the  point  Is  closer  than  the  polygon,  that  point 
is  not  hidden  by  the  polygon  being  examined.  Appendix  C, 
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"Intersection  of  a  Three  Space  Vector  and  Plane,"  describes  In 
detail  how  to  determine  whether  the  plane  or  point  Is  closer  to 
the  viewpoint. 

The  overlap  on  the  projection  plane  of  a  projected  point  and  a 
projected  polygon  Is  determined  by  subroutine  TPOINT.  TPOINT 
returns  a  flag  called  IFLAG  that  Indicates  overlap  or  no  over¬ 
lap. 

4.14  SUBROUTINE  HYDE 

a.  Purpose 

The  purpose  of  HYDE  Is  to  determine  If  a  point  Is  hidden  by  an 
ellipsoid.  If  the  point  Is  hidden,  variable  IFLAG  Is  returned 
as  one;  If  not  hidden,  IFLAG  Is  set  to  two. 

b.  Subroutines  Called 

ABS,  SQRT,  DOT,  DOTT,  MAT,  XYZ,  YZ,  Z 

c.  Labeled  Common  Blocks  Used 

ELLIPSE 

d.  Input  or  Argument  Parameters 

N  —  Possible  hiding  ellipsoid  number  (I) 

R  —  Vector  to  plotting  point  (I) 

IFLAG  —  Flag  Indicating  hidden  (IFLAG*1)  or  not 
(IFLAG»2)  (0) 

e.  Optional  Output 
None 
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f.  Procedure 


The  equations  used  by  subroutine  HYDE  (and  accompanying 
subroutines  XYZ,  YZ,  and  Z)  can  be  found  in  Appendix  A,  "Hidden 
Line  Problem  Between  Two  Ellipsoids."  Refer  to  this  appendix 
for  further  Information  on  HYDE. 

4.15  SUBROUTINE  INPUT 

a.  Purpose 

INPUT  has  two  separate  functions.  The  first  is  to  read  initial 
data  from  the  input  control  file,  read  polygon  data  from  the  ATB 
input  file,  and  to  initialize  variables.  This  is  all  done  on 
the  first  call.  The  second  function,  performed  on  all  sub¬ 
sequent  calls,  is  to  read  ellipsoid  data  from  the  ATB  input 
file. 

b.  Subroutines  Called 
DOT,  DRCYPR,  GENDCM 

c.  Labeled  Common  Blocks  Used 

ATB,  CONECT,  DBUG,  ELLIPSE,  INTERS,  PLTT,  POLYGON,  REMOVE, 

VIEWP 

d.  Input  or  Argument  Parameters 

CTIME  —  Current  time  of  program  passed  from  main. 

Input  Cards  —  6.0,  6.1,  7.0,  7.1,  7.2,  8,0,  9.0  10.0, 

11.0,  12.0,  13.0,  14.0,  15.0,  15.1 
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e.  Optional  Output 


IDEBUG (3 )  =  1;  print  viewpoint  position  vector,  viewpoint 

direction  cosine  matrix,  NSTEPS  from  cards  7.0 
and  8.0,  segment  Inertial  position  vector, 
segment  Inertial  direction  cosine  matrix,  NSEG, 
NPL,  SFACTR,  OFSETX,  and  OFSETY. 

f.  Procedure 

During  the  first  call,  cards  6.0,  6.1,  7.0,  7.1,  7.2,  8.0,  9.0, 
10.0,  11.0,  12.0,  13.0,  14.0,  15.0,  and  15.1  are  read  from  the 
input  control  file.  Polygon  data  are  read  from  the  ATB  input 
file.  The  direction  cosine  matrix  is  initialized  in  three 
different  ways,  depending  on  the  input  flag  ICODE.  When  ICODE 
equals  zero,  the  roll,  pitch,  and  yaw  angles  are  found  in  the  RA 
array  (card  15.0);  and  subroutine  DRCYPR  calculates  the  direc¬ 
tion  cosine  matrix  When  ICODE  equals  one,  the  direction  cosine 
matrix  is  supplied  as  input  (card  15.1).  When  ICODE  is  equal  to 
two,  subroutine  GENDCM  generates  the  direction  cosine  matrix. 

The  subroutine  exits. 

On  all  subsequent  calls,  INPUT  starts  at  FORTRAN  statement 
#600.  Segment  data  are  read  from  the  ATB  input  file.  This 
reading  is  continued  until  the  time  of  the  data  is  greater  than 
or  equal  to  the  current  time  of  the  program  (CTIME).  When  this 
time  is  reached,  the  current  location  of  the  contact  ellipsoid 
is  calculated.  If  the  time  interval  of  the  ATBM  simulation  data 
is  greater  than  the  DT  of  the  view  run,  variable  IFLAG  is  set  to 
ten.  This  is  done  to  signal  the  program  to  continue  Incre¬ 
menting  CTIME  without  making  plots  until  CTIME  reaches  the  next 
available  time  point  in  the  simulation  data.  The  subroutine 
then  exits. 
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4.16  SUBROUTINE  LSEGINT 


a.  Purpose 

LSEGINT  Is  cal  led  by  subroutine  OVERLAP  to  check  for  two  lines 
Intersecting.  The  line  segments  are  represented  by  two  end 
points  for  each  line. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

PI  --  X  and  Y  coordinates  for  line  No.  1,  end  point 

No.  1 

P2  --  X  and  Y  coordinates  for  line  No.  1,  end  point 

No.  2 

R1  --  X  and  Y  coordinates  for  line  No.  2,  end  point 

No.  1 

R2  —  X  and  Y  coordinates  for  line  No.  2,  end  point 

No.  2 

IFLag  —  Flag  passed  back  to  OVERLAP  to  show  intersection 
(IFLAG=1)  or  no  Intersection  (IFLAG=0) 

e.  Optional  Output 
None 
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f.  Procedure 


LSE3INT  checks  for  six  separate  cases: 


Case  1 
Case  2 
Case  3 
Case  4 
Case  5 
Case  6 


Regular  configuration 

One  line  Is  vertical 

Both  lines  are  vertical 

Both  lines  are  horizontal 

Both  lines  have  the  same  nonzero  slope 

One  line  is  vertical,  the  other  Is  horizontal 


It  first  checks  for  Cases  3  and  4.  If  either  of  these  cases  are 
true,  the  subroutine  returns.  For  Cases  1,  2,  5,  and  6,  refer 
to  Appendix  0,  "Intersection  of  Line  Segments  In  a  Plane."  The 
code  of  subroutine  LSEGINT  follows  the  equations  and  text  found 
in  the  appendix. 


4.17  SUBROUTINE  HAT 

a.  Purpose 

Performs  the  matrix  multiple  C  c  AB. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

A  —  Matrix  of  Size  (LL.MM) 

B  —  Matrix  of  Size  (MM,NN) 

C  --  Product  Matrix  of  Size  (LL,NN) 
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LL,  MM,  NN  —  Sizes  of  Matrices  A,  B,  C 

JA,  JB,  JK  —  First  Dimension  of  A,  B,  C  In  Calling 

Subroutine 

e.  Optional  Output 
None 

f.  Procedure 

Each  element  (C(I,J)  of  the  product  matrix  C  is  computed  by: 

MM 

=  Z  A(I,K)  *  B(K,J)  for  I  =  1,  LL  and  J  =  1,  NN 
K  =  1 

4.18  SUBROUTINE  NFRAME 

a.  Purpose 

NFRAME  performs  all  end  of  frame  operations  necessary  for  a 
Grinnell  graphics  system.  The  code  furnished  in  this  version  of 
the  VIEW  program  is  for  a  Perkin-Elmer  3242  running  under 
OS/32. 

b.  Subroutines  Called 
DOLWH,  PLOTS 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 
None 
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e.  Optional  Output 
None 

f.  Procedure 

Again,  this  subroutine  Is  designed  for  a  Perkln-Elmer  3242 
computer  and  a  particular  graphics  system,  a  Grinnell  GMR-27. 

If  a  different  graphics  system  Is  to  be  utilized,  this  subrou¬ 
tine  will  need  to  be  changed.  Immediately  upon  entering  NFRAME, 
subroutine  DOLWH  is  called.  DOLWH  stands  for  digital  output 
with  handshake.  An  end-of-frame  halfword  (FFFF)  Is  output  to 
the  Grinnel  system.  Subroutine  PLOTS  is  called  to  initialize 
the  next  frame.  The  subroutine  then  exits. 

4.19  SUBROUTINE  OVERLAP 

a.  Purpose 

This  subroutine  is  called  by  8UILDIE  to  check  two  objects  for 
overlap.  The  objects  are  defined  by  the  first  two  arguments  in 
the  call  to  OVERLAP.  At  the  point  where  OVERLAP  is  called,  all 
objects  have  a  projected  polygon  representation,  and  OVERLAP 
determines  if  these  two  polygons  overlap. 

b.  Subroutines  Called 
LSEGINT,  TPOINT 

c.  Labeled  Common  Blocks  Used 

POLYGON 
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d.  Input  or  Argument  Parameters 

III  --  Object  No.  1  to  be  Tested 

KKK  --  Object  No.  2  to  be  Tested 

MFLAG  --  Flag  indicating  overlap  or  not.  MFLAG  =  0, 

no  overlap;  MFLAG  =  1,  overlap. 

e.  Optional  Output 
None 

f.  Procedure 

There  are  two  basic  checks  employed  by  this  subroutine  to 
determine  object  overlap.  First  is  the  POINT  INSIDE  A  POLYGON 
TEST.  This  test  starts  with  a  point,  usually  a  corner  of  a 
polygon,  and  tests  are  made  to  see  if  that  point  lies  inside  or 
outside  the  polygon.  Once  a  point  is  found  to  be  within  the 
other  polygon,  MFLAG  is  set  equal  to  1  to  indicate  object  over¬ 
lap,  and  OVERLAP  returns  to  the  calling  program.  If  no  points 
lie  inside  the  polygons  being  tested,  the  INTERSECTING  LINE 
SEGMENT  TEST  is  used.  This  test  checks  for  line  segments  of  one 
polygon  intersecting  line  segments  of  the  other  polygon.  Again, 
if  intersection  is  found,  MFLAG  is  set  equal  to  1,  and  OVERLAP 
returns  to  BUILDIE.  If  no  intersection  is  found,  the  two 
objects  must  not  overlap,  and  MFLAG  is  set  equal  to  0  and  OVER¬ 
LAP  returns  to  BUILDIE. 

4.20  SUBROUTINE  PLPLN 

a.  Purpose 

This  subroutine  is  called  by  the  Main  program  to  set  up  arrays 
containing  polygon  data  for  subroutine  PNTPLT  to  plot. 
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b.  Subroutines  Called 
NEWPEN,  PNTPLT 

c.  Labeled  Common  Blocks  Used 

DBUG,  ELLIPSE,  PLTT,  POLYGON,  REMOVE 

d.  Input  or  Argument  Parameters 

SEG  --  Array  containing  vectors  representing  sides  of 

the  polygons. 

INDEX2  —  Array  size  for  SEG. 

e.  Optional  Output 
None 

f.  Procedure 

Array  SEG  must  have  the  vectors  that  represent  the  sides  of  the 
polygons.  Each  side  of  the  polygon  Is  represented  by  a  series 
of  short  vectors,  even  though  one  long  vector  could  be 
plotted.  A  series  of  short  vectors  Is  used  so  current  algo¬ 
rithms  can  be  used  for  the  hidden  line  segment  problem.  The 
vectors  in  array  SEG  are  In  the  Inertial  reference  frame.  After 
loading  data  Into  SEG,  It  Is  sent  to  subroutine  PNTPLT  to 
plot.  After  doing  this  for  all  planes,  PLPLN  exits. 

4.21  SUBROUTINE  PNTPLT 

a.  Purpose 

The  purpose  of  PNTPLT  Is  to  plot  all  the  points  passed  to  It  In 
array  SEG.  PNTPLT  Is  called  by  subroutines  PSE  to  plot 


ellipsoids  and  by  PLPLN  to  plot  polygons.  PNTPLT  Is  where  all 
calls  to  the  Calcomp  plotting  subroutines  are  made. 

b.  Subroutines  Called 

CLIP,  EXTEND,  HIDE,  HYDE,  PLOT,  PREPLT,  TRANS1 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  INTERS,  PLTT 

d.  Input  or  Argument  Parameters 

SEG  —  Array  containing  points  to  be  plotted. 

IPEN  —  Calcomp  pen  control  variable. 

INDEX2  --  Maximum  number  of  points  to  plot. 

NPTS  —  Number  of  points  to  plot. 

CARD  16.0  --  XMIN,  XMAX  from  input  control  file. 

e.  Optional  Output 
None 

f.  Procedure 

PNTPLT  takes  each  point  and  checks  to  see  if  the  point  is 
hidden.  This  is  accomplished  by  calling  subroutine  HYDE  to 
check  for  possible  blocking  ellipsoids  and  by  calling  subroutine 
HIDE  to  check  for  possible  blocking  polygons.  After  it  is 
determined  that  the  point  is  or  is  not  hidden,  PNTPLT  checks  to 
see  if  a  boundary  is  crossed  by  the  vector  originating  on  the 
last  point  and  ending  on  the  current  point.  *A  boundary  is  indi¬ 
cated  whenever  the  state  of  the  last  point  and  the  current  point 
is  different.  If  the  last  point  was  not  hidden  and  the  current 
point  is  hidden,  the  vector  goes  from  a  not  hidden  zone  into  a 
hidden  zone.  The  boundary  between  the  two  zones  is  the  boundary 
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that  Is  checked  for  at  this  stage  of  subroutine  PNTPLT.  If  a 
boundary  Is  present,  subroutine  EXTEND  Is  called  to  extend  a  new 
vector  from  the  not  hidden  point  up  to  the  boundary.  This  new 
vector  Is  In  the  same  direction  as  the  vector  that  crosses  the 
boundary. 

Once  a  nonhidden  vector  Is  established,  subroutine  TRANS1  Is 
called  to  transform  the  vector  Into  the  viewpoint  reference 
frame.  Subroutine  PNTPLT  projects  the  vector  that  Is  now  in  the 
viewpoint  system.  This  projection  Is  similar  to  the  projection 
of  a  lens  onto  a  projection  plane  and  Is  represented  by  the 
following  equations: 

v,  _  (SFACTOR)  *  (X) 
a  -  - - 

yl  =  (SFACTOR)  ♦  (Y) 

where  X'  and  Y'  are  the  plotting  coordinates  on  the  projection 
plane. 

X,  Y,  Z  are  the  coordinates  of  the  position  vector  in  the 
viewpoint  coordinate  system. 

SFACTOR  is  the  scale  factor  for  the  plot. 

Once  the  plotting  coordinates  of  the  vector  have  been 
calculated,  the  X  and  Y  coordinates  are  checked  to  see  if  the 
plot  move  will  take  the  pen  off  the  boundaries  of  the  plotter. 
The  X  coordinate  is  compared  to  the  variables  XMIN  and  XMAX  and 
the  Y  coordinate  Is  compared  to  YMIN  and  YMAX.  These  variables 
represent  the  bottom,  top,  left,  and  right  boundary  values  of 
the  plotting  region  and  are  defined  In  subroutine.  If  X  *  Y 
coordinate  Is  off  the  plotter,  subroutine  CLIP  is  called  to 
compensate  for  what  the  pen  was  to  have  done.  The  X  and  Y 
coordinates  of  this  intended  pen  move  are  saved  and  the  pen 
variable  Is  set  to  lift  (IPEN  *  3)  to  prevent  any  drawing  until 
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the  pen  returns  to  the  plotting  region.  If  the  X  and  Y  coordi 
nates  were  In  the  plotting  region,  subroutine  PREPLT  is 
called.  PREPLT  completes  any  pen  moves  made  during  clipping. 
The  Calcomp  subroutine  PLOT  Is  then  called  to  perform  the  pen 
move  or  draw. 

4.22  SUBROUTINE  POLYO 

a.  Purpose 

The  purpose  of  POLYD  is  to  generate  direction  cosine  matrices 
for  all  polygons. 

b.  Subroutines  Called 
CROSS,  SQRT 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  POLYGON 

d.  Input  or  Argument  Parameters 
None 

e.  Optional  Output 
None 

f.  Procedure 

POLYD  generates  a  direction  cosine  matrix  for  all  polygons  in 
View.  The  approach  taken  here  is  simply  to  cross  two  adjacent 
sides  of  the  polygon  to  obtain  the  normal  to  the  polygon  sur¬ 
face— this  is  one  coordinate  vector.  One  of  the  sides  used  in 
the  cross  product  is  picked  as  another  coordinate  vector.  The 
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third  coordinate  vector  Is  obtained  by  crossing  the  first  two. 
A  direction  cosine  matrix  Is  obtained  by  placing  these  three 
vectors  In  a  matrix.  (For  the  following  equations,  see 
Figure  4.3.) 


n  =  P2  -  Pi  x  ?3  -  pi 


~  ajX  +  b i Y  +  CiZ 


-  PT^I 

r2  = -  =  a2X  +  b2Y  +  c2Z 

I^tti 

AAA 

r3  =  ri  x  r2  =  a3X  +  b3Y  +  C3Z 
The  direction  cosine  matrix  is  given  by 


A  A  A  A  A  A 

x  •  n  y  •  n  z  •  ri 

A  A  A  A  A  A 

D  =  X  •  r2  Y  •  r2  Z  •  r2 


A  A  A  A  A  A 

X  •  T3  Y  •  T3  Z  •  r3 
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Figure  4.3  Pictorial  Representation  of  Equations  in  POL  YD 
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and  Y  •  Z l  •  0 

ai 

bi 

Cl 
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a2 

b2 

C2 

i  ' 

a3 

b3 

C3 

i 

• 
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This  direction  cosine  matrix  can  be  used  to  transform  from  the 
inertial  reference  system  to  the  local  reference  system. 


4.23  SUBROUTINE  PREPLT 

a.  Purpose 


The  purpose  of  this  subroutine  is  to  position  the  pen  before  a 
call  to  PLOT  with  pen  down  if  that  move  would  cause  the  pen  to 
exceed  the  specified  plotting  region.  The  correct  position  for 
the  pen  is  at  the  point  denoted  by  the  coordinates  of  an  X  or  Y 
intercept  between  saved  and  present  points  at  the  X  or  Y  bound¬ 
ary,  and  the  X  or  Y  boundary  (Figure  4.4). 

b.  Subroutines  Called 
PLOT,  XINTCP,  YINTCP 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

X  —  X  coordinate  of  present  point 

Y  —  Y  coordinate  of  present  point 


47 


I 


I 

I 

! 

i 

> 

•r 


I 


f 


I 

( 

< 

J 


t 

I 

! 

t 

.  ► 

' 


^  SAVED  X,  Y  COORDS. 
\ 

\ 

\ 

\ 


Figure  4.4 


Figure  for  Variables  in  PREPLT 
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XSAV 
YSAV 
XMIN 
XMAX 
YMIN 
YMAX 
IPEN 
NEWPEN 

e.  Optional  Output 
None 

f.  Procedure 

The  method  in  which  PREPLT  knows  If  the  previous  draw  was  across 
the  border  is  a  combination  of  two  reasons.  The  first  is  the 
fact  that  the  subroutine  is  being  called.  That  is,  the  program 
flow  of  subroutine  PNTPLT  (the  only  subroutine  that  calls 
PREPLT)  will  only  call  PREPLT  if  the  present  X  and  Y  coordinates 
are  on  the  plotter.  The  second  reason  is  determined  by  checking 
flag  NEWPEN.  NEWPEN  is  a  variable  that  has  saved  the  previous 
pen  move.  It  is  also  set  to  a  negative  value  if  the  pen  move 
was  out  of  the  plotting  region.  So,  if  NEWPEN  is  equal  to  -2 
(negative  meaning  clipped,  2  meaning  pen  down),  then  we  want  to 
continue.  If  it  is  not  equal  to  -2,  the  subroutine  exits.  If 
continued,  PREPLT  determines  to  which  coordinates  to  move  the 
pen.  If  the  X  coordinate  is  outside  the  valid  plotting  region, 
the  Y  coordinate  Is  set  to  the  Y  Intercept  between  the  present 
and  saved  coordinates  at  the  X  boundary.  If  the  X  coordinate  is 
inside  the  valid  plotting  region,  the  Y  coordinate  is  set  to  the 
Y  boundary  value.  If  the  Y  coordinate  Is  outside  the  valid 
plotting  region,  the  X  coordinate  is  set  to  the  X  intercept 
value  between  the  present  and  saved  points  at  the  Y  boundary. 

If  the  Y  coordinate  was  Inside  the  valid  plotting  region,  the  X 
coordinate  is  set  to  the  X  boundary  value.  PREPLT  moves  the  pen 


X  coordinate  of  previous  point 

Y  coordinate  of  previous  point 

X  coordinate  of  left  side  of  plotting  region 
X  coordinate  of  right  side  of  plotting  region 

Y  coordinate  of  bottom  of  plotting  region 

Y  coordinate  of  top  of  plotting  region 
Calcomp  pen  control  value 

Saved  value  of  previous  IPEN  move 
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(up)  to  this  spot  denoted  by  the  just  calculated  X  and  Y  coordi¬ 
nates.  Variable  IPEN  Is  set  to  2  and  the  subroutine  exits. 

4.24  SUBROUTINE  PRJELR 

a.  Purpose 

This  subroutine  projects  ellipsoids  onto  the  projection  plane 
and  circumscribes  the  projected  shadow  of  the  ellipsoid  with  a 
rectangle.  This  rectangle  is  used  later  in  the  program  by  the 
BUILDIE  routine  which  checks  for  overlaps  of  objects  on  the 
projection  plane.  The  information  obtained  from  these  opera¬ 
tions  is  used  by  the  hidden  line  routines  to  decrease  computa¬ 
tion  time. 

The  PRJELR  routine  assumes  that  an  ellipsoid  projects  onto  the 
projection  plane  as  an  ellipse,  but  in  general  this  is  not  the 
case.  The  assumption  used  in  writing  the  VIEW  program  is  that 
when  the  viewpoint  is  sufficiently  far  away  from  the  subject  and 
the  viewpoint  coordinate  system  is  looking  almost  directly  at 
the  subject,  all  ellipsoids  will  project  approximately  as 
ellipses. 

b.  Subroutines  Called 

ABS,  SQRT,  DOT,  DOTT,  MAT,  SOLVA,  SOLVR 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  POLYGON 

d.  Input  or  Argument  Parameters 
None 
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e.  Optional  Output 
None 

f.  Procedure 

Immediately  upon  entering,  PRJELR  checks  the  number  of 
segments.  If  zero,  the  subroutine  exits.  If  nonzero,  It  con¬ 
tinues.  The  following  procedure  Is  done  for  each  segment.  The 
A  array  is  transformed  Into  the  viewpoint  reference  frame. 
Subroutine  SOLVR  is  called  three  times  to  get  the  vectors 
n,  r2,  and  ?3  as  described  In  Appendix  B.  Vectors  r\,  r2. 
and  r3  are  projected  onto  the  projection  plane.  Subroutine 
SOLVA  is  called  to  get  the  coefficients  of  the  ellipse  matrix. 
Eigenvalues  LAMDA1  and  LAMDA2  are  solved.  The  two  Eigenvalues 
are  constructed  and  normalized.  The  SIGN  and  CONVEC  arrays  are 
set  up.  After  doing  this  procedure  for  each  segment,  the  sub¬ 
routine  exits. 

4.25  SUBROUTINE  PRJPLY 

a.  Purpose 

PRJPLY  is  called  by  the  Main  program  to  project  polygons  in 
three  space  to  the  two-dimensional  space  of  the  projection 
plane. 

b.  Subroutines  Called 

MAT 

c.  Labeled  Common  Blocks  Used 
ELLIPSE,  POLYGON 
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d.  Input  or  Argument  Parameters 
None 

e.  Optional  Output 
None 

f.  Procedure 

PRJPLY  works  with  one  plane  at  a  time  and  fills  the  POS,  CONVEC, 
and  SIGN  arrays  for  that  plane;  then  the  next  plane  is  used 
until  all  planes  have  been  converted.  A  pictorial  representa¬ 
tion  of  the  meaning  of  the  vectors  in  the  arrays  is  given  in 
Figure  4.5.  The  ?  vectors  are  given  in  the  inertial  reference 
frame.  These  vectors  define  the  corners  of  the  polygon  in  three 
space.  These  ?  vectors  are  projected  onto  the  projection  plane, 
and  the  CONVEC  array  is  generated.  The  CONVEC  array  consists  of 
contour  vectors  of  the  projected  polygon.  The  POS  array  con¬ 
tains  a  vector  that  defines  one  corner  of  the  projected  polygon 
from  a  coordinate  system  aligned  with  the  viewpoint  coordinate 
system.  The  SIGN  array  contains  the  result  of  CONVEC ( 1 )  crossed 
with  CONVEC (2) .  The  CONVEC  and  SIGN  arrays  are  used  by  the 
BUILDIE  block. 

4.26  SUBROUTINE  PSE 
a.  Purpose 

This  subroutine  is  called  from  the  Main  program  to  set  up  arrays 
containing  semiellipsoid  data  for  subroutine  PNTPLT  to  plot. 
Subroutine  PSE  is  called  twice  for  each  ellipsoid,  once  to  plot 
its  top  half,  the  second  time  to  plot  the  bottom  half. 


b.  Subroutines  Called 


PNTPLT 

c.  Labeled  Common  Blocks  Used 
None 

d.  Input  or  Argument  Parameters 

XI  —  Semiellipsoid  contour  array. 

IN  --  Number  of  points  in  each  quarter  contour  saved  in 

array  XI. 

SEG  —  Array  containing  a  complete  contour. 

INDEX  —  Number  of  steps  plus  one. 

INDEX2  —  Maximum  number  of  points  any  complete  contour  can 
have. 

IHALF  —  Flag  controlling  which  ellipsoid  half  is  to 
be  plotted. 

IHALF  =  1,  semiellipsoid  with  X  >  0  is  plotted. 
IHALF  =  2,  semiellipsoid  with  X  <  0  is  plotted. 

e.  Procedure 

Contours  are  plotted  in  sequence  starting  with  the  contour  that 
is  represented  by  just  a  point  (i .e.,  Y  *  0,  Z  =  0,  X  >  0)  until 
all  contours  have  been  plotted.  The  last  point  plotted  is 
represented  by  a  point  (1 .e.,  Y  =  0,  Z  =  0,  X  <  0). 

4.27  SUBROUTINE  ROT 

a.  Purpose 

Computes  rotation  matrix  A^  for  angle  TH  about  X,  Y,  or  Z  axes  as 
L  =  1,  2,  or  3. 
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b.  Subroutines  Called 


COS,  SIN 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

A  —  3x3  rotation  matrix  to  be  computed. 

L  —  1,  2,  or  3  indicating  rotation  about  X,  Y,  or  Z 

axes,  respectively. 

TH  —  Angle  of  rotation  8,  in  radians. 

e.  Optional  Output 
None 


f.  Procedure 


1.  For  L  =  1,  computes 


0  0 

cos  8  sin  8 

-sin  8  cos  8 


2.  For  L  *  2,  computes 


cos  0 

0 

-sin 

0 

1 

0 

sin  0 

0 

cos 

0  ) 
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For  L  =  3,  computes 


0 

0 

0 


(cos  6  sin  6 

-sin  0  cos  0 

0  0 

Note:  Special  tests  are  performed  to  insure  that  cos  0  and 
sin  0  are  exactly  0  or  ±1  for  values  of  0  that  are  multiples  of 
^  to  correct  for  small  errors  introduced  by  the  SIN  and  COS 
routines. 

4.28  SUBROUTINE  SOLVA 

a.  Purpose 

SOLVA  solves  three  equations  simultaneously  and  returns  the 
components  of  [a].  See  Appendix  B,  "Discussion  of  Equations 
Used  by  PRJELR"  for  more  information. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 


R 

Contains  values  of  Ti 

the  projection  plane. 

AA11 

an  component  of  [a] 

AA22 

a22  component  of  [a] 

AA12 

a12  component  of  [a] 

r2,  and  r3  on 
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e.  Optional  Output 
None 

f.  Procedure 

See  Appendix  B,  "Discussion  of  Equations  Used  by  PRJELR"  for 
procedural  Information. 

4.29  SUBROUTINE  SOLVR 

a.  Purpose 

SOLVR  solves  a  set  of  simultaneous  equations  to  find  the 
components  of  vector  r  that  satisfy  the  properties  needed  to 
determine  the  equation  of  the  projected  ellipse.  For  more 
information,  refer  to  Appendix  B,  "Discussion  of  Equations  Used 
by  PRJELR." 

b.  Subroutines  Called 
SQRT 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 

CALL  SOLVR  (Ai,  A2,  A3,  A4 ,  A5,  A6,  A7,  A8,  SS,  Rl,  R3) 


Case  No.  1 

Case  No.  2 

Case  No.  3 

< 

II 

<C 

Ai  =  A'  12 

Ai  =  A'  n  +  A'  12 

CM 

< 

It 

CM 

C 

CM 

CM 

II 

CM 

A2  =  A'  21  +  A'  22 
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A3  =  A' 31 

A3  a  A  *  32 

-T 

II 

> 

to 

A4  =  A* 13 

CD 

CM 

II 

in 

c 

As  =  A 1 23 

A6  =  A 1 33 

A6  =  A  *  3  3 

ti 

*—• 

CM 

CM 

<: 

11 

<c 

A8  -  A* 13 

As  -  A 1 23 

SS  =  ss 

ss  =  ss 

R1  =  X  component 
of  ri 

Rl  =  Y  component 
Of  T2 

R3  =  Z  component 
of  n 

Rl  =  Z  component 
of  r2 

e.  Optional  Output 
None 


A3  =  A' 31  +  A'  32 


A4  =  A* 13 


As  =  A 1 23 


As  =  A*  33 


A7  =  A'  u  +  2A *  12+  A' 
As  =  A' 13  +  A'  23 

ss  =  ss 

R1  =  X  and  Y  component 
of  r3 

R3  =  Z  component 
of  r3 


f.  Procedure 

Refer  to  Appendix  8,  "Discussion  of  Equations  Used  by  PRJELR" 
for  procedural  information. 

4.30  SUBROUTINE  TITLE 

a.  Purpose 

The  purpose  of  this  subroutine  is  to  read  the  title  data  cards 
and  write  them  to  the  plot  file. 


22 
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b.  Subroutines  Called 

NEWPEN,  NFRAME,  PLOT,  SYMBOL 

c.  Labeled  Common  Blocks  Used 
DBUG 

d.  Input  or  Argument  Parameters 
Input  cards  2.0  and  2.1  (20) 

e.  Optional  Output 
None 

f.  Procedure 

After  initializing  some  variables,  the  pen  is  moved  (up)  to 
point  0,0  and  is  defined  as  the  origin.  Data  card  2.0  is 
read.  If  number  of  title  frames  equal  to  zero,  the  subroutine 
returns.  The  following  procedure  is  repeated  for  each  title 
frame.  The  set  of  20  data  card  2.1  's  are  read  in  and  SYMBOL  is 
called  to  write  them  to  the  plot  file.  Depending  on  the  value 
of  DEVFLG,  either  PLOT  or  NFRAME  is  called  to  ready  the  plot 
file  for  the  first  flame.  After  doing  all  the  frames,  the  sub¬ 
routine  exits. 

4.31  SUBROUTINE  TPOINT 

a.  Purpose 

This  subroutine  tests  a  point  against  a  polygon,  both  being  on 
the  projection  plane.  The  results  of  the  tests  indicate  If  the 
point  lies  inside  or  outside  of  the  polygon.  IN  is  a  flag  that 
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is  returned  to  the  calling  program  to  Indicate  the  final 
result. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 
POLYGON 

d.  Input  or  Argument  Parameters 

PP2  --  Point  on  the  projection  plane. 

I  —  Polygon  number  on  the  projection  plane. 

IN  --  Flag  returned  telling  if  point  was  inside  (IN=1) 

or  outside  (IN=2)  polygon. 


e.  Optional  Output 
None 

f.  Procedure 

The  test  used  is  called  the  CROSS  PRODUCT  TEST  (see 
Figure  4.6).  This  name  is  appropriate  since  the  test  is  based 
upon  the  sign  of  the  result  of  a  cross  product  between  two  vec¬ 
tors  on  the  projection  plane.  One  vector  represents  a  side  of 
the  polygon,  and  the  other  vector  always  extends  from  the  base 
of  the  side  vector  to  the  point  being  tested.  Examples  of  these 
vectors  are  given  in  Figure  4.6.  The  vector  that  is  represented 
by  points  AB  is  crossed  into  the  vector  represented  by  points 
AE.  This  process  is  continued  until  all  sides  have  been  crossed 
with  a  vector  to  the  point  being  examined.  In  the  case  of  point 
number  2,  the  following  cross  products  would  be  examined. 
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Figure  4.6  Cross-Product  Test  with  Convex  Polygon  in  TPOINT 

Aff  x  A?,  B?  x  B?,  Cft  x  C?,  Bit  x  G? 

Notice  that  the  first  cross  product  gives  a  different  sign  than 
the  second.  If  a  point  is  outside  the  polygon,  there  must  be  a 
change  in  the  resulting  sign  of  the  cross  products.  If  the 
point  is  inside  the  polygon,  the  following  cross  products  need 
to  be  examined. 

A$  x  A?,  B?  x  §?,  x  CT,  DX  x  5? 

Notice  that  all  cross  products  have  the  same  sign.  This  will  be 
true  for  any  polygon  except  for  concave  polygons.  In  the  case 
of  concave  polygons  (see  Figure  4.7),  even  though  the  point  is 
inside  the  polygon,  the  test  would  indicate  that  the  point  was 
outside  the  polygon  because  there  is  a  sign  change  between 

A$  x  a£  and  B?  *  8? 
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Figure  4.7  Cross-Product  Test  with  Concave  Polygon  in  TPOINT 

Therefore,  concave  polygons  are  ruled  out.  Concave  polygon 
shapes  can  be  included  only  when  that  concave  polygon  is  repre¬ 
sented  by  convex  polygons. 

4.32  SUBROUTINE  TRANS1 

a.  Purpose 

The  purpose  of  TRANS1  is  to  transform  the  input  vector  ft  into 
the  viewpoint  reference  frame. 

b.  Subroutines  Called 
DOTT,  MAT 

c.  Labeled  Common  Blocks  Used 

» 

ELLIPSE,  VIEWP 
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d.  Input  or  Argument  Parameters 

R  —  Input  vector. 

P  --  Output  vector  (Input  vector  transformed  to 

Viewpoint  Reference  Frame). 

e.  Optional  Output 
None 

f.  Procedure 

The  vectors  used  by  TRANS1  are  shown  in  Figure  4.8.  The  output 
of  TRANS1  is  vector  ft  in  the  Viewpoint  Reference  Frame.  The 
other  vectors  are  in  the  reference  frame  as  specified  below. 

ft  in  Local  Reference  Frame  of  Ellipsoid. 

This  vector  is  a  position  vector  for  surface  points  of  an  ellip¬ 
soid,  and  it  originates  at  the  center  of  the  ellipsoid. 

SEGLft  in  Inertial  Reference  Frame. 

This  vector  is  a  position  vector  for  the  center  of  the  contact 
ellipsoid  associated  with  ft. 

vft  in  the  Inertial  Reference  Frame. 

This  vector  is  a  position  vector  for  the  origin  of  the  Viewpoint 
Coordinate  System. 

ft  can  be  found  from  vectors  ft,  SEGLft,  and  Vft  once  they  are  all 
in  the  Viewpoint  Reference  Frame.  The  following  equations 
transform  these  vectors  into  the  Viewpoint  Reference  Frame. 
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fEGlP 


CENTER  OF 
CONTACT 
ELLIPSOID 


VIEWPOINT 

COORDINATE 

SYSTEM 


^  INERTIAL  REFERENCE  FRAME 
Z 


Figure  4.8.  Vectors  Used  by  TRANS1 


r2  =  [DVP]  [dt]  ft 
SEGLPft  =  [DVP]  SEGL? 

Vf2  =  [DVP ]  V? 

Note:  [DVP]  is  the  direction  cosine  matrix  that  transforms  from 
the  inertial  to  the  viewpoint  frame  of  reference. 

[D]  is  the  direction  cosine  matrix  that  transforms 
from  the  inertial  to  the  viewpoint  frame  of  reference. 
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Then  P  Is  given  by 


?  =  SEGLP?  +  R?  -  VP? 

IF  IELP  is  greater  than  30,  TRANS1  must  work  with  a  vector  on  a 
polygon  and  not  an  ellipsoid.  The  vectors  that  define  the  sides 
of  the  polygon  are  always  in  the  Inertial  Reference  Frame; 
therefore,  DVP  is  placed  in  DD_  matrix,  and  the  same  equations 
are  used  to  find 

4.33  FUNCTION  XINTCP 

a.  Purpose 

The  purpose  of  XINTCP  is  as  follows:  given  two  sets  of  X  and  Y 
coordinates  and  a  Y  coordinate  that  falls  between  them,  XINTCP 
calculates  the  X  coordinate  at  the  given  Y  value. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 

d.  Input  or  Argument  Parameters 


X 

X  coordinate  of  end  point 

No.  1. 

Y 

Y  coordinate  of  end  point 

No.  1. 

XSAV 

X  coordinate  of  end  point 

No.  2. 

YSAV 

Y  coordinate  of  end  point 

No.  2. 

YTEMP 

Y  coordinate  of  the  point 

wants  the  X  coordinate. 

at  which  the  caller 
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e.  Optional  Output 
None 

f.  Procedure 

XINTCP  has  two  input  pairs  of  coordinates.  In  this  discussion 
they  will  be  referred  to  the  saved  (XSAV,  YSAV)  and  present 
(X,Y)  coordinates  (see  Figure  4.9).  XINTCP  first  calculates  the 
differences  between  the  saved  and  present  X  and  Y  coordinates. 
These  differences  are  called  XI  and  Yl.  It  then  calculates  the 
scaling  factor  SFACTR  as 

XI 

YT 

If  Yl  is  equal  to  zero,  PFACTR  is  set  to  zero.  The  difference 
between  the  saved  and  intercept  Y  value  is  calculated  (Y2). 

This  difference  is  multiplied  by  PFACTR  and  added  to  the  saved  X 
value.  This  is  the  X  intercept  value.  It  is  set  equal  to 
XINTCP  and  the  function  returns. 

4.34  SUBROUTINE  XYZ 

a.  Purpose 

XYZ  solves  the  equation  for  an  ellipse  to  find  a  point  on  that 
ellipse  that  will  return  a  vector.  It  is  used  when  and  M2 
are  not  equal  to  zero.  Refer  to  Appendix  A,  "Hidden  Line 
Problem  Between  Two  Ellipsoids." 
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Figure  4.9.  Variables  Used  by  XINTCP 

Subroutines  Called 
SQRT 


Labeled  Common  Blocks  Used 


d. 


Input  or  Argument  Parameters 


MU 

A 

B 

C 

S 

M 

JFLAG 


L  See  Appendix  A  equations  for  definitions. 


Flag  passed  back  as  either  zero  when 
equation  solved  or  one  when  not  solved. 


e.  Optional  Output 
None 

f.  Procedure 

Follow  equations  in  Appendix  A  for  procedural  information. 

4.35  FUNCTION  YINTCP 

a.  Purpose 

The  purpose  of  YINTCP  is  as  follows:  given  two  sets  of  X  and  Y 
coordinates  and  an  X  coordinate  that  falls  between  them,  YINTCP 
calculates  the  Y  coordinate  at  the  given  X  value. 

b.  Subroutines  Called 

None 

c.  Labeled  Common  Blocks  Used 

None 


-f* *7  ' 
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d. 


Input  or  Argument  Parameters 


X 

X  coordinate  of  point  No. 

Y 

Y  coordinate  of  point  No. 

XSAV 

X  coordinate  of  point  No. 

YSAV 

Y  coordinate  of  point  No. 

XTEMP 

X  coordinate  of  the  point 

wants  the  Y  coordinate. 


1. 

1. 

2. 

2. 

at  which  the  caller 


e.  Optional  Output 
None 

f.  Procedure 

YINTCP  has  two  input  pairs  of  coordinates.  In  this  discussion 
they  will  be  referred  to  as  the  present  (X,Y)  and  saved  (XSAV, 
YSAV)  coordinates  (see  Figure  4.10).  YINTCP  first  calculates 
the  differences  between  the  present  and  saved  X  and  Y  coordi¬ 
nates.  These  differences  are  called  XI  and  Yl.  It  then  calcu¬ 
lates  the  scaling  factor  SFACTR  as 


Yl 

71 

If  Yl  is  equal  to  zero,  PFACTR  is  set  to  zero.  The  difference 
between  the  saved  and  Intercept  X  values  is  calculated  (X2). 
The  difference  (X2)  is  multiplied  by  PFACTR  and  added  to  the 
saved  Y  coordinate.  This  is  the  Y  intercept  value.  It  is  set 
equal  to  YINTCP  and  the  function  returns. 
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^'SAVED'  X  +  Y  COORDS. 
v  \ 

s 


\ 


\ 


\ 


KNOWN  Y  COORD. 


‘PRESENT’ 
X  +  Y 


"xT 


j  COORDS. 


Figure  4.10.  Variables  used  by  YINTCP 


4.36  SUBROUTINE  YZ 

a.  Purpose 

YZ  solves  the  equation  for  an  ellipse  to  find  a  point  on  that 
ellipse  that  will  return  a  vector.  Is  used  when  Mi  =  0 
and  M2  =  0.  Refer  to  Appendix  A,  "Hidden  Line  Problem  Between 
Two  Ellipsoids." 

b.  Subroutines  Called 
SQRT 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  Parameters 


MV 

A 

B 

C 

S 

M 

JFLAG 


l  See  Appendix  A  equations  for  definitions. 


Flag  passed  back  as  either  zero  when  equation 
solved  or  one  if  not  solved. 


e.  Optional  Output 
None 

f.  Procedure 

Follow  equations  in  Appendix  A  for  procedural  information  (Case 
No.  2). 


4.37  SUBROUTINE  Z 

a.  Purpose 

Z  solves  the  equation  for  an  ellipse  to  find  a  point  on  that 
ellipse  that  will  return  a  vector.  Is  used  when  both  Mi  and  M2 
are  equal  to  zero.  Refer  to  Appendix  A,  "Hidden  Line  Problem 
Between  Two  Ellipsoids." 

b.  Subroutines  Called 
SQRT 

c.  Labeled  Common  Blocks  Used 

None 
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d.  Input  or  Argument  ’arameters 


MV 

A 

B 

C 

S 

M 

JFLAG 


See  Appendix  A  equations  for  definitions. 


Flag  passed  back  as  either  zero  when  equation 
solved  or  one  if  not  solved. 


e.  Input  or  Argument  Parameters 


None 


f.  Procedure 


Follow  equations  in  Appendix  A  for  procedural  information  (Case 
No.  1). 


5.0  VIEW  PROGRAM  FLOWCHARTS 


Contained  In  this  section  are  the  flowcharts  for  selected  program 
modules.  They  are  not  direct,  line-by-line  flowcharts,  but  more  of  a 
general  flow  description  of  critical  program  modules. 


5.1  VIEW  MAIN 


© 

© 

© 

© 

© 

© 

© 

© 
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5.2  BUILDIE 


RETURN 
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1 


PIP 


5.3  CLIP 
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4 


CONVREC 


RETURN 
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5.5  ELIPSN 


5.6 


EXTEND 
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FILL  IN  THIRD  ROW  OF  DIRECTION 
COSINE  MATRIX 


E 


YES 


EZ 

INDICATE  POINT 
NOT  HIDDEN 
I  FLAG  =  2 


TRANSFORM  POSITION  VECTOR 
P3  TO  THE  VIEWPOINT 
REFERENCE  FRAME 

PROJECT  POSITION  VECTOR 

ON  TO  THE  PROJECTION  PLANE 

p 

CALL  TPOINT 

TEST  POINT  AGAINST  POLYGON 
REPRESENTED  BY  KK 

IS  POINT  HIDDEN  ? 

NO 

DOES  IFLAG  =  1  ? 

f  UNDETERMINED 

VlETURN/ 

IS  POINT  CLOSER  TO  VIEW- 

NO 

\  / 

POINT  THAN  THE  POLYGON  ? 

V 

INDICATE  POINT 
HIDDEN 
I  FLAG  =  1 


IF  FIRST  CALL 


READ  CARDS  6.0,  6.1  FROM  LUN  #5  AND  POLYGON  DATA  FROM  LUN  #1 
READ  ELLIPSOID  DATA  FROM  LUN  #1  TILL  ZTIME  2  CTIME 
IF  NPL»0,  READ  CARDS  7.1  AND  7.2  FOR  EACH  PLANE 
INITIALIZE  ARRAYS  A  AND  SEGLP 
IF  NSEG>0,  READ  CARD  8.0  FROM  LUN  #5 
IF  NPL>0,  READ  CARDS  9.0  AND  10.0  FROM  LUN  #5 
READ  CARDS  11.0,  12.0,  13.0,  14.0  AND  15.0  FROM  LUN  #5 
IF  ICODE  =  0,  CALL  DRCYPR 
IFICODE  =  2,  CALL  GENDCM 

IF  !COO"  *  0  AND  ICODE  *  2,  READ  CARD  16.0  FROM  LUN  #5 
STORE  INITIAL  VALUES  OF  VP  +  DVP  ARRAYS 


RETURN 


IF  NOT  FIRST  CALL 


READ  SEGMENT  DATA  FROM  LUN  #1  TILL  ZTIME  S  CTIME 


WHEN  CTIME  REACHED,  CALCULATE  LINEAR  POSITION  OF  CONTACT 
ELLIPSOID  IN  INERTIAL  REFERENCE  SYSTEM 


IF  THIS  MATRIX  MULTIPLY  ATTEMPTED  TWICE  WITHOUT  AN 
INTERVENING  READ  FOR  LUN  #1,  SET  (FLAG  =  10 


RETURN 
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5.10  LSEGINT 
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5.11  OVERLAP 


RETURN 


PLPLN 


r 


L. 
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INITIALIZE  last  AND  CUNRfcNT  POINT  FLAGS, 

lflaq  and  iflag 


SAVE  X+Y  COORDS. 


5.14  PREPLT 


5.15  POLYD 


RETURN 


RETURN 
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5.16  PRJELR 


RETURN 


RETURN 
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5.17  PRJPLY 


5.18  PSE 
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5.19  XINTCP 
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5.20  YINTCP 
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T 
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6.0  VIEW  PROGRAM  VARIABLE  GLOSSARIES  BY  PROGRAM  MODULE 


This  section  contains  lists  of  selected  variables  contained  In  VIEW.  The 
lists  contain  all  pertinent  variables  and  exclude  temporary  variables  such 
as  DO  loop  Indices  and  Intermediate,  calculated  variables.  This  section  Is 
organized  In  the  following  manner;  common  block  variables,  main  program 
variables,  and  variables  for  each  subroutine,  subroutines  In  alphabetical 
order.  The  variables  for  all  common  blocks  have  been  combined  and  listed 
In  alphabetical  order.  For  which  common  block  the  variable  belongs,  see 
Section  7,  cross-reference  Chart  7.3.  The  variable  descriptions  are  alpha¬ 
betized  within  each  subroutine.  The  variable  description  consists  of  the 
variable  name;  Its  dimension,  If  any;  FORTRAN  data  type;  and  a  short 
description. 


6.1  COMMON  BLOCK  VARIABLE  DEFINITIONS 


A(3,3,30) 


Real  Array.  Contains  the  ellipsoid 
matrices  used  In  the  following  equation: 

r  •  [A]  r  »  1 


BDRS 

CONVEC (2,5,90) 


D(3,3,90) 


DEVFLG 


Integer  scalar.  Constant  that  represents 
a  particular  graphics  output  system. 

Used  with  DEVFLv  .  Initialized  in  Main. 

Real  array.  Contains  vectors  on  the 
projection  plane  that  represent  the 
projected  polygons  and  rectangles  that 
circumscribe  projected  ellipsoids  on  the 
projection  plane.  Initialized  In  sub¬ 
routine  PRJPLY. 

Real  array.  Contains  directional  cosine 
matrices  for  segments  (third  subscript 
value  1-30)  and  polygons  (31-90).  Data 
for  the  segments  Is  read  from  the  ATB 
Input  file  In  subroutine  INPUT  and  data 
for  polygons  Is  Initialized  In  subroutine 
POLYD.  Transforms  from  the  inertial 
coordinate  system  to  the  local  geometric 
coordinate  sytem  of  the  segment. 

Integer  scalar.  Represents  VIEW  program 
output  device  number.  Read  from  Input 
control  file  In  Main. 
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DVP(3,3) 

DVPd (3,3) 

ICOLOR  (91) 

IDEBUG(80) 

IE (90,90) 

IELP 

I  FLAG 

I  NT 

I REMOV ( 30 ) 


Real  array.  Contains  the  direction 
cosine  matrix  for  the  viewpoint  coordi¬ 
nate  system  that  transforms  from  the 
Inertial  reference  system  to  the  view¬ 
point  coordinate  system.  Initialized  In 
subroutine  INPUT  and  modified  In  Main. 

Real  array.  Contains  the  original 
direction  cosine  matrix  that  transforms 
from  the  segment  reference  system  to 
which  the  viewpoint  Is  attached,  to  the 
viewpoint  coordinate  system.  Initialized 
In  subroutine  INPUT. 

Integer  array.  Contains  color  numbers 
for  each  ellipsoid  (subscript  values  1- 
30),  polygon  (31-90 )  and  title  (91). 

Read  from  input  control  file  In  sub¬ 
routine  INPUT. 

Integer  array.  Contains  debugging 
printer  control  flags.  Read  from  Input 
control  file  in  Main. 

Integer  array.  Contains  object  numbers 
that  overlap  with  a  given  object.  Second 
subscript  values  1-30  represent  ellip¬ 
soids,  31-90  represents  polygons. 
Initialized  in  subroutine  BUILDIE . 

Integer  scalar.  Represents  ellipsoid  or 
polygon  that  is  currently  being 
plotted.  Initialized  and  modified  In 
Main  and  subroutine  PLPLN. 

Integer  scalar.  Flag  Indicating  state  of 
subroutine  INPUT.  Initialized  and  modi¬ 
fied  In  Main  and  INPUT. 

Integer  scalar.  Number  of  iterations 
used  In  subroutine  EXTEND  for  finding  the 
edge  of  a  hidden  line  segment.  I  NT  =  4 
Is  recommended.  Read  from  Input  control 
file  In  subroutine  INPUT. 

Integer  array.  Contains  numbers  that 
represent  polygons  to  be  removed  from 
View  program  output.  Read  from  Input 
control  file  In  subroutine  INPUT. 
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IVP 


MPL(3,5,60) 


N IE (90) 


NISG 

NP 

NPLANE 


NPPP ( 90 ) 


NPREM 


NSEG 

NSTEPS (90 ) 


Integer  scalar.  Segment  number 
representing  ellipsoid,  vehicle,  or 
ground  to  which  the  viewpoint  and  focal 
point  are  attached. 

Integer  array.  Contains  segment  and 
ellipsoid  Identification  numbers  for  each 
plane  segment  contact  as  defined  in  the 
ATB  program.  MPL  is  read  in  from  both 
the  ATB  Input  file  and  the  input  control 
file. 

Integer  array.  Contains  the  number  of 
overlapping  ellipsoids  and  polygons  for  a 
given  ellipsoid  or  polygon.  Subscript 
values  1-30  represent  ellipsoids,  31-90 
represent  polygons.  Initialized  In 
subroutine  BIJILDIE. 

Integer  scalar.  Defines  the  reference 
frame  for  unattached  planes.  Read  from 
Input  control  file  in  subroutine  INPUT. 

Integer  scalar.  Represents  the  number  of 
ATB  input  polygons.  Read  from  ATB  Input 
file  In  subroutine  INPUT. 

Integer  scalar.  Represents  total  number 
of  polygons  (ATB  input  plus  VIEW  input 
planes).  Initialized  in  subroutine 
INPUT. 

Integer  array.  Contains  number  of  points 
per  polygon.  This  equals  number  of  sides 
plus  one.  Initialized  in  subroutine 
PRJELR  for  rectangles  about  ellipsoids 
and  read  from  input  control  file  In 
subroutine  INPUT  for  polygons. 

Integer  scalar.  Represents  the  number  of 
polygons  to  be  removed  from  plot  out¬ 
put.  Read  from  Input  control  file  In 
subroutine  INPUT. 

Integer  scalar.  Total  number  of  segments 
read  from  ATB  Input  file  in  subroutine 
INPUT. 

Integer  array.  Contains  number  of  steps 
(grid  points  along  an  axis)  for  each 
ellipsoid  and  polygon.  In  the  case  of 
the  polygons,  NSTEPS  represents  the 
number  of  vectors  on  any  one  side  of  the 
polygon.  Read  from  Input  control  file  In 
subroutine  INPUT. 
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OFLINE 

OFSETX 

OFSETY 

ONLINE 

P (3,4,60) 

PL ( 1 7 ,30 ) 

POS (2 , 90 ) 

P0(3,4,60) 

RA  (3 ) 


Integer  scalar.  Constant  that  represents 
off-line  condition.  Used  with  DEVFLG. 
Initialized  In  Main. 

Real  scalar.  Offset  variable  to  move 
plot  In  X  direction  on  the  plotting 
area.  Read  from  Input  control  file  In 
subroutine  INPUT. 

Real  scalar.  Offset  variable  to  move 
plot  In  Y  direction  on  the  plotting 
area.  Read  from  Input  control  file  In 
subroutine  INPUT. 

Integer  scalar.  Constant  that  represents 
online  condition.  Used  with  DEVFLG. 
Initialized  in  Main. 

Real  array.  Contains  vectors  that 
represent  the  polygon  sides  in  the  iner¬ 
tial  reference  frame.  Initialized  in 
subroutine  CONVREC. 

Real  array.  Contains  polygon  parameters 
for  ATB  Input  polygons.  Read  from  ATB 
input  file. 

Real  array.  Contains  position  vectors 
for  the  polygons  represented  by  the 
CONVEC  array.  These  position  vectors 
originate  from  the  viewpoint  origin  on 
the  projection  plane  to  a  corner  of  the 
polygon  stored  in  the  CONVEC  array.  POS 
vectors  are  In  the  viewpoint  coordinate 
system.  POS  is  initialized  in  subroutine 
PRJPLY. 

Real  array.  Contains  position  vectors  of 
the  polygon  vertices  in  reference  to  the 
MPL  coordinate  system.  ATB  input  poly¬ 
gons  are  Initialized  in  subroutine 
CONVREC,  and  VIEW  input  polygons  are  read 
from  Input  control  file  In  subroutine 
INPUT. 

Real  array.  Contains  three  different 
kinds  of  data  depending  on  the  value  of 
ICODE.  These  data  are: 

ICODE  =0  Roll,  pitch,  and  yaw 

angles  are  supplied  In  RA 
array.  Initialized  In 
subroutine  DRYCYPR. 


97 


4 


ICOOE  *  1 


Direction  cosine  matrix 
supplied  as  Input.  The  RA 
array  Is  the  first  row  of 
the  direction  cosine 
matrix.  The  second  and 
third  rows  are  on  the  next 
record.  This  data  is  read 
from  the  Input  control 
file  In  subroutine 
INPUT. 

ICODE  =  2  Point  at  which  viewpoint 
Z-axis  to  aim  Is  supplied 
In  the  RA  array.  Read 
from  input  control  file  In 
subroutine  INPUT. 

SEGLP(3,90)  Real  array.  Contains  position  vectors 

for  ellipsoids  and  polygons.  These 
vectors  go  from  the  inertial  reference 
origin  to  the  center  of  the  contact 
ellipsoid.  In  the  case  of  a  polygon, 
this  vector  points  to  one  corner  of  the 
polygon.  SEGLP  is  in  the  inertial  refer¬ 
ence  frame.  It  is  read  from  the  ATB 
Input  file  In  subroutine  INPUT. 

SFACTR  Real  scalar.  Reprsents  scale  factor  for 

plot.  Is  the  Z  coordinate  of  the  pro¬ 
jection  plane  in  the  viewpoint  reference 
frame.  Read  from  input  control  file  in 
subroutine  INPUT. 

SIGN(90)  Real  array.  Contains  cross  product  of 

side  1  and  side  2  of  a  polygon.  Initial¬ 
ized  In  subroutine  PRJPLY. 

TERM  Integer  scalar.  Constant  that  represents 

terminal  as  output  device.  Used  with 
DEVFLG.  Initialized  in  Main. 

TIME  Real  scalar.  Time  (seconds)  of  the 

current  segment  data  set.  Read  from  ATB 
Input  file  In  subroutine  INPUT. 

VP(3)  Real  array.  Position  vector  for  the 

viewpoint  coordinate  system.  This  vector 
Is  In  the  Inertial  reference  frame.  Read 
from  Input  control  file  In  subroutine 
INPUT  and  modified  In  Main. 
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VP0(3) 

Real  array.  Contains  the  coordinates  of 
the  viewpoint  In  the  coordinate  system  to 
which  the  viewpoint  Is  attached.  These 
are  the  original  values  read  Into  array 
VP.  Initialized  In  subroutine  INPUT. 

VP2 ( 3 ) 

Real  array.  Contains  the  viewpoint 
position  vector  In  the  viewpoint  coordi¬ 
nate  system. 

ZTIME 

Double  precision  scalar.  Represents  time 
of  current  data  set  (TIME)  rounded  to  the 
nearest  microsecond.  This  was  done  in 
order  to  get  equal  comparisons  between 
Input  data  time  and  program  calculated 
time.  Modified  in  subroutine  INPUT. 

MAIN  PROGRAM 

CTIME 

Double  precision  scalar.  Represents 
current  time  of  the  program  rounded  to 
the  nearest  microsecond. 

DTIME 

Double  precision  scalar.  Represents 
delta  time  or  what  time  step  value  CTIME 
will  be  incremented.  Read  from  record 

4.0  of  input  control  file. 

ETIME 

Double  precision  scalar.  Represents  end 
time  of  the  program.  Read  from  record 

4.0  of  input  control  file. 

ID(IO) 

Real  array.  Contains  40  character  title 
strip.  Read  from  record  3.0  of  input 
control  file. 

LUPLOT 

Integer  scalar.  Represents  logical  unit 
number  for  the  output  plotting  file. 

NF 

Integer  scalar.  Frame  counter  that  is 
written  every  frame  to  output  list 
file. 

STIME 

Double  precision  scalar.  Represents 
start  time  of  the  program  or  when  to 
begin  plotting.  Read  from  card  4.0  of 
input  control  file. 

XTIME 

Real  scalar.  Temporary  variable  that  is 
the  current  program  time  expressed  In 
milliseconds  for  plotting. 
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WORK (10000) 

6.3  SUBROUTINE  BUILD IE 
I  OBJ 

KPLANE 

MFLAG 

MPLANE 

6.4  SUBROUTINE  CLIP 

IPEN 

I  PLOT 

LCALL 

X 

XLIMIT 

XLSAV 

XLTEMP 


Real  array.  Temporary  array  space 
allocated  for  subroutines  ELIPSN,  PSE, 
and  PLPLN. 


Integer  scalar.  Equal  to  present  object 
number  plus  one. 

Integer  scalar.  Equal  to  total  number  of 
polygons  plus  30. 

Integer  scalar.  Flag  passed  back  from 
subroutine  OVERLAP  telling  whether  or  not 
there  is  a  overlap  condition. 

Integer  scalar.  Equal  to  total  number  of 
polygons  plus  29. 


Integer  scalar.  Calcomp  pen  value  passed 
through  argument  list  from  subroutine 
PNTPLT. 

Integer  scalar.  Flag  that  tells  CLIP  if 
last  pen  move  was  clipped  (IPL0T=0),  or 
plotted  ( I  PLOT =1 ) . 

Integer  scalar.  Flag  that  tells  CLIP  if 
last  call  to  CLIP  was  a  pen  up  to  first 
point  in  line.  LCALL=0  means  it  was, 
LCALL=1  means  that  it  was  not. 

Real  scalar.  X  coordinate  of  point  to  be 
plotted  passed  through  argument  list  from 
PNTPLT. 

Real  scalar.  If  X  coordinate  off 
plotting  region,  this  variable  contains 
value  of  which  boundary  (left  or  right) 

It  Is  off. 

Real  scalar.  X  coordinate  of  first  point 
In  line  that  was  clipped  In  last  call  to 
CLIP. 

Real  scalar.  Value  of  X  coordinate  to 
move  to  If  previous  call  to  CLIP  was  a 
pen  up  to  first  point  In  segment. 
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XMAX 


XMIN 

XOFF 

XSAV 

XTEMP 

Y 

YLIMIT 

YLSAV 

YLTEMP 

YMAX 

YHIN 

YOFF 

YSAV 


Real  scalar.  Constant  that  represents 
right  side  X  limit  of  plotting  region. 
Read  from  record  16.0  In  input  control 
file. 

Real  Scalar.  Constant  that  represents 
left  side  X  limit  of  plotting  region. 

Read  from  record  16.0  In  Input  control 
file. 

Logical  scalar.  Is  .FALSE.  If  X 
coordinate  within  valid  plotting  region, 
.TRUE.  If  outside  plotting  region. 

Real  scalar.  X  coordinate  of  last  plot 
move.  Passed  through  argument  list  from 
PNTPLT. 

Real  scalar.  X  coordinate  of  end  point 
of  line  segment  to  be  clipped. 

Real  scalar.  Y  coordinate  of  point  to  be 
plotted.  Passed  through  argument  list 
from  subroutine  PNTPLT. 

Real  scalar.  If  Y  coordinate  of  plotting 
region,  this  variable  contains  value  of 
which  boundary  (top  or  bottom)  It  is 
off. 

Real  scalar.  Y  coordinate  of  first  point 
In  line  that  was  clipped  In  last  call  to 
CLIP. 

Real  scalar.  Value  of  YMIN  or  YMAX, 
depending  on  if  first  point  in  line 
clipped  was  off  top  or  bottom  of 
plotter. 

Real  scalar.  Constant  that  represents 
upper  (top)  limit  of  the  plotting 
region.  Defined  in  subroutine  PNTPLT. 

Real  scalar.  Constant  that  represents 
lower  (bottom)  limit  of  the  plotting 
region.  Defined  in  subroutine  PNTPLT. 

Logical  scalar.  Is  .FALSE,  if  Y 
coordinate  within  valid  plotting  region, 
.TRUE.  If  outside  plotting  region. 

Real  scalar.  Represents  Y  coordinate  of 
previous  point  plotted.  Passed  through 
parameter  list  from  subroutine  PNTPLT. 
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YTEMP 


Real  scalar.  Value  of  YMIN  or  YMAX, 
depending  on  If  point  clipped  Is  off  top 
or  bottom  of  the  plotter. 


SUBROUTINE  CONVREC 

ODD 

Real  scalar.  Intermediate  variable  used 
In  calculating  data  for  P0  array. 

OX  ( 3 ) 

Real  array.  Contains  results  from 
determinant  function  DET  for  solving  the 
three  equations  CONVREC  uses. 

ISG 

Integer  scalar.  Segment  number  that 
defines  the  coordinate  system  for 
plotting  that  particular  plane. 

NPPPP 

Integer  scalar.  Number  of  points  In  the 
polygon  Loaded  from  array  NPPP. 

R  (3 ) 

Real  array.  Contains  matrix  multiple  of 
portions  of  arrays  DVP  and  P0. 

SUBROUTINE  CROSS 

A  (3 ) 

Real  array.  Vector  In 

B  (3 ) 

Real  array.  Vector  In 

C  ( 3 ) 

Real  array.  Result  vector  In 

FUNCTION  DET 

All  - 

A12 

A13 

A21 

A22 

A23 

A31 

A32 

A33  -  - 

Real  scalars. 

Values  Representing  3x3.  square  array. 

OET 

Real  scalar.  Determinant  of  Input 
array. 
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6.8  SUBROUTINE  DOT  and  DOTT 

• 

A(L,3) 

Real  array.  Array  A  in  matrix  multiply 
C_=AB_. 

B(l,3) 

; 

Real  array.  Array  B  in  matrix  multiply 
C=AB_. 

C(N,M) 

Real  array.  Output  of  DOT(T),  in  array 
matrix  multiply  C=AB. 

C 

L 

Integer  scalar.  First  subscript  value 
for  arrays  A  and  B. 

1  M 

Integer  scalar.  Second  subscript  value 
for  array  C. 

N 

Integer  scalar.  First  subscript  value 
for  array  C. 

6.9  SUBROUTINE  ORCYPR 

A  (3 ) 

Real  array.  Contains  rotation  angles 
(degrees). 

0(3.3) 

Real  array.  Output  of  DRCYPR,  contains 
direction  cosine  matrix. 

11 

Integer  scalar.  Axis  of  rotation  for 
first  angle  (X=l,  Y=2,  Z=3). 

12 

Integer  scalar.  Axis  of  rotation  for 
second  angle  (X=l,  Y=2,  1=3). 

13 

Integer  scalar.  Axis  of  rotation  for 
third  angle  (X=l,  Y=2,  1=3). 

’  M 

Integer  scalar.  Constant  of  6,  size  of 
array. 

x 

N 

Integer  scalar.  Constant  of  3,  size  of 
array. 

x 

.  P 

Real  scalar.  Pitch  angle  in  radians. 

R 

Real  scalar.  Roll  angle  in  radians. 

,  RADIAN 

Real  scalar.  Constant  for  degrees  to 
radians  conversion. 
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6.10  SUBROUTINE  ELIPSN 

DELTAX 

DELTAY 

IN(INDEX) 

INDEX 

N 

SIMP1 

SIMP2 

SIMP3 

TEMP 

TEST 

X 

XI (3,  INDEX,  INDEX2) 


Real  array.  Temporary  buffer  space  used 
for  matrix  multiple. 

Real  scalar.  Yaw  angle  in  radians. 


Real  scalar.  Semiaxes  length  in  X 
direction  divided  by  number  of  steps. 

Real  scalar.  Semiaxes  length  in  Y 
direction  divided  by  number  of  steps. 

Integer  array.  Number  of  points  in  each 
contour  array. 

Integer  scalar.  Number  of  steps  plus 
one. 

Integer  scalar.  Number  of  points  in  this 
half  contour  that  is  about  the  X  axis. 

Real  scalar.  —  where  "a"  is  the 
2 
a 

semiaxes  length  in  the  X  direction. 


Real  scalar.  —  where  "b"  is  the 
2 
b 

semi  axes  length  in  the  Y  direction. 


Real  scalar.  —  where  "c"  is  the 
2 
c 

semiaxes  length  in  the  Z  direction. 
Real  scalar.  The  value  of  1-X  *SIMP1. 


Real  scalar.  The  value  of  TEMP-Y  *SIMP2. 

Real  scalar.  Current  X  coordinate, 
starts  at  zero. 

Real  array.  Semiellipsoid  contour  array. 
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Y 


Real  scalar.  Current  Y  coordinate, 
starts  at  zero. 


Z 

6.11  SUBROUTINE 

I 

I  FLAG 

J 

KK 

N 

NUM 

P (3,2) 

P  3  ( 3 ) 

6.12  SUBROUTINE 
CAMERA (3) 

0(3,3) 


Real  scalar.  Current  Z  coordinate, 
starts  at  length  "C." 


EXTEND 


Integer  scalar.  Points  to  unhidden  point 
location  In  P  array. 

Integer  scalar.  Flag  passed  back  from 
subroutines  HIDE  and  HYDE  telling  whether 
or  not  the  point  is  hidden. 

Integer  scalar.  Points  to  hidden  point 
location  in  P  array. 

Integer  scalar.  Overlapping  object 
number  that  overlaps  with  present  object, 
from  array  IE. 

Integer  scalar.  Denotes  midpoint 
coordinates,  either  I  or  J,  depending  on 
hidden  or  not. 

Integer  scalar.  Represents  the  number  of 
overlapping  objects  with  a  particular 
object. 

Real  array.  Contains  X,  Y,  and  Z 
coordinates  of  the  end  points  of  the 
line. 

Real  array.  Temporary  buffer  containing 
coordinates  of  the  midpoint  of  the  line. 


GENDCM 


Real  array.  Position  vector  for  the 
viewpoint  in  the  reference  sysem  of  the 
segment  to  which  the  viewpoint  is 
attached. 

Real  array.  Direction  cosine  matrix  for 
viewpoint.  This  direction  cosine  matrix 
transforms  from  the  inertial  reference 
frame  to  the  local  reference  frame. 


t 

J 


■ffl  mJm 
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FOCUS (3) 

Real  array. 

Position  vector  for  point 

which  viewpoint  Z  axis  Is  aimed. 

SUM 

Real  scalar. 
In  Z  array. 

Running  summation  of  values 

XNDRM 

Real  scalar. 

Normalization  factor. 

Z(3) 

Real  array. 

Contains  differences  between 

matching  values  in  CAMERA  and  FOCUS 
arrays. 


6.13  SUBROUTINE  HIDE 


IFLAG 

Integer  scalar.  Flag  passed  back  to 
caller  Indicating  hidden  (IFLAG-1),  or 
not  hidden  ( I FLAG*2 ) . 

KK 

Integer  scalar.  Polygon  number  to  check 
blocking  on. 

NPRIME (3) 

Real  array.  Contains  matrix  multiply  of 
arrays  DVP  and  PPRIME. 

P3(3) 

Real  array.  Contains  position  vector  of 
point. 

P4(3) 

Real  array.  Contains  P3  transformed  to 
viewpoint  reference  system. 

P5(3) 

Real  array.  Contains  matrix  multiple  of 
NPRIME  and  PPRIME  arrays. 

P6(3) 

Real  array.  Contains  matrix  multiple  of 
NPRIME  and  P4  arrays. 

P7(2) 

Real  array.  Contains  position  vector  P4 
projected  on  the  projection  plane. 

PP  (3 ) 

Real  array.  Contains  polygon  sides  P 
minus  the  viewpoint  array  VP. 

6.14  SUBROUTINE  HYDE 


IFLAG  Integer  Scalar.  Flag  passed  back  to 

caller  to  indicate  whether  or  not  the 
point  was  hidden.  IFLAG=1  means  hidden. 
IFLAG»2,  not  hidden. 
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N 


Integer  scalar.  Possible  hiding 
ellipsoid  number. 


r (3 )  Real  array.  Vector  to  plotting  point  In 

local  reference  frame  of  ellipsoid  or 
polygon  of  which  It  Is  a  part. 


6.16  SUBROUTINE  INPUT 


BD(24,  40) 

Real  array.  Contains  contact  ellipsoid 
parameters  for  each  contact  ellipsoid. 

The  rotation  of  the  contact  ellipsoid 
relative  to  the  segment  coordinate  system 
Is  already  Incorporated  In  the  values  of 
this  array.  Read  from  ATB  Input  file. 

CTIME 

Double  precision  scalar.  Current  program 
time  (seconds)  calculated  in  Main  pro¬ 
gram. 

DO  (3) 

Real  array.  Contains  the  offset  vector 
(In  the  inertial  coordinate  system)  of 
the  contact  ellipsoid  from  the  segment 
c.g. 

I  CODE 

Integer  scalar.  Input  flag  controlling 
the  generation  of  the  direction  cosine 
matrix.  Read  from  Input  control  file. 

ISW1 

Integer  scalar.  Flag  controlling  whether 
plots  are  to  be  made  when  the  time  of  the 
VIEW  run  is  less  than  the  time  Interval 
of  the  simulation  data.  If  this  Is  true, 
only  the  available  simulation  data  will 
be  plotted. 

NFAST 

Integer  scalar.  Represents  number  of 
segments  to  be  removed  used  In  ATB  simu¬ 
lation.  Read  from  Input  control  file. 

NSIDES 

Integer  scalar.  Number  of  sides  In  a 
polygon.  Loaded  from  NPPP  array. 

NSP 

Integer  scalar.  Defines  number  of 

supplemental  planes  being  Input  from 
Input  control  file. 
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6.16  SUBROUTINE  LSEGINT 


I  FLAG 

Integer  scalar.  Flag  Indicating 
Intersection  ( I FLAG-l ) .  or  no  Inter¬ 
section  (IFLAG-0). 

M(2) 

Real  array.  Contain  the  slopes  of  the 
two  lines. 

P(4,2) 

Real  array.  Contains  coordinates  of 
lines  P1-P2  and  R1-R2. 

Pl(2) 

Real  array.  Contains  X  and  Y  coordinate 
of  one  end  point  of  line  P1-P2. 

P2(2) 

Real  array.  Contains  X  and  Y  coordinate 
of  the  other  end  point  of  line  P1-P2. 

Rl(2) 

Real  array.  Contains  X  and  Y  coordinate 
of  one  end  point  of  line  R1-R2. 

R2(2) 

Real  array.  Contains  X  and  Y  coordinate 
of  the  other  end  point  of  line  R1-R2. 

T  (2) 

Real  array.  Contains  T  ratio  values, 
where 

T  XO  -  XI 

T '  xri  xi 

SUBROUTINE  MAT 

A(JA,1) 

Real  array.  Array  in  matrix  multiply 
C_«AB. 

B(JB.l) 

Real  array.  Array  In  matrix  multiply 
C-AB. 

C(JC.l) 

Real  array.  Array  in  matrix  multiply 
C-AB_. 

JA 

Integer  scalar.  First  subscript  value 
for  array  A. 

JB 

Integer  scalar.  First  subscript  value 
for  array  B. 

JC 

Integer  scalar.  First  subscript  value 
for  array  C. 

LL 

Integer  scalar.  Size  of  array  A. 
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6.18  SUBROUTINE  NFRAME 


Integer  scalar.  Size  of  array  B. 

Integer  scalar.  Size  of  array  C. 

Real  scalar.  Running  summation  of  matrix 
multiple. 


ENOFRA 


Integer  scalar.  End-of-frame  halfword  to 
Grlnnell  graphics  system. 

Integer  scalar.  Dumny  argument  for 
PLOTS. 

Interger  scalar.  Durrwqy  argument  for 
PLOTS. 

Integer  scalar.  Halfword  mask  defining 
what  corresponding  bit  In  ENOFRA  to 
change. 

Integer  scalar.  Dumnjy  argument  for 
PLOTS. 


STATUS 


Integer  scalar.  Contains  return  status 
of  the  request. 


6.19  SUBROUTINE  OVERLAP 


Integer  scalar.  Denotes  object  No.  1  to 
be  tested. 

Integer  scalar.  Same  as  variable  I. 

Integer  scalar.  Denotes  object  No.  2  to 
be  tested. 


MFLAG 


NPTS1 


Integer  scalar.  Same  as  variable  K. 

Integer  scalar.  Indicates  overlap  or 
not.  MFLAG=0  indicates  no  overlap  and 
MFLAG=1  indicates  overlap. 

Integer  scalar.  Number  of  sides  plus  one 
for  object  No.  1. 


NPTS2 


Integer  scalar.  Number  of  sides  plus  one 
for  object  No.  2. 
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Pl(2) 

Real  array.  Contains  X  and  Y  coordinates 
on  projection  plane  of  first  end  point  of 
first  line  segment. 

P2(2) 

Real  array.  Contains  X  and  Y  coordinates 
on  projection  plane  of  second  and  point 
of  first  line  segment. 

Rl(2) 

Real  array.  Contains  X  and  Y  coordinates 
on  projection  plane  of  first  end  point  of 
second  line  segment. 

R2(2) 

Real  array.  Contains  X  and  Y  coordinates 
on  projection  plane  of  second  end  point 
of  second  line  segment. 

6.20  SUBROUTINE  PLPLN 


A 

Real  scalar.  The  fraction  of  a  contour 
line  that  a  single  step  represents. 

11 

Integer  scalar.  Starting  point  number 
for  the  current  side. 

12 

Integer  scalar.  Ending  point  number  for 
the  current  side. 

INDEX2 

Integer  scalar.  Variable  telling  how 
large  to  dimension  SEG  In  subroutine 
PNTPLT. 

I  PEN 

Integer  scalar.  Calcomp  pen  control 
variable. 

NSIDES 

Integer  scalar.  Number  of  points  for 
present  polygon. 

NUM 

Integer  scalar.  Total  number  of  points 
that  make  up  the  contour  of  a  plane. 

SEG(3,3000) 

Real  array.  Array  of  vectors  that  make 
up  the  sides  of  a  polygon. 

SUBROUTINE  PNTPLT 

I  FLAG 

Integer  scalar.  Flag  passed  back  from 
subroutines  HYDE  and  HIDE  telling  whether 
or  not  the  line  is  hidden. 
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INDEX2 

INUM 

I  PEN 

I  PLOT 

KK 

LFLAG 

NEWPEN 

NPTS 

P(3) 

PP(3,2) 

PPP (3) 

SEG ( 3 , I NDEX2 ) 

X 


Integer  scalar.  Input  argument  telling 
how  large  to  dimension  array  SEG.  Calcu¬ 
lated  In  Main  as:  (NSTEPS  for  that 
object  *  4)  ♦  1. 

Integer  scalar.  Represents  the  number  of 
objects  which  overlap  with  a  particular 
object. 

Integer  scalar.  Pen  value  sent  to  the 
Cal  comp  subroutine  PLOT.  Either  pen  up 
(3)  or  pen  down  (2). 

Integer  scalar.  Flag  sent  to  subroutine 
CLIP  to  plot  only  first  line  segment 
after  clipping  determined. 

Integer  scalar.  Overlapping  object 
number  that  overlaps  with  present  object, 
from  array  IE. 

Integer  scalar.  Flag  set  which  tells 
whether  or  not  last  point  was  hidden. 
LFLAG«1  means  point  was  hidden  and 
LFLAG-2  means  not  hidden. 

Integer  scalar.  Variable  contains  value 
of  last  pen  move.  The  number  Is  positive 
If  the  pen  move  was  performed  and  nega¬ 
tive  if  the  pen  move  was  clipped. 

Integer  scalar.  Number  of  points  to 
plot,  from  subroutines  PLPLN  and  PSE. 

Real  array.  Contains  coordinates  of  end 
point  of  line  returned  by  subroutine 
EXTEND. 

Real  array.  Contains  coordinates  of  end 
points  of  line  sent  to  subroutine  EXTEND. 

Real  array.  Coordinates  of  point 
currently  being  plotted,  converted  to 
viewpoint  reference  frame. 

Real  array.  Contains  coordinates  of 
points  to  be  plotted,  passed  through 
argument  list  by  subroutines  PSE  and 
PLPLN. 

Real  scalar.  Final  X  coordinate  of  point 
to  be  plotted. 


Ill 


XMAX  Real  scalar.  Constant  that  represents 

the  right  limit  of  the  plotting  region. 
Defined  In  data  card  16.0  In  Input  con¬ 
trol  file. 

XM1N  Real  scalar.  Constant  that  represents 

the  left  limit  of  the  plotting  region. 
Oeflned  In  data  card  16.0  In  Input  con¬ 
trol  file. 

XSAV  Real  scalar.  Represents  X  coordinate  of 

previous  point  plotted. 

Y  Real  scalar.  Final  Y  coordinate  of  point 

to  be  plotted. 

YHAX  Real  scalar.  Constant  that  represents 

upper  (top)  limit  of  the  plotting  region. 

YMIN  Real  scalar.  Constant  that  represents 

lower  (bottom)  limit  of  the  plotting 
region. 

YSAV  Real  scalar.  Represents  Y  coordinate  of 

previously  plotted  point. 


6.22  SUBROUTINE  POL YD 


Dl(810) 

INDEX(G) 

J 

NUM 

SUMD1 

SUMD2 


Real  array.  Equivalent  with  array  D, 
contains  directional  cosine  matrices  for 
the  polygons. 

Real  array.  Contains  constants  used  for 
transposing  data  elements  In  D  to  get 
correct  direction  In  matrix  D. 

Integer  scalar.  Calculated  variable  that 
converts  a  3-dlmensional  subscript  value 
to  a  single  dimensional  one. 

Integer  scalar.  Offset  used  to  put  data 
from  P  array  Into  SEGLP  array. 

A 

Real  scalar.  Magnitude  squared  of  the  x 
coordinate  axis  of  .the  polygon. 

A 

Real  scalar.  Magnitude  squared  of  the  y 
coordinate  axis  of  the  polygon. 
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6.23  SUBROUTINE  PREPLT 


IPEN 

NEWPEN 

X 

XLIMIT 

XMAX 

XMIN 

XOFF 

XSAV 

XTEMP 

Y 

YLIMIT 

YMAX 

YMIN 

YOFF 


Integer  scalar.  Variable  that  denotes 
Calcomp  pen  move  value  from  PNTPLT. 

Integer  scalar.  Saved  value  of  the 
previous  pen  (IPEN)  move. 

Real  scalar.  X  coordinate  of  first  end 
point  of  line  X,Y  -  XSAV,  YSAV. 

Real  scalar.  If  X  coordinate  outside 
valid  plotting  region  Is  equal  to  bound¬ 
ary  value  (left  or  right),  pen  Is  off. 

Real  scalar.  Constant  that  represents 
right  limit  of  valid  plotting  region. 
Defined  In  subroutine  PNTPLT. 

Real  scalar.  Constant  that  represents 
left  limit  of  valid  plotting  region. 
Defined  In  subroutine  PNTPLT. 

Logical  scalar.  Is  .FALSE.  If  X 
coordinate  In  plotting  region,  .TRUE.  If 
outside. 

Real  scalar.  X  coordinate  of  second  end 
point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  X  coordinate  of  the 
Intercept  point  at  YTEMP  between  X,  Y  and 
XSAV,  YSAV. 

Real  scalar.  Y  coordinate  of  the  first 
end  point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  If  Y  coordinate  outside 
valid  plotting  region,  is  equal  to  bound¬ 
ary  value  (top  or  bottom)  pen  Is  off. 

Real  scalar.  Constant  that  represents 
upper  (top)  limit  of  the  plotting 
region.  Defined  in  subroutine  PNTPLT. 

Real  scalar.  Constant  that  represents 
lower  (bottom)  limit  of  the  plotting 
region.  Defined  in  subroutine  PNTPLT. 

Logical  scalar.  Is  .FALSE.  If  Y 
coordinate  In  plotting  region,  .TRUE.  If 
outside. 
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YSAV 


Real  scalar.  V  coordinate  of  the  second 
end  point  of  tine  X,  Y,  and  XSAV,  YSAV. 


YTENP  Real  scalar.  Y  coordinate  (either  YM1N 

or  YMAX)  where  pen  went  out  of  the 
plotting  region. 


6.24  SUBROUTINE  PRJELR 


All 

Real  scalar, 
component  of 

Represents  the  aM 

To]. 

A12 

Real  scalar, 
component  of 

Represents  the  a12 
[a]. 

A22 

Real  scalar, 
component  of 

Represents  the  a22 
[a]. 

DD (3.3) 

Real  array, 
intermediate 
array  to  A'. 

Array  containing 
values  for  transposing  A 

ODD (3.3) 

Real  array, 
array. 

Contains  values  of  the  A1 

LAMDA1 

Real  scalar, 
circumscribe 

Eigenvalue  used  to 
ellipse  with  rectangle. 

LAMDA2 

Real  scalar, 
circumscribe 

Eigenvalue  used  to 
ellipse  with  rectangle. 

Ml 

Real  scalar, 
axis  vector. 

Value  representing  major 

M2 

Real  scalar, 
axis  vector. 

Value  representing  minor 

R(3.3) 

Real  array.  Contains  X,  Ywand  Z 
components  of  r?,  rj,  and  r 3. 

R2(2,3) 

Real  array, 
and  rf  on  the 

Contains  values  of  r7,  ri 
•  projection  plane. 

RX1 

Real  scalar. 

rx  value  for  rj. 

RX2 

Real  scalar. 

rx  value  for  r2. 

RY1 

Real  scalar. 

ry  value  for  r£. 

RY2 

Real  scalar. 

ry  value  for  r£. 
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S  ( 3 ) 

Real 

array. 

Contains 

matrix  multiple  of 

DVP 

and  SS. 

SS(3) 

Real 

array. 

Contains 

position  vectors 

for  individual  ellipsoid  minus  inertial 
reference  frame. 


6.25  SUBROUTINE  PRJPLY 


NPTS 


Integer  scalar.  Number  of  points  in  the 
polygon.  Taken  from  array  NPPP. 


PP 1  ( 3 ) 

PP2(3) 

PPP2(3) 


Real  array.  Vector 
origin  to  the  point 
inertial  coordinate 

Real  array.  Vector 
origin  to  the  point 
inertial  coordinate 


from  the  viewpoint 
(vertex)  in  the 
system. 

from  the  viewpoint 
(vertex)  in  the 
system. 


Real  array.  Vector  from  the  viewpoint 
origin  to  the  point  (vertex)  in  the 
viewpoint  coordinate  system. 


6.26  SUBROUTINE  PSE 


I  HALF 

Integer  scalar.  Flag  that  indicates 
which  half  of  the  ellipsoid  to  plot. 

I HALF  =  1 ,  semiellipsoid  with  X  >  0 
plotted,  IHALF=2,  semiel 1 ipsoid  with 

X  <  0  plotted. 

IN 

Integer  scalar.  Number  of  points  in  each 
quarter  contour  saved  in  array  XI. 

INDEX 

Integer  scalar.  Number  of  steps  plus 
one. 

INDEX2 

Integer  scalar.  Maximum  number  of  points 
any  complete  contour  can  have. 

IPEN 

Integer  scalar.  Calcomp  pen  control 
variable. 

KK 

Integer  scalar.  Variable  in  DO  60  loop 
which  runs  the  loop  backwards. 

LINE 

Integer  scalar.  Varaible  in  00  100  loop 
which  runs  the  loop  backwards. 
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NPT 

NPTS 

SEG(3,INDEX2) 

XI (3, INDEX, INDEX) 


Integer  scalar.  Number  of  points  in 
semiellipsoid  to  transfer  to  array  SEG. 

Integer  scalar.  Number  of  points  in 
semiellipsoid. 

Real  array.  Array  containing  a  complete 
contour. 

Real  array.  Semi  ellipsoid  contour  array. 


6.27  SUBROUTINE  ROT 


A(M,3) 


Real  array.  Output  of  ROT,  rotation 
matrix  to  be  computed. 


C 


Real  scalar.  Cosine  of  angle  of 
rotation. 


L 

M 

S 

TH 


Integer  scalar.  Variable  telling  which 
axis  to  rotate  about  (X=l,  Y=2,  Z=3). 

Integer  scalar.  First  subscript  value 
for  array  A. 

Real  scalar.  Sine  of  angle  of  rotation. 

Real  scalar.  Angle  of  rotation 
(radians). 


6.28  SUBROUTINES  SOLVA,  SOLVR 

See  Appendix  B,  Discussion  of  Equations  used  by  PRJELR.  The 
variables  correspond  to  the  equations  found  in  this  appendix. 


6.29  SUBROUTINE  TITLE 


10(10,20)  Integer  array.  Contains  text  for  title 

frame. 

NFRAME  Integer  scalar.  Number  of  title  frames 

to  plot. 
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NUNC 


Integer  scalar.  Constant  that  represents 
number  of  lines  In  each  title  frame. 


6.30  SUBROUTINE  TPOINT 


Real  scalar.  X  coordinate  of  where  to 
start  plotting  title  frame. 

Real  scalar.  Height  of  characters  In 
title  frame. 

Real  scalar.  Y  coordinate  of  where  to 
start  plotting  title  frame. 


NPTS1 


PP 1 ( 3 ) 


PP2(3) 


Integer  scalar.  Denotes  polygon  on 
projection  plane. 

Integer  scalar.  Flag  passed  back  to  I N= 1 
caller,  IN»1  means  point  Inside  polygon, 
means  outside  the  polygon. 

Integer  scalar.  Number  of  points  for 
polygon  I. 

Real  array.  Contains  position  vector  for 
polygon  I. 

Real  array.  Point  on  the  projection 
plane  to  be  tested. 

Real  array.  Contains  differences  between 
PP1  and  PP2. 


SIGN2 


Real  scalar.  Test  variable  created  like 
array  SIGN. 


6.31  SUBROUTINE  TRANS 1 


DD ( 3 , 3 ) 


Real  array.  Contains  matrix  multiple  of 
OVP  and  parts  of  D. 

Real  array.  Output  of  TRANS1,  it  is 
position  vector  of  an  ellipsoid  trans¬ 
formed  to  the  viewpoint  reference 
frame. 

Real  array.  Input  to  TRANS1,  contains 
position  vector  for  surface  points  of  an 
ellipsoid  In  the  segment  coordinate 
system. 
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R2(3) 

SEGLP2 ( 3 ) 

6.32  FUNCTION  XINTCP 

PFACTR 

X 

XI 

XINTCP 

XSAV 

Y 

Y1 

Y2 

YSAV 

YTEMP 

6.33  SUBROUTINE  XYZ 

See  Appendix  B, 
Variables  in  XYZ 
appendix. 


Real  array.  Contains  the  position  vector 
of  an  ellipsoid  in  the  viewpoint  coordi¬ 
nate  system. 

Real  array.  Contains  the  location  of  the 
segment  c.g.  in  the  viewpoint  coordinate 
system. 


Real  scalar.  The  slope  of  the  line 
X,  Y  -  XSAV,  YSAV. 

Real  scalar.  X  coordinate  of  first  end 
point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  Difference  between  X  and 
XSAV. 


Real  scalar.  X  coordinate  value  of  YTEMP 
in  line  X,  Y  -  XSAV,  YSAV. 


Real  scalar.  X  coordinate  of  second  end 
point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  Y  coordinate  of  first  end 
point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  Difference  between  Y  and 
YSAV. 


Real  scalar.  Difference  between  YTEMP 
and  YSAV. 


Real  scalar.  Y  coordinate  of  second  end 
point  of  line  X,  Y  -  XSAV,  YSAV. 

Real  scalar.  Y  coordinate  value  (between 
Y  and  YSAV)  for  which  the  caller  wants 
the  X  coordinate. 


'Hidden  Line  Problem  Between  Two  Ellipsoids." 
correspond  directly  to  the  equations  in  this 
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6.34  FUNCTION  YINTCP 


See  Appendix  B,  "Hidden  Line  Problem  Between  Two  Ellipsoids." 
Variables  in  YZ  and  Z  correspond  directly  to  the  equations  in  this 
appendix. 
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7.0  SUBROUTINE,  COMMON  BLOCK,  AND  VARIABLE  CROSS-REFERENCE  CHARTS 


This  section  contains  three  cross-reference  charts:  subprograms 
called  by  other  subprograms,  common  blocks  used  by  subprograms,  and 
variables  contained  within  each  common  block. 


I 


7.1  Subprogram  Cross-Reference  Chart 


I  CALLING  !  H  0rC  C  C  0  0  0  0  E  E  5  H  H  I  L  ft  N  0  P  P  P  P  *  P  P  A  S  5  T  T  T  X  X  Y  Y  2! 

i  rout  ine  miionoomniyiinvutommoiMiui  .• 

f  jiiinotttcitnddpetreptlejjetlltoanzn  ! 

J  I  HLM  S  TVPEDIEUG  ARLPYPEP  VVUNT  T  ! 

!  !DRS  PSNC  T  1  MLNLDLLL  ARENSCC  ! 

!  !  1  E  R  N  0  M  NEATTRY  TIP  P  ! 


! 

! 

E 

C 

T 

P 

!  CALLED 

! 

!  ROUTINE 

1 

FREQ 

J 

* 

1 

Ibuildie 

1 

""i 

j 

i 

!  CLIP 

I 

i 

3 

!  CONVREC 

1 

i 

i 

•CROSS 

2 

2 

2 

!  DET 

A 

i 

A 

!  DOT 

9 

i 

i 

1 

2 

2 

2 

1 

!  DOTT 

3 

2 

1 

1 

1 

•DRCYPR 

1 

i 

1 

! ELIPSN 

1 

1 

i 

! EXTEND 

1 

1 

1 

! GENDCM 

1 

J 

1 

•HIDE 

2 

I 

1 

1 

!  HYDE 

2 

1 

1 

1 

JINPUT 

1 

1 

i 

• LSEGINT 

1 

1 

1 

SMAT 

16 

2 

2 

2 

2 

A 

3  1 

2 

•NFRANE 

2 

i 

1 

1 

! OVERLAP 

3 

i 

3 

! PLPLN 

1 

i 

1 

! PNTPLT 

2 

i 

1 

1 

! POLYO 

1 

i 

1 

! PR  EPLT 

3 

i 

3 

! PRJ  ELR 

1 

j 

1 

JPRJPLY 

1 

i 

1 

!  PSE 

2 

i 

2 

!  ROT 

3 

i 

3 

iSOLVA 

1 

i 

1 

! SOLVR 

3 

i 

3 

! TITLE 

1 

i 

1 

! TPOINT 

2 

i 

1 

1 

! TRANS1 

3 

i 

1 

2 

!  XINTC  P 

3 

» 

2 

1 

!  XYZ 

1 

i 

1 

!  YINTCP 

3 

i 

2 

1 

!  YZ 

1 

i 

1 

!Z 

1 

i 

i 

1 
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7.2  Common  Block  Cross-Reference  Chart 


!  COMMON 

m  mm  m 

!  A 

C 

■  m 

D 

C 

I 

p 

m  m 

P 

R 

V! 

!  3LCCK 

!  T 

0 

3 

L 

N 

L 

0 

e 

I! 

!  H 

N 

U 

L 

T 

T 

L 

M 

C  1 

i 

• 

c 

G 

I 

E 

T 

Y 

C 

w! 

i 

c 

P 

R 

G 

V 

P! 

; 

T 

S 

S 

0 

E 

i 

• 

i 

• 

c 

'1 

i 

!  CALLING 

i 

j 

!  ROUTINE 

i 

i 

; 

j 

!  MAIN 

!  * 

* 

* 

* 

* 

* 

* ! 

! 3UILDI6 

; 

* 

* 

* 

* 

j 

!  CLIP 

i 

i 

! CONVPEC 

!  * 

* 

* 

* 

* 

t 

! CPOSS 

i 

i 

!  OET 

j 

t 

!  DOT 

i 

i 

!  DOTT 

i 

i 

! DRCYPR 

1 

! SLIPSN 

j 

★ 

i 

! EXTEND 

i 

* 

* 

* 

t 

! SENDCM 

i 

i 

!  HIDE 

i 

* 

* 

i 

!  H  YOE 

• 

* 

i 

! INPUT 

j  * 

* 

* 

* 

* 

* 

* 

* 

♦  J 

! LSEGINT 

; 

i 

!  MAT 

j 

i 

• NFKAME 

i 

! OVERLAP 

i 

* 

; 

! PLPLN 

i 

* 

★ 

* 

* 

* 

; 

! PNTPLT 

• 

* 

★ 

ft 

; 

! POLYD 

i 

* 

* 

i 

! PPEPLT 

I 

i 

! PRJELR 

i 

* 

* 

i 

! PR  J  PLY 

• 

* 

* 

i 

!  PSE 

• 

i 

!  ROT 

• 

i 

! SOLVA 

• 

i 

! 5CLVR 

i 

i 

! TITLE 

; 

* 

i 

! TPOINT 

I 

* 

i 

! TRANSl 

; 

★ 

*  J 

•  XINTCP 

j 

i 

!  XYZ 

; 

i 

!  YINTCP 

j 

i 

!  Y  Z 

j 

j 

!  Z 

i 

i 

; 

i 
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7.3  Variable  Cross-Reference  Chart 


!  COMMON 

!  A 

m 

V 

2 

E 

i 

p 

P 

R 

V! 

!  SLOCK 

!  T 

0 

5 

L 

N 

L 

0 

c 

I ! 

j  a 

N 

U 

L 

T 

T 

L 

M 

E! 

1 

• 

E 

G 

I 

z 

T 

Y 

0 

W! 

1 

a 

C 

P 

a 

5 

V 

P! 

1 

• 

T 

S 

s 

0 

E 

l 

• 

E 

N 

!  VARIABLE 

; 

j 

!  A 

j 

★ 

!  30RS 

j 

★ 

!  CONVEC 

i 

• 

* 

!D 

i 

★ 

12EVFLG 

i 

• 

* 

!  DV? 

; 

* 

!  D  V  3  0 

i 

*  { 

•  I C  OL  0  3 

i 

* 

! IOERUG 

i 

* 

!  IE 

i 

* 

!  I  EL° 

j 

★ 

•IFLAG 

j 

* 

!  INT 

j 

* 

»IRE*OV 

; 

* 

!IV=> 

i 

• 

* ; 

!  i^°L 

; 

* 

!  N  I  £ 

i 

• 

★ 

!  N  I  S  o 

i 

• 

* 

!  N° 

i 

• 

* 

!N PLANE 

; 

* 

JNPPP 

i 

• 

* 

!  m  p  a  e  r 

j 

* 

!  N  S  -  G 

; 

* 

INSTEPS 

j 

* 

JOFLINE 

j 

★ 

JCFSETX 

i 

* 

!  0  c  S  £  T  Y 

j 

* 

!o\l:ne 

j 

* 

1  p 

; 

★ 

!  PL 

J  * 

!  °CS 

j 

* 

!  PC 

j 

★ 

!  RA 

; 

* 

!  SEGLP 

; 

* 

i SFACTR 

j 

* 

!  SIGN 

; 

* 

' 

ITERM 

j 

* 

!  T  I  v  £ 

j 

* 

!VP 

j 

* 

!  VPO 

; 

*  j 

!  VP  2 

i 

• 

* ; 

!  ZTIM£ 

j 

i 

* 
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8.0  VIEW  PROGRAM  SOURCE  LISTING 


PROGRAM  VIEW  1 

COMMON/PLTT/SFACTR,INT,TIM£, I C OLOR < 9 1 > ,0F3ETX ,0F S STY, ZTIME  2 

COMMON/ INTERS/  N I E < 90 ) , I E (90/ 90)  J 

COMMON /£LLI®SE/NSTEPJ(S0)/!£LP/A(3.3/30)/SEGLP(3»90)/VP(3)»  4 

»D< J/3,90>,0VP(3,3>,RA<3>,NSEG  5 

C9MK0N/P0LYG0N/NPLANE, I FLAG,NPPP<93> ,P3( 3, 4,60) ,® < 3, 4,60) ,  6 

♦CONVEC  C:,<.,VO>,POS(2,9C),S:SN(90>  7 

t jMMOS/Arj/PL(17,SO)  3 

COM < ON/VIE WP/VP0(3),9VP0(3, 3) ,:v°,vp2< 3)  9 

C0MM0N/3uJG/I0e3UG(oQ),NISG,9£VFLG, ONLINE, TERM, B3RS*0FLINE  10 

DIMENSION  WORK  ( 1 3010)  11 

DIMENSION  19(10)  12 

INTEGER  DEVFL  G,ONLI  < E , TERM, EDR S, 3 FL INE  13 

DOUBLE  PRECISION  ZTIM£,CTIM£,STIM£,DTIME,ETIME  14 

C  15 

C  PROGRAM  VIEW  VERSION  1.1  MARCH  9,1933  15 

C  WRITTEN  6Y  SRL  IN  SUPPORT  OF  THE  17 

C  MATHEMATICS  AND  ANALYSIS  3RANCH  OF  AMRL  IS 

C  AT  «i  P  A  F  9 «  19 

C  20 

C  STORED  ON  T APE  1  THAT  IS  OUTPUT  FROM  THE  ATOM  VERSION  V5D.  21 

C  22 

C  THIS  PROGRAM  USES  CONTOUR  LINES  TO  REPRESENT  THE  3-D  PROPERTIES  OF  23 

C  THE  DATA  ON  TAPE1.  THE  CONTOUR  LINES  ARE  PLOTTED  ON  PAPER  THAT  24 

C  REPRESENTS  THE  PROJECTION  PLANE.  THE  POINTS  THAT  COMPOSE  A  CONTOUR  25 

C  LINE  IN  3-SPACE  ARE  PROJECTED  THROUGH  A  POINT  ON  TO  THE  PROJECTION  25 

C  PLANE.  27 

C  28 

C  CURRENTLY  THERE  ARE  TWO  CLASSES  OF  OBJECTS  THAT  ARE  PLOTTED  USING  29 

C  CONTOUR  LINES.  30 

C  CLASS  1  -  ELLIPSOIDS  31 

C  ELLIPSOIDS  ARE  USED  TO  REPRESENT  93DY  SEGMENTS.  THIS  PROGRAM  ALLOWS  32 

C  ELLIPSOIDS  TO  BE  IMBEDDED  IN  OTHER  ELLIPSOIDS.  23 

C  34 

C  CLASS  2  -  CONVEX  POLYGONS  35 

C  CONVEX  POLYGONS  ARE  USED  TO  REPRESENT  03JECTS  THAT  CAN  BE  DEFINED  35 

C  3Y  A  SET  OF  PLANES.  ALL  POLYGONS  DEFINED  BY  THE  INPUT  MUST  BE  37 

C  CONVEX  POLYGONS;  CONCAVE  POLYGONS  CAN  3E  OBTAINED  USING  A  3 3 

C  COMBINATION  OF  CONVEX  POLYGONS.  39 

C  40 

C  THE  HIDDEN  LINE  “CJTINES  CHECK  FOR  POINTS  HIDDEN  9Y  41 

C  ELLIPSOIDS  OR  POLYGONS.  THESE  ROUTINES  "UST  CHECK  42 

C  FOR  ANY  POSSIBLE  03JECT  THAT  MAY  BE  BLOCKING  THE  43 

C  CURRENT  POINT  AS  SEEN  FROM  THE  VIcW’OINT.  44 

C  IN  ORDER  TO  ELIMINATE  CHECKING-  ALL  OBJECTS  FOR  45 

C  EACH  POINT,  SUBROUTINES  ARE  INCLUDED  IN  THE  46 

C  VIEW  PROGRAM  THAT  DETECT  03JECT  OVERLAP  ON  THE  47 

C  PROJECTION  PLANE.  BEFORE  THE  PLOTTING  PHASE  OF  THE  43 

C  VIEW  PROGRAM,  OBJECTS  WHICH  OVERLAP  EACH  OTHER  ON  49 

C  THE  PROJECTION  PLANE  ARE  RECORDED  IN  THE  IE  ARRAY.  50 

C  DURING  THE  PLOTTING  PHASE  OF  THE  VIEW  PROGRAM,  THE  51 

C  IE  ARRAY  IS  REORDERED  TO  DECREASE  THE  SEARCH  TI"c  52 

C  FOR  OBJECTS  THAT  MAY  BE  BLOCKING  THE  CURRENT  ®OINT  53 

C  BEING  PLOTTED.  THE  ASSUMPTION  USED  HERE  IS  -  IF  AN  54 

C  OBJECT  BLOCKED  THE  PREVIOUS  POINT  ON  A  CONTOUR  LINE  55 

C  THEN  THAT  OBJECT  PROBABLY  BLOCKS  THE  NEXT  POINT  56 

C  ON  THE  CONTOUR  LINE.  57 

C  59 

LUPLOT=  3  59 

0NLINE=1  60 

TERM*2  51 

NF  =0  52 
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0FL:N£*3  43 

3DRS*4  64 

RE  A9 ( 5, 1 30 )  OEVFLG  65 

VRITEC6/130)  OEVFLG  66 

130  FCRMAT(II)  67 

C  68 

I'  (DEVFLi.EO.1 .03.0EVFLG.EG. 3)  CALL  PLOT S (0/ 3, LUPLDT)  69 

If  (OEVFLG. £4.2.09. OEVFLG. E3.»)  CALL  PLOT $( 0,0, LU»LOT>  70 

C  71 

I  FLAG  *  C  72 

CALL  TITLE  73 

READ(5,2GC)  (ID(I),I*1,10)  74 

WAIT£(6/2J0)  (ID(I)/ 1*1/10)  75 

200  F0RMAT(10A4)  76 

3EADC5/150)  STINE/OTIME/ETIME  77 

150  FORMAT (301 0.0)  78 

CTI'IE  =  STI>;E-DTI»'E  79 

ITIMCsCTir,E*1 003000. 00  80 

CT::iE=:riKE/ioooooo.DO  hi 

R  £ AD ( 5/ 1 25)  IDE3UG  82 

v»KITE(o/125)  IDE8UG  63 

125  FCRMAT(SOII)  84 

100  CONTINUE  85 

I CLA 3  =  I  FLAG  ♦  1  66 

C7IrtE=CTIME+DTIME  87 

ITI"E=CTIM;*1000 jOO.DO  38 

CTI*c=ITlrS/1000000.00  89 

IF(CTIME.GT.ETHE)  CALL  PLOT ( 0. / 0. / 99?)  »0 

tF(CTIME.GT.ETIME)  STOP  91 

CALL  INPUT ( CTI WE )  92 

IF(IFLAG.sa.lO)  GO  TO  100  93 

IFCDEVFLG.  EG.  OFLINE.  OR.  OEVFLG. EQ.aDRS)  CALL  NELP EN ( I COLCP ( 91 )  )  94 

XTIM£  =  ZTI*!£«1000.0  95 

N  F  =  N  f  ♦  1  96 

NA I T£ ( 6/ 1 030) N  F  97 

1000  F  03  HAT  ( *  MAIN  -  PROCESSING  FRAME  #'/I4)  98 

CALL  PLOT(0./C./-3)  99 

CALL  SYM39LC. 5/1 0.0/. 335/ID/C. ,35)  100 

CALL  5Yf'i?0L(.5/9.0/.335/*TIME(MSEC),/0./13)  101 

CALL  NU'53E3(  3.35,9.0/.  335/XTIME/0./-1)  1C2 

CALL  ',AT(?VP0/O(1/1/IV°),DV-’/3/I/3»3/2/3)  103 

CALL  DOT(0(1,1,IVP),V?0,V»,3,1,3>  104 

DO  10  K=1/3  105 

10  VP(<)*VP(0+SSGLPOC,IVP)  106 

CALL  MAT(DVP/VP/VP2/3/I/1/3/3/3)  107 

CALL  CON VR EC  108 

CALL  =RJDLY  109 

CALL  °RJEL3  110 

CALL  PCLYO  111 

CALL  5JILDIE  112 

IFC0£8UC(1).£3.1)  «R  I T  E  (  6/ 3  50 )  (M  E  ( I)  /  1  =  1 ,90)  113 

IF(ISE3U5(2).E3.1)  ViRITc(6/35J)  ( ( I E  ( I,  J  )  / 1  =  1  /  «  0)  /  J  =1  /  90)  114 

350  F0RMAT(Z73(JO(1X,I2)/))  115 

DO  30  IK=1 /NSEG  116 

IFOEVFLG.ES. OFLINE. OR.  OF V  =  L3 . E 0 .3DR S)  CALL  NEWPEN  (  I  COLOR ( I K ) )  117 

:el°=ik  ns 

IN0EX  =  N$TcPS(I>()+1  119 

IiO£X2  =  4*INDEX-3  120 

1X1=1  121 

IIN=3*INDEX»INDEX+IX1  122 

I5E3=IIN ♦INDEX  '  123 

CALL  iLI°SN(INDEX/W0®<(IX1),WD°<(IIN))  124 

CALL  PSE(W0R<(IX1)/WORX<:iN)/'WO3X(ISEG)/INDEX/INDEX2/1)  125 
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CALL  PSE(W2PX(!X1)/WCRK(!IN)/W0RX{ISEG)/IN0£X*ISDeX2/2)  126 

JO  CONTINUE  127 

CALL  PLPLM(W0R</INDEX2)  128 

IFOEVFLG.E3.30AS>  CALL  NFRA'IE  129 

IFOEVFLC.£a.OFLIN£.OR.DEVFLG.E3.3\LINE>  CALL  PLOT  (12.,0.0/-3>  130 

IFOEVFLG.e3.TEPi*)  CALL  PLOT  (0./0./-3)  131 

GO  TO  100  132 

SO  133 


) 
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SUBROUTINE  3UILDIE  174 

C  13$ 

C  OSCE  THIS  SUBROUTINE  IS  CALLED  ALL  OBJECTS  ARE  REPRESENTED  3Y  136 

C  POLYGONS  PROJECTED  ON  THE  VIEWPOINT  PROJECTION  PLANE.  13? 

C  THIS  SUBROUTINE  wIlL  BUILD  THE  It  AND  NIE  ARRAYS.  133 

C  NS£<<)  REPRESENTS  TmE  SEMPER  OF  ENTRIES  IN  THE  I£(I,K)  ARRAY  FOR  139 

C  OBJECT  X .  1U 

C  THE  I  £ (  I  / < )  ARRAY  CONTAINS  OOJFCT  SJNjERS  FO»  OBJECTS  THAT  OVERLAP  U1 

C  I,N  THE  PROJECTION  PLANE.  1 42 

C  FOR  EXAMPLE/  ltd  .2)  ¥I5HT  CONTAIN  A  3  WHICH  MEANS  03JECT  3  US 

C  OVERLAPS  WITH  OBJECT  2.  144 

C  OMMON/POLTSON/NPLANE, IFLAG,NP«»P<9C),P0<3,4,6C>,PC3,4,63>,  14$ 

•C JN VS C< 2,4,90) ,POS<  2,90), SI 3N<  90)  146 

C9M*lON/£LLIRSE/NjTEPS(90),lELP,A<3,3,30),SE3LP(3, 90), VP (3),  147 

*D(3,3,90),DVP(3,3>,*4(3),NSEG  143 

COMHON/INTERS/  NIt(90),IE ( 90, 9  J)  149 

COMMON/ REPO  VE/NRHEM/  I  RE, *13  VC  30)  150 

DO  5  I»1,9C  151 

NIE(I)  *  0  152 

00  5  X*1,90  153 

5  IECI/K)  «  0  154 

I F ( NS  EG  .E3.  0)  33  TO  60  155 

03  55  :=1/N3£G  156 

IMNISU)  .NS.  T)  G"  TO  10  157 

M  £  ( I )  =  1  158 

I E ( 1 / I )  *  I  15° 

10  I  03 J  =  1  +  1  160 

IF CI.E0.N5S5)  30  TO  31  151 

DO  30  K*IO?J,NS£G  162 

CALL  OVERLARCI/X/HFLAG)  163 

IF (MFLA3  .El.  0)  SO  TO  30  164 

C  165 

C  YES,  THERE  IS  OVERLAP  BETWEEN  I  AND  K  166 

C  167 

IF  ( S'  I H  ( K  >  .NE.  0)  50  T3  20  168 

N I  c  (  K  )  =  1  169 

IE(1,K)  =  X  170 

20  NIE(K)  =  MECX)  ♦  1  171 

Man  =  n i e < x )  ♦  i  172 

IE(NIc(I)/I)  =  <  173 

I£(N!£(K),K)  =  I  174 

30  CONTINUE  175 

!1  CONTINUE  176 

IF (N°LANE.E3.0)  SO  TO  55  1 7? 

1 03 J  =  MPLANE  ♦  20  17- 

DO  50  r.  =  J1,I03J  179 

DO  2CC  LT  =  1,N?PEM  UO 

IF(X“3C.E'i.IRE¥.CV(LT))  30  TO  50  131 

200  CCNTI.Ut  182 

CALL  OVERLAPCI,X,NFLAG)  183 

IFCNFLA5  .EQ.  0)  GO  TO  50  134 

C  165 

C  YES,  THERE  IS  OVERLAP  3ETWEEK  I  AND  X  186 

C  137 

NIECK)  =  ME  (  K )  ♦  1  130 

NIE(I)  =  NIE(I)  ♦  1  139 

:e<sieci>,i)  =  <  190 

IE(NIE(<),0  =  I  191 

50  CONTINUE  .  192 

55  CONTINUE  193 

C  194 

C  NOW  CHE CX  PLANE  AGAINST  PLANE  195 
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60  1HKPLASE  .LI.  1)  RETURN  jqf 

NPLANEtNPLANEPIV  1(H 

<PLANE»NPLAN«pJ9  19; 

95  IOC  I  * J1 #NPL  *NE  100 

30  ICJ  LT«1#NPRSN 

If  <I-2».E3.lRE»i:V<LT>>  50  T  3  10U  20' 

300  CONTINUE  2ol 

IOEJ  ■  I  ♦  1  jr,* 

93  *03  K»I03J/Kf»LANE  ?05 

CALL  0VcRL4A(I*k/NfLAG)  ?)6 

If (  If LA 5  . £ 3.  0)  50  TO  403  207 

c  :Di 

C  T:5»  THERE  IS  OVERLAP  2ETWEEN  I  AND  K  2 09 

c  213 

N  1 1  (  < )  *  NIE(K)  ♦  1  >n 

Mill)  *  NIE(I>  ♦  1 

IE(NIKI)sI)  ■  K  jjj 

:£(NI£(K)/<)  »  I  »u 

400  CONTINUE  ->15 

103  CONTINUE  il* 

RET  URN  2 1 7 

END  218 
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iun«ouf  ;ns  :lip(«/ Y*xiAv*YiAv»xNiN*xMAx*y*iN* ypui* ip(n*i»lot) 


t  . . . . . 

c 

c  t •« i s  su:»ouiine  cl:»s  ^lotting  off  *otm  ends  of  ime  ttuour  ohum 

c 

t  . . . . . 

c 

logical  AOff#YO»f 

OATA  LC  AlL / J/ 

c 

C  OETCM*lNt  If  X  AND/OR  T  CLIPPING  AND  I  f  Off  TOf/'j  )ITO*.LI  F», 

C  0*  <»I3«I  Of  PL  3  T  T  t  * 

C 

KOMi.MLSE, 
rof *  » . f  als  e  . 

If  (x.LT  .X*IN.C».X. it.  X MAX)  X3FF«.T>U!. 

If  (Y.LT.VpIN.OX.Y.GT.YHAA)  YOff-.TRuE. 

If  <X.ir.X*I%)  XLI»’T«X»IN 
If  (X.3T.XMAX)  <l!,*IT»X«UX 
If  (Y.LT.Y*IN>  YLl»IT«Y4IN 
If  (Y.GT.  Y*AX>  YLIP2T*VNAX 
C 

C  IF  PREVIOUS  CALL  YO  C  L I  •*  WAS  »  PE  H  UP  TO  1ST  POINT  IN  St  3ME  NT/ 

C  INTERPOLATE  USING  i:«  AND  SAV'D  C03»5.‘S/  MOVE  PEN*  RESET  SAVE 
£  VALUES  fOP  X ♦  Y  o0HTi*  and  CONTINUE 

c 

If  (LCALL.EJ.O)  30  TO  10 

if  <x:ff)  yltenp«yikt:p(a/Y/asav*yspv**lip;t) 

If  (  .  NOT  .  X  Of  f  )  YL  T  £  “P*  YL  I  *f  I  T 

If  (  V  0  F  f )  XLTEMP’x INTCP{X,Y*X3AV*YsAV*YLI*IT) 

If  (.NOT. TOFF)  XLTE'<PjXLINI  T 
I*  <ALTEAB.LT.X*lr,)  XL  T  £  *«* »  X  *  I N 
If  <XLTEMP.GT.X«'AX)  XLTEN"«XMAX 
If  (VLTENP.LT. YMM  YlTE*p*YPIN 
If  (YLTEMP.OT.Y-iX)  YLTE  MP.y.-AX 
CALL  PLOKXLTfPB.YLTi-P*  S> 

XifcV=XLTEMP 
YS‘  V  =  YLT£r.° 

LC ALL-0 
C 

C  IF  1ST  POINT  Of  SEGMENT  AND  P£N  JP,  SAVE  THESE  COORD. 'S*  SET 
C  FLA3*  ANC  EXIT 
C 

10  CONTINUE 

If  <!?LOT.N-.1.0».!?£N.NE.T)  30  TO  23 
x  L  S  A  v  *  X 
y  L  *  *  V  8  Y 
LCALL^I 
I PLC  T  *3 
RETURN 
C 

C  03  * l  rant  TC  PLOT? 

C 

20  CONTINUE 

IF  (irLOT.NE.1)  30  TO  30 
C 

C  DETERMINE  X  AND  Y  COORDINATE:  TO  PLOT  TO 

c 

IF  (  X  0  F  F  )  YTCM?8YINTCP(X/Y*XSAV/YSAV*XLI'»IT> 

IF  (.NOT.XCff)  YTENP*YL!VIT 

Ic  (YOFf)  XTEMP8XINTCP(X,Y*XSAV» YSAV*YLI*fIT) 

If  (.NOT. YOFF)  XT'MOrXLI»IT 


21  * 
no 

f!  1 
:?? 
:?  i 

.*.** 
*  i  • 

2?S 

:;t 
2 : ! 

:  jo 

?ji 

2  J  J 
?H 
?  5  S 
?  54 
?.»7 
:  j  j 
2  3* 
?43 
?;i 
21  ? 
24  3 

244 

245 

24  6 
21  7 
245 
2  *  •* 
:so 
251 
25? 
?  5  3 
251 

25  5 
256 

255 
>5  / 

26  3 
?i1 
262 
263 
?64 
?65 
266 
26  7 
?6S 
26« 
?  7  3 
271 
?7? 

273 

274 

275 

276 

277 
2  7  a 
275 
2  ?  3 


c 

c 

c 

( 

S3 


C*U. 

(  }‘.f  t’.Ui 
I 3T 
«(IU*S 
l«9 


?»* 

t’i 

?!i 

Ut 

S4* 


su"*outin'  convrsc  ?90 

t  251 

c  this  subroutine  converts  rectangles  is  the  ah  29: 

C  SIMULATION  FORMAT  TO  PClT/ONi  IN  THE  VX:m  PLOTTING  FORNAT.  293 

£  294 

C)*V3*,/ah/»l(12,!3)  2  95 

JIXfiSIJi.  fc(!)  2  9 i 

9!aE'<in<  9>(!)  2  97 

CO'-CV/ POL  VGON/NP..ANE,  If  LAE./NPF  ’<  )0)  ,P  0<  X,4,S>0),3<  3,4,63),  2  93 

•  cosvf  c  /  90)  /®ds(  s !  j’a  50)  ??? 

£0"* 1N/C0NECT /  NP/NPL ( J,S,fcO)  3C0 

COrHON/COUL/IOCBUOISO),  *  I G 3, J : VF LG, OSL I \ E, T E R“, ?5 R 5, OFL I NE  301 

C  J  <''3s/-:u:»SE/NSTt?S<9J)»i:LP»  4(3, 3/  30),:EGlP(3,?0),V»(3>,  302 

•0<!,:,?L>,5VP<3,3>,pA(3),NSEC  30  3 

INTEGER  OEVflii, ONLINE, TERri, 30*5, 3FLINE  304 

IF(NPLANE.EJ.O)  RETURN  305 

IF  C10E?v,CU>.N£.0)  lpITE(6,50>  30s 

50  F0R“AT(lHl,' PLANS  I  Sf  OR  7  A  T  I  ON  '  /  /  ,  1  H  ,17(1*0)  307 

00  1  30  JO/NPLANE  303 

IF(J.GT.NP)  GO  TO  15  30? 

IFtlFLAG.NE.I)  GO  TO  15  310 

DDO*0£T(PL<1,J),PL<2,J),PLC3,J),PL<(5,J  >  ,PL<9, J)#»L(15.J)*  311 

PL(13/J)/Pl(l4/J),°L(15/J))  312 

0X(1)  r  dET£PL(4,J>,PL<2,J),pl(3,J),PL<11,J),PL(?,J>,PL(10,J),  313 

Pl(1  j»J),PL(14,J),Pw(15,J))  314 

OX (2)  =  DET(Pl(1/J)/PL(4,j),PL(3,J)/PL(<3,J)/PL(11/J)/PLC1Q/J)/  315 

PL<13,J),PL<16,J),»L<15,J)>  516 

D  X  <  3 )  =  DET(PL(1/J)/PL(2/J)/PL{4,J),PL{8/J)/3L(9/J)/0L{11/J)/  317 

PL(13/J)/?l(14/J)/?L(14/J))  318 

00  10  1=1,3  31? 

ts.x^dxc  : wood  320 

PjCI/1/J)sT?h?  321 

PO(I/2/J)=pl(IP7/J)«PL(12/J)*Tf*“  322 

P0(I/3/J)iPL<I+12/J)*PL(17/J)*p.(I/2/J)  323 

10  ?3(;/A,j)*pl(Ip12/J>*PL{17/J>*T-v?  324 

15  CONTINUE  325 

1  3  3  =  "  p  _ ( 1 /I / J  )  326 

IF(ISG.EO.O)  ISC*NISG  327 

IF(IOIP'Jj(4).M.O)  7E  <  6/200 )  1  33  323 

200  F0RS4T  (IX, ' I  5  5  =  ' / 1 2i  329 

00  20  l=1/4  353 

CALL  O0T(0(1/1/I'iG)/PC{1/L/J)/°/3/1/3)  731 

00  20  I  - 1  /  .7  332 

p<:,l,j)*r<i)P3sgl’(:,:»3)  333 

20  CONTINUE  334 

if(j.l:.n»)np?pu*30)=4  335 

I  F  <  IDEOL'.U)  .  SE.O)  WRIT  i(. 6 ,3000)  336 

3000  F3RNATC1 X, 30(1 H*) >  337 

IF(IOE8^G(4).NE.O)  »R  I  TE{ 6, 2000)  J  333 

2000  f:R1AT(3X/ ' PLANE  NUwc:R  =  ’,13)  33? 

JJ=J  340 

N»POP  =  .NPPP(J*30)  341 

IFdDEOUSUKNE.O)  URITE(6/1200)((P(I/X/JJ)/I=1/3)/X=1/NDPPP)  342 

1000  FC«MAT(3X/F7.2,3X/F7.2/3X/F7.2)  343 

100  CONTINUE  344 

RETURN  345 

END  346 
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SUBROUTINE  £  POS  S  ( k,  !l#C  )  347 

C  COMPUTES  VECTOR  CROSS  PRODUCT  C  =  »<3  34* 

DIMENSION  »(3)/T(3>/C<3)  34V 

C  (1  >*»(2>»3<  35D 

Ctl)««m>  J(1>-A(1W(5>  351 

CC3)*A(1)»2C?>“*(2)»311>  352 

RETURN  353 

E*D  354 
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FUNCTION  DET<A11/A12/A13/A21/A22/A23/A31/A32/A23) 
DET=A11*(A22»A33-A23*A32)-A12*(A?1*A33-A23*A31> 
1+A13*(A21*A32-A22»A31) 

PE  TUAN 
END 


355 

356 

357 
353 
359 


SUBROUTINE  OCT<»/3,C,N/M/L> 

C  REtf  03  OS/31/73  361 

C  PERFORMS  MATRIX  MULT I°L I C  AT  ION  :  =  A*3.  36? 

C  If  A  AND  i  ARE  VECTORS/  C  IS  THE  DOT  9R0DUCT  A. 3  343 

C  364 

C  ARGUMENTS:  *65 

C  A:  1  A T  9 1 X  OF  SIZE  <L/N).  Z6& 

C  3:  v  A  T  R  I X  0‘  SIZE  (L/.M).  347 

C  C:  PRODUCT  MATRIX  Of  SIZE  <n»M).  36.1 

C  *»'/►./ L :  SIZES  Of  "aTRICES  A,3,C.  369 

‘  370 

C  (NOTE:  S  JO  ROUTINE  ASSUMES  THAT  THE  FIRST  DIMENSION  371 

C  OF  A/6  ANO  C  IN  THE  ZALLINS  PROSRAM  IS  L/L  AND  N. )  372 

C  373 

DIMENSION  A(L/1)/3(L/1)/C(N/1)  374 

DO  10  1*1, N  375 

DO  13  J  s  1  /  M  376 

C(I/J)s0.0  377 

DO  10  K=1/L  373 

10  cc/j)  =  :<i,j>  ♦  a t </ : > » 3 (<., j )  379 

RETURN  330 

END  3S1 
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SUBROUTINE  K)TT(;,3,t/N,l’/L) 

P£HFCI>y,S  MATRIX  "HILT  I  PL  H  A  T 1 0  fi  C  *  A3* 

WHc«£  DIMENSIONS  A » E  AC\,l>  ,  3(*,L>  A.\i  C(,N,M). 

SION  4<N/1  ),3  (X,  1  ),c(N,1  ) 

DO  10  1  =  1, N 
: c  io  j  =  i , m 
CCI/J)=0. 

00  5  K- 1 / L 

5  C(i/J3=  A(I,<)»P(J,K)+CC:#J) 

10  continue 

SiTUR.N 

END 


RtV  01  11/20/72 


3?2 
3?3 
334 
3  3  5 
3  35 
337 
333 
33® 
393 
3  91 

392 

393 

394 

395 
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SU3  ROU  TI NE  DRCYPR  (D/A,  II  /  12/ 13)  396 

C  REV  03  07/03/74  397 

C  SETS  UP  3X3  DIRECTION  COSINE  MATRIX  FOR  GIVEN  YAW/PITCH  AND  ROLL.  393 

C  399 

C  ARGUMENTS:  400 

C  D :  2X3  DIRECTION  COSINE  MATRIX  TO  3£  COK°UTED.  431 

C  A:  ARRAY  0'  LENGTH  3  CONTAINING  S0T  A7ATI0N  ANGLES  (DEGREES).  402 

C  II:  AXIS  OF  POTATION  FOR  10T  ANGLE  (1/3/3  =  X/Y/2)  403 

C  12:  AXIS  OF  ROTATION  FOR  2ND  ANGLE  (1/2/3  =  X,Y,Z>  434 

0  13:  AXIS  OF  ROTATION  FOP  3RD  ANGLE  (1/2/3  =  X,Y,Z)  405 

C  406 

DIMENSION  D(3/3)/A(3J/T(6/3)  407 

RAD IAN=. 01 745  329251  99*33  403 

Y  =  AC1)* RADIAN  409 

?  =  A (2) ‘RADIAN  413 

R  =  A  (3)*RADHN  411 

"  *  6  412 

N  =  3  413 

00  10  1-1/3  414 

DO  5  J=1 / 3  415 

D(I/J)=0.  416 

5  T(I/J)=C.  417 

T(I/I)=1.  413 

10  DC/!)»1  .  419 

IF(Y.EG.0.)30  TO  20  420 

CALL  A0T(T,I1/Y/M)  421 

DC  15  1=1/3  422 

DO  15  J  =  1  / 3  423 

15  0(1/ J )=T(I,J>  424 

20  IF(P.EG.O.O)GO  TO  30  425 

CALL  ROT(T(4/1)/l2/?/P)  425 

CALL  MAT(T(4/1)/T(1/1),D(1/1)/3/3/3/r/M/N)  427 

DO  25  1=1,3  423 

DO  25  J  =  1 / 3  429 

25  T(I,J)=D(i,J>  430 

30  If CR.E5.C.0)  GO  10  40 

CALL  RD:(T(4/1)/I3/R/N)  437 

CALL  MAT  (T  (4/1  )/T(  1/1)  ,0(1/1  )/3/2/3/M/'1/N)  433 

40  CONTINUE  434 

RETURN  435 

END  436 
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SU3R0UTINE  ELIPSSUNDEX/XI/IN)  *57 

c  . . . . . . .  433 

C  43? 

C  THIS  SUBROUTINE  3  £NE  R  ATE  S  1/4  OF  THE  CONTOUR  LINES  FOR  AN  ELLIPSOID  440 
C  441 

C  * . . . •*. ••*•*****•»**. ••»»»•*  442 

DIMENSION  XI  <3/IND£X/I'lDiX)/IN<!NDEX)  443 

C0*1'ion/3LL:pS£ /NSTE»S  (90)  #IELS/A(  3/3/305  /SE£LP(3/P0)#VP(3)/  444 

0(3/3/9j)/DVP(3/3)/k4<3)/NS£3  445 

5ZLTAX*SSST(1/A(1/1»IELp))/NSTEPS(IELP)  445 

D'LTAY=3SPT(1/AC2/2/ISt?) ) / NS  TEAS < I £L°>  4  47 

Sir°1  =  A(1/1/I£LP)  443 

SIVP?  =  ,4(2/2/I£L?)  44? 

a  I M  ?  3  *  A(3/3/IEL?)  430 

s:  101  L=1/IN0'X  451 

X=(.--1)*DSLTAX  452 

N-3  453 

DO  50  <=1/ I  NO  2  X  45* 

Y*(K-1)*0ELTAY  455 

TEX’  =  1-X*X«SI'1P1  456 

TEST  =  TENP  -  Y*Y*SirP2  457 

\ = N  * 1  453 

IKTEST.LT. 0.0)  GO  TO  100  459 

I  =■  SQRTUSST/5INP3)  460 

XI  (1/\/L)  =  X  461 

XI  C#K/L)=Y  462 

XI  (3/N/D-Z  463 

50  CONTINUE  464 

GO  TO  101  465 

100  Y  =  SiU(TSMB/SI'l!02)  466 

XI  (1/N/L)*X  467 

X1(2/\,L)=Y  463 

XI  (3/N/l)*0.0  469 

101  IN(L)=N  470 

RETURN  471 

END  472 
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SUBROUTINE  EXT£ND(P/I/J)  473 

C  DMMON/ INT  ER  S/  N I E (90) / I E (93/ 90)  4 74 

C0MMCN/ELL:PSc/NSTEPS(93)/ISL»/A(3/3/3a)/$E3LP(3/90)/VP(3)/  A  75 

»D(!/3/50)/DVP(3/3)/HA(3)/NStC-  4  76 

DIMENSION  ?(3/2>/?3(3)  477 

COMMCN/PLTT/Sf A C TP / I  NT/ TI ME/ I  COLOR  (  91 > / 0 F S ETX/ 0 FSETY / ZT IME  4  76 

C  4 79 

C  4?0 

C  THIS  SUBROUTINE  FINOS  A  MIDPOINT  FOR  A  LINE  THAT  481 

C  BEGINS  ON  P(1>  AND  ENDS  ON  ?(2>.  482 

C  THIS  NEW  POINT  IS  CHECXED  BY  SUBROUTINE  HYDE .  483 

t  IF  IT  IS  HIDDEN  THEN  484 

C  IF  IT  IS  NOT  HIDDEN  THEN  R(J)*?3  *65 

C  THIS  ALGORITHM  IS  ITERATED  INT  TIMES.  436 

C  UPON  LEAVING  EXTEND-  =>(1)  .IILL  CONTAIN  THE  RESULT.  487 

C  488 

C  NOTE:  P  ARRAY  IS  CHANGED  PY  THIS  SUBROUTINE.  489 

C  490 

C  .........................................  491 

DO  3  IN=1 / INT  492 

DC  1  L=1/3  493 

1  P3(L)*(P(L/2)+P(L/1 >)/2.0  494 

N'JM*NI  c  ( I  EL?)  495 

DC  2  IM  =  1/,NUR!  496 

XK*I£(IM/ltLp)  497 

Ir(XK.LE.JO)  CALL  H YDE  ( XK/P3/ 1  FLA G)  498 

IF  UK. GT. 30)  CALL  HI  D  £  (  KK/P3/ I  FUG)  499 

IFUFL4G.E9.1)  GO  TC  10  500 

2  CONTINUE  501 

10  N=I  502 

IFdFLAG.EO.I)  N  =  J  503 

DO  3  L=1 / 3  504 

3  P(L/N)=B3(L)  505 

RETURN  506 

END  507 
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SU° ROUTINE  Sc  NO CM( CAM:  PA/ FOCUS/ 01 

503 

c  * 

509 

c 

THIS  SU3R0UT INE  SeNERAT E S  A  DIRECTION  COSINE  MATRIX. 

510 

c 

511 

DIMENSION  CAMERA(3)/'0CJS(3)/Z  (31/0(3/3) 

512 

S Jr  =  0.0 

513 

IF  ( FC:u:(1).Ni.O.O.OR.FOC:jS(:>.NE.C.0.3». 

514 

*  CAN£PA(1)..\ 5.0. C.OV. CAMERA!?). NE.O.O)  GO  TO  10 

515 

00  AO  1-1/3 

514 

30  40  J=1/3 

517 

D(I/J)S0.3 

513 

40 

CONTINUE 

519 

0(1/2)*1 . 

523 

0(2/11*1. 

521 

D(3/3)*-1. 

5  22 

GO  TO  999 

523 

;o 

CONTINUE 

524 

C 

525 

DO  100  1=1/3 

526 

Id)  *  FOCUS  ( I )  -  C  A**E  9  A  (  I ) 

527 

130 

SUN  =  SUM  ♦  ZCIW(I) 

523 

SUM  =  SORT (SUM) 

529 

DO  200  1=1/3 

530 

200 

Id)  =  I  (  I  )  /  SUM 

531 

C 

532 

XNORM  =  $4RT(Z(1)*Z(1)  ♦  Z(2)«Z(2)> 

533 

C 

534 

C  FILL  IN  FUST  »OW  OF  0 

535 

C 

536 

0(1/1)  =  1(2)/ X NORM 

537 

0(1/2)  =  -Id)  /  XNORM 

538 

3(1/3)  =  0.0 

539 

C 

540 

C  FILL  IN  SECOND  R0.  OF  0 

541 

C 

542 

0(2/1)  =  1(1  )*Z(3)/XN0RM 

543 

0(2/2)  =  I  (  2  )  *Z  (  3)  /  <*iOR  M 

544 

0(2/3)  =  -XNORM 

545 

c 

546 

C  FILL  IN  THI73  ROW  OF  0 

547 

C 

543 

00  300  :=1/3 

549 

300 

DC/I)  =  Z(I) 

55C 

9  99 

CONTINUE 

551 

RETURN 

552 

e:o 

553 
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SUBROUTINE  HID£<<*/P3/IFLAG)  554 

COMMON/POLYGON/ NPLANE>KFUG/N«PP<90)/POC 3,4,60) /P< J/4/60)  »  555 

•CJNVECC^^O)  ,PO$<2,9C),5I3N<90)  5  5  6 

COMMON /£LLIP$£/,iSTEP$(90),IEL?/A<3/3,30)/$E3LP(3/9Q),VP(3)/  55  7 

*D(3,3,V)),DVP(3,5>,R4(3),NS£S  55? 

0 1 IE  NS  I  ON  *7(2)  /P0(3)  559 

DIMENSION  ®3(3>/°;(3)  56C 

REAL  PPRX.M£(3)/N*RIM£<3)  561 

CALL  TRANS1 <?J,P4)  562 

»7<1 >«?4<1 )/?4<3)  563 

»7<t>*<*4<2>/94  (!)  564 

CALL  TP3INT(P7,K<,IFLAG)  565 

If  (IFLAG.E3.2)  RETURN  566 

C  557 

C  POINT  IS  IN3IDS  POLYGON  CHECK  TO  SE£  566 

C  If  POLYGON  09  POINT  IS  CL3SEP  TO  VIEWPOINT.  569 

C  CALCULATE  TA'J.  570 

:  571 

DO  3  1=1/3  572 

PPAXKE(X)*D(1/I/KK)  573 

If  IP  PRIME (I). £0.0. 3)  PPRIME(I)*.000001  574 

5  CONTINUE  575 

C4LL  ¥AT(DV0/PPRI.ME/NPRIME,2,j,1,3/3/3)  576 

CC  1C  J=1,  3  577 

10  PP(J)*p{J/1/KK-33)-V0(J>  573 

CALL  PAT10V?,PP,PPRIME,3,3,1,3,3,J)  579 

CALL  DOT (N®RIME,PPRIM£/?5/1/1/3)  580 

CALL  D0T(NPRIK£,P4,?6,1 ,1 ,3)  531 

I f L A3  =  2  552 

If (P5/PS.GE.. 99999959)  RE  TURN  533 

I F  L AG*1  554 

RETURN  535 

END  536 
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SUBROUTINE  HYO£(N,R, IfLAG)  587 

C  585 

C  . .  589 

C  590 

C  SUBROUTINE  HYDE  DETERMINES  If  A  POINT  IS  HIDDEN  SY  591 

C  ANOTHER  ELLIPSOID.  592 

C  593 

C  **..*»•*.*********••**•*•••**»*•*••♦**♦•  594 

C  . . . .  595 

C  N*  POSSIBLE  HIDING  ELLIPSOID  NJMiER.  596 

C  8*  VECTOR  TO  PLOTTING  POINT.  597 

C  If LAG*  fLAG  THAT  INDICATES  HIDDEN  LINE  OR  NOT.  598 

C  IfLAG  *  2  =  NOT  HIDDEN  59V 

C  IfLAG  *  1  *  MIDOEN  600 

C  . . *•»•»»»»**»*»***  601 

C0PM0N/tLLIPS£/NSTE?S(9O),IcLP,AA(3,3,3O),S£GLP{3,90),VR(S>,  602 

»D(3,3,93>,DVP(3/3) ,RA(3),NSEG  603 

DIMENSION  Pi  ( j),P2(3),R2(3),S(3>,V(3)  604 

DIMENSION  MU(3),M(3,2),P(3>  605 

DIMENSION  OD(3,3>, VP1 (3)  606 

DIMENSION  »(3)  607 

REAL  M  ,  MU  ,  PAG  608 

C  ASSUME  NOT  HIDDEN.  609 

I  FLAG  *2  610 

C  611 

C  512 

C  613 

C  °UT  Sc  CLP  ( N)  IN  N  *  S  FRA**.  614 

C  PUT  SEGLP(M)  IN  N'S  FRA-E.  615 

C  °UT  R2  IN  N'S  FRAME.  616 

C  PJT  VIEW  POINT  IN  N'S  FRAME.  617 

C  618 

C  619 

CALL  MAT(D(1,1,N>,$cGLP(1/N),?1,S/3,1/3/3/3>  620 

CALL  KATID(1,1,N),SESlP(1,m?>,?:,3,3,1,3,3,3)  621 

I  f ( I  EL3  .LE.  30)  GO  TO  55  622 

00  56  1*1/3  623 

DO  56  J *1  / 3  624 

56  DOU/J)  *  0(1, J/N)  525 

GO  TO  57  626 

55  CONTINUE  627 

CALL  C0TT(D(1,1/N),D(1/1/IcL3)/0D/J,j/3)  628 

57  CONTINUE  629 

CALL  M«T(DD,R, 82,3,3,1,3, 3, 3)  630 

CALL  !'AT(0(1,1,N),VP,V31,3,3,1,3,3,3)  631 

HAS  =  0.0  63? 

C  633 

C  634 

C  635 

C  FIN?  VECTORS  3/ V/4ND  *"J.  636 

0  MU  HILL  3  c  C  0  M  £  A  UNIT  VECTOR  IN  VECTOR  M'S  DIRECTION  637 

C  JR  IN  p'S  OPPOSITE  DIPECTION.  633 

C  639 

C  641 

DO  1  1=1,3  642 

SU)  =  P3(  I)  ♦  R  2  ( I )  -  PI  (I)  643 

V<t)  *  VPKI)  -  »1  (I)  644 

MUCH  =  SC)  -  V(I)  645 

1  MG  =  MAG  ♦  MUC>*»2  646 

MAS  *  SORT (MAG)  •  647 

C  ‘48 
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C  . 64? 

C  650 

C  '1AKE  MJ  4  UNIT  VECTOR.  651 

C  652 

C  . . 653 

C  656 

DO  2  1*1/3  655 

2  'J(l)  *  “JO  /  “AO  456 

4  *  44(1/1/ M)  657 

3*AA(«./2/N)  656 

C  *  44(3/3/;;)  65« 

lt(49S(K'J(1)).GT..0003C0i101)  GO  TO  10  640 

:F(43S(“j(2)).GT..OOOOC2j01  )  GO  TO  20  661 

C*LL  Z(KU/A/6/C/S/M/JfL4C)  662 

30  IFtJFLAG.EQ.I )  RETURN  663 

C  664 

C  . . . . »•*»••*»••.•»*•••»»*••*•*•••»»»•••»•»«»«  665 

C  664 

C  FIND  ®  AND  COWARE  *1  TO  P  TO  DETERMINE  WHAT  POINT  667 

C  IS  CLOSER  TO  THE  VIEW  POINT.  666 

C  669 

C  . •*•* . . . . . . . ********  670 

C  671 

DO  3  1*1/3  672 

3  P(I)«  $<I>  -  V(I)  -  .1(1/1)  673 

CALL  00T(P/“(1/1)/RE3LT1/1/1/3)  674 

CALL  30T(P/(1(1/2)/RESLT2/1/1/3)  675 

I‘(N.ES.!£LP)  GO  TO  400  676 

If (RESLT1  .GT.0.0030G0001 )  I FLA3*1  577 

41  IFIRE3LT2.  ST . 0. 0000 J0001 )  IFL4G*1  678 

RETURN  679 

10  CALL  XYZ(MU/A/tl/C/$/v/JFLAS)  630 

GO  TO  30  631 

20  CALL  YZ(f!'J/A,3/C/S/r/JFLAG)  632 

GO  TO  30  6? 3 

400  IF(A3S(BESLT2>.3T.43$(»ESLT1>)  30  TO  41  434 

RE  SLT2*?£  SLT1  635 

GO  TO  A 1  436 

RETURN  637 

END  63  3 
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SJiADUTISC  INPUTtfTJPE)  64  9 

C)",'ON/PLTT/SFACT9,!NT,TIME,lCClOR<>1>,OFS£TX,0*ScTY,2TI*t  490 

C3.W(/IKTE»S/  NIE(V3),IE<?J,90>  691 

£  jrNCN/'Lll’S£/NiTEaS<90>,IEL*», »(3,  3,  30 1  *  SS  SLP<  5# 90)  ,  V* ( 5)  ,  69  2 

•  D<3,:,*?>,DV»0,!>,Ri(3),.NSEG  6’3 

C'PMCN/»0lYi3N/N?LANf,IFCA5,N*P»C»0),P0(3,4,60>,P<.',4,40>,  4*4 

•f  :,;vcc<:,i,90>,»o»<:,  ?:>,s:c..(  ?o>  4*5 

C0*13\/4T./PLC1’,30>  641 

DIMENSION  j0  (  24/40)  '?7 

DIMENSION  ?J<3)  694 

CDM*GN/D  lUGMOS  JU3<  aO/M  »S  ,  0 E V FL i, ONL INE,  T £ »M, 3  DPS, OFl INE  69  9 

C:'«MDN/VIS*P/VP3  0),CVP  j(3,  3),I  VB,V?2<3>  TOO 

CJMMO.N/COJSCT/  NP, **1(3,5,60)  »C1 

COMMON  /9£MDVc/NP9£9'/IPcvOV(30)  B02 

30U4LE  PRECISION  CTIME/ETI-E  703 

INTEGER  DEVILS  704 

IFC’ltG.MS.I)  GO  TO  400  735 

°cAD(5,70)  NF  4ST,1IBRE*',NISG  706 

aRIT£(6/7C)  Nf  AST/NPJ’EM/NISG  707 

Pc 43 (  5/ 72)  CREMjvd  )  ,1*1,  NP?£-1>  70S 

W<IT£(6,72)  (HcMDV(I),I*1  /NPRi.l)  7J9 

72  F09P,AT(3(4  0I2/))  710 

PEA»(1,SND*300)  N$£S/NP/°L/9D/  (((MPL(I,J/K),!*1/3),  0*1, 5), K  711 

•*1/30)  712 

£  7)5 

39  R£AD(1,END*7J0)  TIM-/  ( ( S £ 3LB ( I , J  ) /  1*1/3)/  J*1/30)/  714 

*  (((DC2/J/K)/  1*1,3)/  J  *  1 / 3 ) /  4*1/53)  715 

IT!Mc»TIPi»1393000.P.5  716 

ZTIMS*ITIM£/1003D03.00  717 

C  71 8 

70  FORMAT(3I2)  71? 

IFCZT1ME.LT.CTIME)  GO  TO  39  720 

SEAD(5,40)  NS?  721 

43  f 3»VAT ( I  2)  722 

IF(NSB.:C.3>  30  TO  46  7 23 

DO  45  L  *  1  /  N  3  B  724 

<  *  *1  °  ♦  L  725 

It*I3*NB*L  724 

3tAD(5/41  )  NP?P(II),rBL(1,1,lO  727 

41  F0«MATU1,I2>  72? 

NS!9£S*N??BCX)  72? 

DO  45  0*1/ N SIDES  730 

?cAD(5/A:)  (P3(I/J/0/I*1/3)  731 

42  fORMATOMO.O)  732 

45  CCNTI.NUt  733 

46  N?LASE*N®*K’SP  734 

S3  10:  J *1/95 £3  735 

DO  mo  1*1  /3  736 

ICO  A(I,I,0>*1.O/9D(I,O>**2  737 

00  203  J*1/\Se3  733 

CALL  D0T(0(1,1,J)/4t(4,j),0C, 3/1/3)  739 

DO  203  1*1/3  740 

200  S£GLP(I/J)*SEgL®(I/J)*DD(X)  741 

IF(IDE3o5(i>.E9.1>  .PIT£<6,6)  N3E3,NBLANE  742 

6  FJPMAT(lHl/*NUK9E9  CF  SEGMENTS  «  '/1 2/'/  NUMBER  Of  BL ANES  ■  ',I2>  743 

NS  EG  *  NSEG-NFAST  744 

READ( 5/ 331 )  ( I C DLOP ( I ) , I *1 , 30)  745 

301  F3«NAT(2(5X,I5>>  746 

1 I  *NPL AN£*  JO  747 

PSAD(5,:01>  (ICDLD?C)/I*31/90)  745 

PEAD(5, 331  )  Z t  OLD  R ( 91 )  749 

:F(I0EdJG(3).£J.1 )  dPI TE ( 6/ 71 )  N5EG  750 
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71  FCR“*T(1X, '  T MS  NJKOcR  Of  SEGMENTS  TO  3s  PLOTTED  «  ’,12)  701 

R. *0(5,1)  (NST£RS(IPP>,IPP*1,N$£G>  752 

Rc*D(5,1)  (NSTePS(IPRMC),IPP«1,NPLANE>  755 

1  fJ*rAT(I3I2)  754 

If(I5En'JG(3).E9.1)  WRIT£(9,2>  ( NS TEP S < I PP> , I®P«1 ,N SE 5)  755 

2  f  Ok  VA  T  (13x,’NUK3£R  Of  DIVISIONS  ALONG  *  P  AD  US  •  ,  /  2X,  3  J I  5)  754 

I'(IDE3JG(S).F0.1>  wPITE  (6,  55)  (NSTE "G ( I P »♦ 3G> , I *»*  1  , HPL ANE )  757 

55  fO*S*T(1SA,' KUK3ER  JF  DIVISIONS  AL0N3  A  SIDE  ’,/2X,30I3>  753 

Ri*D(i,1 1 )  IM,S‘ACTP  759 

1 1  f  JRMATC  3,7X,F1C.  2)  760 

kR’TE(6,11)  INf,SFACTR  741 

®E  AD ( 5, 751 >  Cf«ETX,Of SETY  762 

W»!Ts(5,9Jl>  0‘SCTX,Cf SETT  763 

901  F0R*AT(’F10.0)  764 

If C0E3UG(3).ES.1)  WRITE(6,90?)  Of SETX, OF S F TT  765 

902  f  34MAT(1X,'0fStTX*  ’ , f  1  J . 3, 4X, ’ 0 F S E T Y  =  ’,F10.3)  7e>6 

If (I0EEUS(3).£3.1>  URIT£(6,1 ?)  CFACTR,INT  767 

12  MRKATdX, ’SCALE  FACTOR  *  ’,F10.2,2X,*  ITERATION  NUMBER  *  *,I3)  768 

R£AC(5,13)  VP,RA,IVP,ICODE  769 

13  F3RKAT(6F1J. 0,2113)  770 

C  771 

C  I  CODE  *  0  :  ROLL,  PITCH,  AND  YAW  ANGLES  ARE  SUPPLIED  IN  RA  ARRAY.  772 

C  773 

C  ICODE  =  1  :  DIRECTION  COSINE  MATRIX  SUPPLIED  AS  INPUT. RA  APPAY  IS  1ST  774 

C  ROW  OF  MATRIX. THE  NEXT  CARD  CONTAINS  THE  2ND  AND  3RD  ROWS  775 

C  776 

C  ICODE  *  2  :  POINT  AT  WHICH  VIEWPOINT  2-AXIS  IS  TO  AIK  IS  SUPPLIED  777 

C  IN  RA  ARRAY.  773 

C  779 

IfdCODE  .HE.  0)  GO  TO  500  780 

If (IDE3UG(3).CS.1)  WRIT£(6,4)  VP.RA  7S1 

4  FOR -MAT(1X, ’VIEWPOINT  VECTOR  ( ’ , f 1 0. 1 , • , ’ , F 1 0. 1 , * ,  ’  , f  1  0 . 1  , * ) ’ 71 X  732 

1 ' ROT  A  T I  ON  OF  VIEWPOINT  RELATIVE  TO  SEGMENT  COORDINATE  SYSTEM’/IX  733 

2'TmESE  ROTATIONS  MUST  RE  DETERMINED  3Y  PERFORMING  ROLL  NOTION  FIRS  734 

JT’/IX'THEM  A  PITCH  “OTION  ANO  THEN  THE  YAW  MCTION’/IOX  735 

4’RCLL  *  ’, F13.1, IX, ’DEG. ’,5X, ’PITCH  *  ’  , FI  0 . 1 , 1 X, ’ DE 5. * ,  7P5 

55X, ' YAk  *  Sfn.lwIXw'PsS.')  737 

CALL  P?CYPR(CV»,RA,1,2, 3)  783 

GO  TO  550  78> 

50C  IFUCODE  .ED.  2)  CALL  >,E  HOC «  ( VP,  ?  A,  DV  P)  790 

IFCCOD:  .£0.  2)  GO  TO  5  50  791 

DO  501  J J J  =  1 , 3  992 

501  DVP(1 , J  J  J )  *  R  A ( J  J  J )  7°  3 

PEAD(5,13)  ((DVP(I,J),J=1,3),I=2,3)  794 

W>  ITS  (6,14)  <OVP(I,J>,J  *1,3 >,1*1,3)  795 

14  FORM,  AT(3(7*  ’,3(1X,F10.0) ))  796 

Ir(IDEELG(3).E».1)  l»*IT  E  <0,1  5)  (  V»  (  I )  ,  I  *1 , 3)  797 

15  fORKAK’  VIEW  POm  VECTOR  (',  F 1  0. 1  FI  0 . 1 ,’,  1 ,  F 1  0.  1  ,').•, /,  79  3 

-  '  VIEW  POINT  ORIENTATION  DEFINED  IN  DIRECTION  COSINE  MATRIX  FORK  799 

-.’)  303 

550  CONTINUE  501 

DO  3C  J *1 , 3  302 

V’C  ( J  )  *'/p(  J )  333 

DO  3C  1*1,3  334 

30  DVP3(J,I)*0VP(J,I)  805 

IFUFLA3  .NE.  1)  RETURN  334 

I  F (  I  DE3UG ( 3 ) .£0.1)  w®ITE(6,40)  307 

60  FSRMATdX, •••**••••••**••««•»*•'/, IX,  SC* 

1 ' V I EWPOINT  DIRECTION  COSINE  MATRIX')  PO® 

IF(IDEbJ5(3).ED.1)  RITE(6,53)((DVP(I,J),J=1,3),I*1,3)  813 

If ( I0EeuG(3).EJ.1)  WRITt<6,54)  •  811 

CO  20  IXE«1,NSE3  51? 

:c(IDEElUG(3).£0.1>  WRITE  (6,50)  IX  E,  (  SE  SL®  ( I ,  I XE  > ,  I  *  1 , 3 )  813 
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50  FO^^IAT  C 1 *  *»**•»•*••*******»**  * /I  X • SEGMENT  *  •/ I 3/ 5X/ ' S E CLP  *  <*  814 

-Fl0.3/2<,/*F10.3)/,/'///1X/,A  MATRIX* /T50/* DIRECTION  COSINE  MATRIX  815 

-  */)  816 

DO  400  111*1/3  81? 

I F ( 1 0E0U5 ( 5 ) . E  9. 1 )  WRITE! 4/ 51)  <  A  <  1 1 1  /  J  /  IX  E  >  /  J  »1  /  3  >  /  <  D  <  1 1 1  /  J  /  I X  E)  ,  818 

♦J  =  1  /  3)  519 

400  CONTINUE  820 

51  t  oaMATO(2X/F3.5)/T50/3<2X/ f  5.6) )  821 

53  F  0  'AT13I2X/F9.6))  822 

IF  {  I JE3U3<3).ti.1)  w°ITE(c/54)  323 

54  F3AVAT(1X///1X' ***»*****•**•**•****')  824 

20  CONTINUE  325 

RETURN  326 

600  CONTINUE  8  27 

IFtZTIKE.LT.CTINE)  60  TO  675  328 

IF (ISW1.EQ.3)  50  TO  630  329 

CO  650  J=1/NSES  830 

CALL  D0T(D(1/1/J)/3D(4/J)/0D/3/1/3)  831 

00  650  1=1/3  832 

650  SE jLPCI/ J)*SEGL?(I/J)*DD(I)  833 

I  S  W 1  =  0 

I F  L A3  =  5  835 

RETURN  336 

675  REAOd/S  NO  =  700)  TI  HE  /  (  <  S  EGLPC  I  /  J  )  /  I  =  1  /  3)  /  J  *1  /  30)/  337 

*  UOCI/J/X)/  1  =  1/3)/  J  =  1  /  3)  /  K  =  1  /  30)  833 

ITI*E=T:k:*103C300.+.5  839 

:TINE  =  ITI'!:/1000000.D0  340 

C  341 

13*1=1  842 

30  TO  600  843 

700  WRIT£(6/720)  344 

720  F  OR  TAT ( 1 X/ *  END  OF  DATA  REACHEO.*)  845 

I S  W 1  * 1  346 

STOP  347 

300  WRITE<6/820)  343 

E20  FOANATUX/’SO  DATA  ON  TAPE.*)  ?49 

STOP  350 

630  IFLA5=10  851 

RETURN  852 

END  353 
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SUBROUTINE  L$EGINT(P1/P2/R1/R2/IFLAG) 


35  4 

C  855 

C  THIS  S'JEROUTINE  DETERMINES  IF  TWO  LINE  SEGMENTS/  PIPE  AND  R1R2,  ’56 

C  INTERSECT.  857 

C  ALL  PARALLEL  LINE  SEGMENTS/ WHETHER  COINCIDENT  OR  NCT/  ARE  35? 

C  CONSIDERED  TO  3E  NON-INTERSECTING.  859 

C  CASE  1  IS  CONSIDERED  TC  5E  THE  REGULAR  CONFIGURATION.  860 

C  THE  SPECIAL  CASES  ARE  AS  FOLLOWS:  861 

C  2)  ONE  LINE  IS  VERTICAL  562 

C  3)  20 Th  LINES  ARE  VERTICAL  863 

C  4)  BOTH  LINES  ARE  HORIZONTAL  864 

C  5)  20TH  LINES  HAVE  THE  SAME  NON-ZERO  SLOPE  865 

C  6)  ONE  LINE  IS  VERTICAL/  THE. OTHER  IS  HORIZONTAL  866 

C  567 

C  S63 

C  I  FLAG  =  1  INDICATES  INTERSECTION,*  IFLAG  =  3  INDICATES  NO  INTERSECTION.  369 

C  370 

DIMENSION  PI (2)/P2(2)/R1  (2)/R2(2)/P(4/2)/T(2)  871 

REAL  N ( 2 )  372 

I FL  A3  =  0  373 

C  374 

C  SET  UP  ARRAYS  875 

DO  1  1=1/2  376 

PCI/I)  =  P 1 ( I )  877 

P  (  2  ,  I )  =  ?2(I)  378 

°  (  3  /  I )  =  Rid)  3*’9 

1  P  (  4 /  I )  =  R2C)  880 

C  831 

C  DETERMINE  IF  CASE  3  832 

IF(A3S(P<1,1)-P(2,1)).LT.1  .=-11  .AND. AES < P C 2 / 1 )~p 1 4, 1  )  )  .LT.  583 

♦l.E-11)  RETURN  834 

C  355 

C  DETERMINE  IF  CASE  4  386 

IF(A3S(P(1/2)-P(2/2)).LT.1  .E-11.AN0. A3S<P(3/2)-»{4/2)).LT.  837 

♦l.E-11)  RETURN  883 

C  839 

C  DETERMINE  Ic  CASE  6  890 

DO  2  1=1/2  891 

J  =  3-1  392 

I?(A2S(P(1/I)-?<2/I)).LT. l.E-11. AND. ASS(P(3/J)-r{4/J)).LT.  89  3 

♦l.E-11)  GO  TC  10  894 

2  CONTINUE  895 

C  S96 

C  DETERMINE  IF  CASE  2  £97 

IF(49S(?(1/1)-°<2,1)).LT. l.E-11. OR. ASS(P(3/1)-P<4,1)).LT.  893 

♦l.E-11)  GC  TO  6  899 

GO  TO  5  900 

6  DO  3  1=1/4  9C1 

TEW.P  =  p(!,1)  902 

P(I/1)  =  “(1/2)  903 

3  P ( I / 2)  =  TEMP  904 

C  905 

C  REGULAR  RROCEDU’E  906 

5  DO  4  1=1,2  907 

11  =  2*1  933 

12  =2*1-1  939 

4  MCI)  =  (P(I1,2)  -  R(I2,2))/<P<:1,1)  -  P ( 12/1 > )  910 

C  911 

C  CHECK  FOR  CASE  5  «12 

IF (A3S(M(1 )-M(2) >.LT.  l.E-11)  RETURN  913 

X  =  ( P  ( 3/  2 )  -  P  ( 1  /  2 )  ♦  M(1)*P(1/1)  -  M(2)*p(3,1))/(v.  (1)  -  M(2))  914 

DO  7  1=1/2  ®15 
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7  TO)  =  (X  -  ?(2*I*1/1))/(P<2*I/1 )  -  P(2*I-1/1>)  916 

20  ITtTMJ.GT.O  •  AN5 .  T(2>.3T.3  .AND.  TtD.LT.I  .  AND.  T(2).LT.1>  917 

I  FIAG*1  918 

RETURN  919 

920 

CASE  6  PROCEDURE  921 

10  IF(A3S(9(1/1)”°(2/1)).LT.1.:*11)  GC  TC  11  923 

J=1  923 

1=2  924 

GO  TO  12  925 

11  1=1  926 

J  =2  927 

12  T ( 1 )  =  (P(3,J>  -  "(1 , J ) ) / (P(2, J )  -  PC1,J>)  923 

T ( 2 )  =  (P(1, I)  -  P(3/I))/(P(4,I)  -  p  ( 3, 1 ) )  929 

GO  TO  29  930 

END  931 


SUBROUTINE  IMAT<A/9/C/LL/MM,NN/JA,JS,JC)  992 

C  REV  03  05/31/73  933 

C  PERFORMS  MATRIX  MULTIPLICATION  C  «  A9.  934 

C  935 

C  ARGUMENTS:  V34 

C  A:  MATRIX  Of  SIZE  <L/M>.  937 

£  9:  MATRIX  0*  SIZE  <V.,N>.  935 

C  C:  PRODUCT  MATRIX  Of  SIZE  <L,N>.  93? 

£  L t M/ N j  SIZES  Of  MATRICES  A/3/C.  940 

C  LA,l2*LC:  1ST  DIMENSION  Of  A/j#C  IN  CALLING  PROGRAM.  941 

C  942 

DIMENSION  ACJA/1)/P(je,1)/C(JC,1)  943 

DC  20  L=1/LL  944 

00  10  N * 1 / N N  945 

S  ^  0.0  946 

DO  5  M*1 /MM  947 

5  SsSMA(L/',)*3(M/N)  948 

C<L/N)*S  949 

10  CONTINUE  950 

20  CONTINUE  951 

RETURN  952 

END  95J 
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SUBROUTINE  NFRAHE  954 

C  955 

C  This  ROUTINE  PERFORMS  THE  END  OF  FPANE  HANDLING  FOR  THE  3DRS  95& 

c  957 

INTEGER*!  ENOFRA/NASK/STATUS  958 

DATA  £N3FRA/HAJK/IFFFF,ZFfFF/  959 

CALL  DOLWHCS/I/ENDFRA/HASK/STATUS)  960 

CALL  PLOTSCH/N/LU)  941 

RETURN  962 

END  965 
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SCOUTIN'  OV£RLAP(III/«K/P!‘LAG)  964 

DIMENSION  p1(2>/P2(2)/»1(2>/R2(2)  965 

DI«EN3I0'I  PP 2(2)  966 

C3“NON7POLYGON/NPLAN-/IFLAG/N»PP(9O>/PO(3/4/6O)/P(3/4/60)/  967 

•CCNVE C ( 2/4/90) /R03(2/«0>/ SIGN (93)  968 

C  960 

C  0VE9L4°  TAKES  OBJECTS  I  AND  K  AND  TESTS  FOR  ANY  OVERLAP  ON  THE  970 

C  PROJECTION  PLANE.  071 

C  “FLAG  WILL  3E  RETURNED  TO  INDICATE  IF  OVERLAP  OR  NCT.  972 

C  M  F  L  A  3  *  3  MEANS  NO  OVERLAP  973 

C  MFLAG*1  MEANS  OVERLAP  974 

C  975 

I  *  III  976 

K  *  KICK  977 

5  CONTINUE  973 

DJ  10  J  *1 / 2  979 

13  »°2 ( J )  =  POS(J/K)  930 

C  931 

C  GO  AROUND  THE  RINGS  9<32 

C  933 

NPTS1  *  NPPP(I)  964 

NPTS2  *  NPPP(K)  985 

DO  200  J*1  /NPT32  936 

CALL  TPCINT(PP2/I/HFLA3)  987 

IF (MFLA3. E  3.1 )  RETURN  9?3 

DO  200  N *1/2  9?9 

200  PP2(N)  *  P?2(N)  ♦  CONVEC (N/J/K)  990 

C  991 

C  CHECKED  ALL  POINTS  AND  FOUND  THEY  HERE  ALL  OUTSIDE.  992 

C  993 

C  NEXT/  CHECK  FOR  INTERSECTING  LINE  SEGMENTS .  99« 

C  995 

DO  60  11*1/2  996 

?1CI)  =  POSdI/I)  097 

60  RKII)  *  ?OS(II/K)  993 

DO  61  L  =  1/r»?TS1  999 

DC  62  11=1/2  1000 

62  P2(  1 1 )  *  RKII)  ♦  CONVEC  (II /L/  I )  10C1 

DO  63  J  *  1 / NPTS2  1002 

DO  64  11=1/2  1003 

64  R2(II)  =  RKII)  ♦  CONVEC  (II/J/K)  1004 

CALL  LSEGINT(P1/P2/R1/R2,yFLAG>  1005 

IKMFLAG  .E«.  1)  RETURN  1006 

PK1)  *  R2(1)  1007 

63  A1(2)  =  R2(2)  1003 

P1(1)  =  ?2(1)  1009 

ol  PIC)  =  P»(2)  1310 

I F  ( I  .NE.  Ill)  RETURN  1011 

I  *  KKK  1  Cl  2 

K  *  III  1013 

33  TO  '  1014 

END  1015 
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SUBROUTINE  PLPLN ( S£ 3/ INDEX 2 )  1016 

C  ..................... . .  1317 

c  101? 

C  THIS  SUBROUTINE  PLOTS  THE  PLANES.  131? 

:  iooo 

c  . .  1021 

COM'10K/ELLI?$£/NST£S»S<90)/IEL?/AA(3/3/3  3)/:EGLP<3/90)/ Va<3)/  102  2 

»D(3/3/92)/0VP(3/3)/S*(3)/NS£G  1023 

CO'17ON/P0LY6ON/N?LAN£/IFLAG/N?PP(90)/P0(3/4/6O)/a{3/4/6Q)/  1024 

•C0NVEC(7/4/?Q)/a33(2/9C)/$ISN(90)  1025 

CO«*HON/03UG/IOE3Jj(80)/N:SS/OiVFLS/ONLINs/T£RP./231S/ OPLINE  10  26 

5 1  PENSION  S EG (  3/  333  3)  1  02  7 

INTEGER  02VFLG/0NLINE/TeRK/35RS/0FLINE  1023 

C3fiHON/PLTT/SFACTR/INT/TI*1E/IC0LCia<91)/OF3ETX/O«S  =  TY/2TIPE  1029 

COTPON  /R£M0VE/N?R?FvIR£!»3V<30>  1030 

I F (NPLANE  .£3.  0)  RETURN  1031 

SEGC1/1)  *  0.0  1032 

S£G(2/1)  *  0.0  1033 

SE3<3/1)  *  0.0  1034 

DO  500  LL  =  1  / NPLANE  1035 

DO  100  1*1 /NPREM  1035 

IF  (LL.EQ.IREIOVtl))  GO  TO  500  1037 

100  CONTINUE  1033 

L  *  LL  ♦  30  1039 

IFOEVFLG.E3.OFLINE.OR.DEVFLG.E3.3DRS)  CALL  NEWPEN C I C OLOP CL ) )  1C40 

NUP  =  NaPP(L)*N3TEPS(L)  ♦  1  1041 

A  *  1 . /N  STEPS (L)  1042 

NSIDES  =  NPPP(L)  1043 

DO  400  K*1/NSID£S  1044 

<K  *  <  *  1  1045 

I F ( X  . E 3.  NSIDES)  KX  =  1  1046 

11  *  (K~1) *N$TE°S(L)  ♦  2  1047 

12  =  II  ♦  NSTEPS(L)  -  1  1343 

DO  400  1*11/12  1049 

DO  40C  J*1/3  1350 

400  S£G(J/I)  *  SEG(J/I“1)  ♦  A*(?(j,kY/LL)”3(J/K/LL)>  1051 

I°HN  *  3  1052 

IELP  =  L  1353 

CALL  ?.NTPLT<$E3<1 /I  )  /  I  «>E\/ I  ND  5X2/ NUN)  1054 

500  CONTINUE  1055 

RETURN  1056 

END  1057 
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SU9R0UTINE  PNTPLT < S £0/ I Pc N/ I NO SX2/NPTS > 

1058 

V 

t 

1059 

C 

1063 

c 

POINT  PLOT  SL3R0UTINE. 

1061 

c 

1062 

c 

« 

1363 

COMMON/ PLTT/SF  AC  TR/I.NT/TIME/ICOLOR  (91  > / OF SS TX/ OF  SET Y/ ZT IME 

1064 

10  M.  MON/ELLIPSE/ NSTEPS(93)/IELP/A(3/3/30)/»EGLP(3/ 9 0)/VP(J)/ 

1065 

♦DE3/3/93)/DVOE3/3)/RAEl)/Ni£G 

1066 

COMMON/ INTERS/  N  I E  E  90 )  /  IE  (9 3/ 90 } 

1367 

DIMENSION  PE3)/PPE3/2)/P?PE3) 

1063 

DIMENSION  SEC (3/3333) 

1069 

DATA  Y“IN/0.0//YMAX/11. 3//IPL0T/1/ 

1373 

DATA  IFIRST/O/ 

1C71 

DATA  ITW0/2//ITNREE/3/ 

1072 

IF  (IFI3ST.E9.0)  REA0<5/1000)XPIN/X«AX 

1373 

1300 

FORMAT ( 2F1 0. 2) 

1074 

IF  EIFI.RST.E3.0)  WR I TE ( o, 1 001 ) XMIN/XMAX 

1075 

1001 

FORMATE*  XMIN/XMAX*' /2E1X/F10. 3)) 

1076 

IFI RST*1 

1077 

LFLA3*2 

1073 

IFLA3*2 

1079 

NtU?EN*0 

1080 

DO  100  IPNT*1/NPTS 

1031 

INU«*NI£(IELP) 

1332 

IF (INUM.53.0)  GO  TO  61 

103J 

DO  60  <*1 / I NUM 

1034 

KX*I£ EK/IELP) 

1085 

IFEXK.LE.30)  CALL  HY?EEKK/S£jE1/IPNT)/IFLA3) 

1036 

IFEKK.GT.30)  CALL  HIDE E <K/ SEG £1 / I?NT> / I  FLAG) 

1097 

IFEIFLAG.Ea. 1)  GO  TO  61 

1333 

60 

CONTINUE 

1339 

61 

IFEK.GT. INUM  .OR.  X.fO.1)  GO  TO  62 

139J 

ITEM*  *  IEE1/IELP) 

1091 

IEEI/IcLP)  *  IEEK/IELP) 

1092 

IEEK/ZELP)  *  ITEMP 

1393 

62 

IFEIFLAG  .NE.  LFLAG)  GO  TO  200 

1394 

70 

IFEIFLAG.EQ.1 )  GO  TO  400 

1095 

LFLAG*2 

1096 

CALL  TRANG1 E3EGE1/IPNT)/PPP) 

109  7 

X*-PP°C1)*SFACTR/PP?C3)+0FSETX 

1093 

Y*PPPE2) *SFACTR/»P«E3)*0FSETT 

1090 

IF  EX.G£.XMIN.AND.X.LE.XMAX.AND« 

1130 

1  Y.Gt.YMIN.  ANO.Y.LE.YriAX)  SO  TO  71 

1101 

CALL  CLIPE  X,Y/ XSAV/ TSA Y, X*I N/XMAX/ Y“I  N/ TMAX/ I  PEN/ IPL3T) 

1132 

NEV.P£N*-2 

1133 

IPEN*3 

1104 

GO  TO  75 

1105 

71 

CONTINUE 

1106 

IF  EIPNT.NE.1) 

1107 

1  CALL  PREPLTEX/Y/XSAV/YSAV/XMIN/XPAA/YMIN/YMAX/IPEN/NEUPEN) 

1103 

CALL  PLOT  EX/Y/IPeN) 

1109 

IP^0T*1 

1113 

NEWPEN*2 

1111 

IPEN*2 

111? 

75 

CONTINUE 

1113 

XSA V*X 

1114 

Y3AV*Y 

1115 

103 

CONTINUE 

1116 

RETURN 

1117 

200 

IFEIPNT.E9.1)  GO  TO  70 

1118 

DO  250  IJ*1/3 

1119 

PP<JJ/1)«SEG<IJ/IPNT-1> 

1120 

250 

‘»P<IJ/2)«S£i<IJ/XPNT> 

1121 

CALL  EXTESOIPP/IFLAG/LFLAG) 

1122 

DO  260  ! J*1s3 

1123 

260 

P<IJ)»PP<!J/IFLAG> 

1124 

CAU  TAAHS1(»/P»»> 

1125 

X«-»»o(i)#sFACTA/P*P(J)*CfSETX 

1126 

Y«P*P<2>«SFACTH/PPPO>*0FSETY 

1127 

I  F  (L FLAG. £2.1 )  GO  TO  350 

1125 

IF  IX.CE.XNn.AND.t.LE.XNAX.ANO. 

1129 

y.ge.ymn.ano.t.le  .y/ax)  go  to 

2o1 

1130 

CALL  CLIP(X/Y;XSAV/VGAV*XP[N«X.MAX«V*IN#VNAX/!Trf9#XPL0T) 

1131 

NE«P£M»"3 

1132 

GO  TO  265 

1133 

261 

CONTINUE 

1134 

IF  (IPNT.NE.1) 

1135 

CALL  PAEPLT(X/T/X5AV#TSAV/XNIN/*NtX»TNIN/TNAX/I*EN»NE'WPEN) 

1136 

CALL  PLOKX/T/2) 

1137 

I  *»LOT -1 

1132 

NcWPEN* 3 

1139 

265 

CONTINUE 

1U3 

IPEN«3 

1U1 

X  SAV*  X 

1142 

YSAV«Y 

1143 

LFLA3«1 

1144 

GO  TO  100 

“114  5 

350 

CONTINUE 

1146 

IF  <X.G£.X*!IN.AND.X.L£.XrAX.  AND. 

1147 

1 

Y.GE.YPIN.AND.Y.LE.YNAX)  50  TO 

j  5 1 

114.5 

CALL  CLI°<X,Y,XSAV#Y3AV/XrlN/X«l«X*YPIN»YHAX,ITM?EE» IPLOT) 

in? 

NS«PEN*>-2 

115C 

I?EN*3 

1151 

GO  TO  355 

1152 

351 

CONTINUE 

1153 

IF  (JONT.NE.I) 

1154 

1 

CAlL  ?®EPLT(X/Y,XSAV^Y  JAV»X*I*l,X«'AX,Yrif:,Y’AAX#I»EN»NE,.PEN) 

1155 

CALL  PL0T(X/Y/3) 

1156 

IPL0T«1 

115’ 

NE  <IP;N«2 

115  6 

:pen«2 

115? 

2  55 

CONTINUE 

1160 

XS  AV*X 

1161 

YS  AV-Y 

116? 

GO  TO  20 

116  5 

*00 

I ®EN*  3 

1164 

L  FLAG  *  1 

1165 

CO  TO  100 

116  6 

END 

1147 

SUMOUMNE  POL  TO  11j8 

C  1169 

(  *OLTO  &ENE*AT£J  DIAEC  T I  ON  COSINE  N6Tf»ICI*S  9  09  TnE  POLYGONS.  1170 

C  THE  •  MIS  09  TNI  POLYGON  COO*DINATI  SYSTEM  IS  THE  1171 

C  NO«N*L  '/E  C  TO*  TO  THE  POLYGON  SUMAC*.  1177 

:  Y  VECTOA  IS  ALIGNED  PITH  ONE  C r  THE  POLYGON  SIDES.  117J 

C  1176 

CO?'’.SN/£LLI»SE/NSTEPS(9')),IEL*,A(3,3,!3),S£GLP(5,73),VP<3>,  1175 

0(3/3/93), D'/p<S/3)/AA<3)/NSEG  1176 

CONr.ON/POLYJON/*>lANE,  IfL*G,NPM(90),P3(  J*6,60>,P(),6,60),  1177 

•  CONViC  (7/6/90) /*0S(2/9L)/S!G  1(93)  1174 

DIMENSION  IND£X(6),D1(  ■»1C)  1179 

EOJI V6LENC £  (0/D1)  1103 

OAT*  INDEX/ 3/6/ 7/1 ,2/5/  1131 

OS  130  L«1/NPLANE  1162 

J*9*L*262  111) 

DO  20  1*1/3  1196 

Dl(J*f62)*P(I/2/L)“P(2/1/L)  ms 

20  Dl(J*I»5)*P(I/3/L)*P(I/1/L>  1136 

CALL  C*0SS(91(JM>,D1(J*6),01(J>)  1187 

SUMD1 *0.0  11P9 

SU*D2«0.0  1197 

DO  30  1*1/3  1190 

SUND1*SUND1*D1(J*I-1)**2  1171 

30  SUM32«SUHC201  (J*l«2)**2  1192 

SUND1 •  S3XT (SUMDl  )  1193 

SUH92*S3»T(SJMD2)  1196 

DO  63  1*1/3  1193 

D1(J6I-1)*D1(J*I-1)/SIMD1  1196 

60  D1(J*I*2)»01(J*I*2)/SUH,32  1197 

CALL  C93SSCD1 (J)/01 (J*J)*01 ( J  ♦  (  ) )  1193 

DO  50  1*1/3  1199 

TE3**D1  ( INDEX ( I)*J)  1200 

01  (INDEX OO) *01 ( INDEX (1*3)0)  1201 

50  D1(INDEX(I*J)6J  J»T!J««  1202 

103  CONTINUE  12?J 

DO  603  J*1/N»LANt  1206 

NU«  «  J  ♦  30  1205 

DO  500  K*1,3  1206 

600  SEGLP(X,NU«)  ■  P(</1/J)  1237 

XETUPN  1203 

END  1209 


SUBROUTINE  P*EPlT(X,Y»XSAV,  Y  S  A  V/ X  »  M/ X  M  AX  ,  YP  IN,  Y«*X, IPEN,NEVP£N>  1210 

. . . . . . * . . . .  1211 

C  1212 

c  this  sujacutins  peapoaps  necesjaay  pen  moves  pop  1st  plot  i2U 

:  APTS*  CLIPPING  121* 

:  i  ?i  s 

. . . . .  1216 

LOGICAL  X3PP/Y0PP  1217 

IP  (N£»»ES.N£.-2>  A£TU«N  1214 

:  1210 

C  OUTSIDE  X  AND/OX  T  -OUNDAP  IS  S?  ?  ?  1220 

:  1221 

X  3  P  P 1 • P  AL4 £.  1222 

T3PP».PALS£.  1223 

IP  (X.LT.<«IN.O».ASAV.LT.XPIN.Oa.  122* 

1  A.ST.(MAX.O».XSAV.ST.Xr.AX)  XCP*«.T»U*.  122S 

IP  (Y.LT.Y.riN.OP.YSAV.LT.YMN.OA.  1226 

1  Y.GT.YNAX.O*.YSAV.CT.Y*AX)  YOFP».T«Ut.  1227 

C  1226 

C  DS1£1«INS  IP  3PP  T0P,33TTCM,*IGHT,3»  L£»T  SIDE  OP  PLCTTE X  1229 

C  1230 

IP  (X.LT.XHIN.O*.XSAV.LT.X«IM)  XL  IPII  T«  X.MIN  1231 

IP  (X.GT.XPAX.OA.XSAV.GT.XNAX)  Xl1P!T*X««X  1232 

IP  (Y.LT.YMIN.01.YSAV.LT.YPIN)  YL:»IT*YMIN  1233 

IP  (Y.GT.YHAX.GH.YSAV.GT.YKAX)  YLIPIT»YNAA  123* 

C  1235 

C  GST  X  AND  Y  POINTS  DEFENDING  ON  WHAT  30U.N9AAISS  OVER  1236 

C  1237 

IP  (XCPP)  YTt*IP«YINTCP{X,Y,XSAV,YSAV,XLi*IT>  1236 

IP  (.NOT.XOPP)  YTENP-YLINIT  1239 

I»  CYCPP)  XTEHP-XINTCPCX, Y,ASAV,YSAV,YL:*IT)  12*0 

IP  (.NOT.YOPP)  XTENP-XLIM*  12*1 

C  12*2 

C  WOVE  PEN  UP  TO  THAT  S*OT  12*3 

C  1?** 

CALL  PLCT(XTEMP,YTt«P,3)  12*3 

IP€N«2  12*6 

«STU»N  12*7 

SND  12*5 


SUJ'OUTINf  MJCL1  12*9 

e  . •••••• . • . ••••***•  1210 

C  1211 

e  this  stmourm  mojccts  fLL!*soi3S  onto  tmi  paojcctiof*  plaic.  1 25 2 

C  1251 

C  •  •••• . . . * . . .  121* 

i,9o,vm),  1255 

•  o(  i,i,  n>,  9vk  :,}>,»*(}),  vie  (  1210 

C9’f1OV/POLr39»/NPLAMC/IFLA6/NPPP(9G)/P0(!/*/60)»P(3/*/S0>/  12  5  2 

•  coNV!cc,<.>n),»oi(2/9'3)»5iof.(?:)  175* 

OIRE'lSIOV  99(3/3>/D9t>(  3/3>/SS(3)/1M>  12*9 

OJIIVJIOA  X(3/3>/*2(2/3)  1 2«3 

Xt*L  LA40A1/LArt942/Hl/*2  12(1 

IKNSfi  .*}.  3)  P  f  T  U  A  ••  1 2o2 

93  100  I  *  1 / N  S  t  J  1261 

CALL  03TT(9(1/1/I)/9VP/00/3/2/  !)  126* 

CALL  '<AT(A(1/1/I>/09/999/3/3/i/3/3/3>  1261 

CALL  0Ct(0(1/1/t)/9OC/93/3/3/3)  1266 

CALL  fm(DVP/DO/99C/3/3/3/3/3/3>  1267 

00  13  <*1/3  1261 

10  SKA)  >  S£  GLP (K/  I  >  -  VACO  126* 

CALL  4AT(0VP/SS/S/3/3/1/3/3/!)  1270 

00  5C  !!*1/3  1271 

IF  (SCI). C9. 0.3)  SIID'U  1272 

SO  C9VT I  M’Jt  1273 

C  1274 

CALL  SOLVA (ODD (1/1) *009(2/1  > / ODD ( 3/ 1 > *000 ( 1 / 3)  / 999  ( 2/ 3 > /D0D( 3/ 3 > /  1271 

•DD0(1/1)/DD9(1/3>/S/P(1/1)/<i(3/1>)  1276 

CALL  SOLVA (900 ( 1  / 2 ) / 090 < 2/ 2 ) / 909 C 3/2 > / 903 C 1 / 3 ) *000(2/ 3)/000C 5/3)/  1277 

*003(2/2)/000(2/3)/S/*(2/2)/A(3/2))  1271 

CALL  SOLVA (030(1/1) *900 (1/2)/ 993(2/1 > *909 ( 2/ 2 ) / 039 ( 3 / 1  ) *t 9 0 ( 3 / 2)/  1279 

*909(1/3)/)99(2/3)/O9D(1/3)/99S(1/1)*2.3*)D9(1/2)*DOO(2/2)/  1243 

•90<X1/3)»3  09(2/3>/S/M1/1)/A(3/3)1  1731 

«(2/1)*3.J  1292 

»(1/2)»J.O  1273 

4(2/3)*4(1,3)  12?* 

93  13  1221 

00  11  IJ«1 ,2  1236 

15  A2(IJ/I<)MS(IJ)«A(IJ/I<)>/(K3)»<(3/IK>)-S(:J)/1(3)  1217 

CALL  S3LVA(*7/A11/«22/A12)  1241 

Ti*P*(Al1«A2;>**2-*.0*(*11*A22-A12**2>  121* 

IF(TE-«.lt.:.C)  T««P*0.0  1293 

TE*?»S:*T(T£«»>  1291 

LA19Al»(Al1*A22*Tt1P)/2.0  1292 

iA’19a;«(A11*A22-tmp)/7.3  1*7* 

A  X 1  *  A 1 2  129* 

ATI •LA40A1-A11  1291 

A42*LADA2-A22  1296 

»V2*A12  1297 

LAK9A1«AaS(LA49*1>  1294 

LAAe*2«A9J(LA)*0A2)  1299 

11  ■  S3  IT  (1.  3/(LA*'9A1*(AXl**2**Y1**2>>)  1300 

12*SJ4T(1.3/(LA,*9*7*(AX2**2*AY2**2)>>  1301 

axl*A*1*Pl  1302 

4X2*ax2**2  1303 

ATI *AY1*H1  130* 

AY2*AY2*r2  1101 

C0fcV£C(1/1/I>»-2. 0*171  1306 

C0NV£C(2/1/I>*-2.0*4V1  1337 

C9\VEC(1/2/I)*-2. 0*9x2  1308 

C0*iV£C(:/2/l)«-2.0*Ar2  13C9 

00  20  :j<1/2  1310 
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C3‘ivtc<:j/ j,i>  •  -c3sv(c(;j/i«n 
(0KVfC(IJ#4»!)  •  -CONVeCUJ*?.  !> 

29  P9uu,n  •  c9Nvtc(ij*s*:>/2.o  *  C9Mvic(u#4«n/2.o»sm>/i(s> 

HPPHU*i 

s::NU>»co><vic<i,i,:>»coimc<r,*,i>- 
ico%vrc(2#i*n*cvivrc(i»2#i) 
no  C3vt;*iu! 

«:ru4\ 

£  NO 


i 

J 
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J^lr*****-'** 


I'JJAOuT  INI  MJ*Lf  1  3?0 

i  di 

. . . id? 

t  wi 

T .  I  J  lU.AOuDS'  .lit  il^u*  T  1  *  COlVtC  AMAT.  D?4 

;•  .ii)  iff.  o"  t*.{  a«3  \i  ii  *  aaati.  idi 

i  D6 

A.  A  AT  A  Ji  TM(  0*1  it  UAL  »3DDC1  *tCT3*J  13*  Nt  •3t»53Si  11*7 

!(  T  •*  A  D‘*T|AL  •A*t*(t:i  U1’l«.  D?4 

A  **AT  ( 3Sv!  (  .ILL  C01TAJ1  T"l  C  3*.TOO*  VCCT3D  1}?4 

>0*  **Cj::T'0  *0L  * i 311 .  1)13 

A.AAT  *14  '.Jtt  CIM  All  *3i:i:0'«  DCTO*S  >»C*  »Wl  1111 

**CJtCT|J'i  *L  All  3«!i:i  f)  *CL  T  ION  *0J1»  *  1.  1 !  5? 

A.. AT  m*  .'LL  CC-.TA11  T*»l  l|i%  T*A»  R.i.AD  M3“  ID) 

l-t  C  *3  D  •  4  33  uC  f  <*?-*)>  <<*»)-*?>.  IDA 

ID) 

. * . 1)3) 

C3«»01/*3LT  i3A/'«*t.A>(  /  I  »  L  AC.  )•>**(  >3>,»3D,4,AG>.»(1,4,63>,  1D7 

•comet  »,%,»:>, *c:t?» »o». m\( »o>  idi 

C0-»3A7lu:*lf/15Tf*$C7j».:CL*#AC  ;,),  )3>,ltiL»(3,  D>,VM3>,  id* 

•3()»l/*C)»3V*t)» ))»**())» 1116  DAO 

Di^mioi  <***?<)>.*»?( ji.mii!)  dai 

JMl*LA1l.(0.:)  •  *  T‘J*1  da? 

00  A 3  1*1, 6*1*11  DA) 

::»:»io  daa 

AM  ).%•>««(  11)  1)A) 

00  ))  * • I » 1*  T 1  1)4) 

03  10  J*1,l  1)47 

***?(J)»*(D*/!)-W*CJ»  DA) 

10  COiDiul  1)4  > 

CALL  XAM3V*,***?,**?»),),1, ),),))  1)10 

I'U.li.D  03  T?  16  Dll 

03  D  J  •  1  #  ?  ID? 

*>cit  j»;  :>•**?<  j  >/**?<  d  i  )i) 

1)  C3NMVJA  IDA 

c.z  to  :i  D)i 

16  33  ?3  •'•!#?  1  Do 

?  3  C3  m  C(4»*-I  I*>!UI/»»!I  )»-»»!  (4)/»*1  D)  1117 

?)  DC  ! j  4*1,1  111* 

);  *»i(j).»*:(j>  Dio 

)>  s:ii(D>»c3S9:C(i»i,m.i3sviCt?«?,j:)-  ido 

icomc  c;,i.: j)»i3sv‘C(i»?.i n  1341 

)>  C3.Tl.Li  136? 

03Aj4*1»?  1)6) 

cciwi  »•*  ;;(4,  n  )-•»?(  j  ;/•*?<))  i.'oa 

ao  cj’.dvj:  i)6i 

••.tj»1  1566 

£•.3  1)67 
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S'JjAOJMNE  P*C(Xl/IS#5E  j/lN5EX,  INDEX?, IHALF)  136  3 

oih-nsisn  xi(3,:n3cx,inuex),:s(:xo£x)/ScG(3,jj3J)  1349 

1  370 

Thu  SU3X0JTINC  *»L 3T  j  A  SSMI  ILLIQUID.  1371 

T » c  M*Lf  3  f  T  H •  *ILI®J3!P  »«.3TT£5  9:®£N3S  U»0N  IHAlf.  H72 

u  HA.f  ■  1  1  0  pld’tes.  1373 

:*  iHAtf  »  ’  x  .lt.  0  ;;  ®worrE3.  1374 

33  1?u  I  ■  Hi'.f  *  INDEX  1375 

L  Ut  ■  l  % 3  i  X- I  ♦  1  1  V6 

lfC<4AL(.C:.C>  L  I N  t  *  I  1377 

NPTi*I'.(LlM  )  1  379 

'jZ  5„  t>1,N®TS  1J79 

33  50  3*1/3  1330 

>UU/0*X1U,K/im)  1331 

I  f  CNAlf  .*  J.7.ANJ.J.E 3.1  >  S-IIJ/Kl-SUO/K)  133? 

50C3'<TI,ijf  1333 

V*\PTS  1 394 

If (LINE. ES. INDEX)  G3  T3  71  1335 

SPT»f.PT*”l  13C6 

03  63  X  « 1 , N®  T  13d7 

*1  1398 

V'NM  1333 

3E0(1/N)«JC3(1,<<)  1390 

i E 5 ( 7/ N ) ■ » £  5 ( It xx )  1331 

60  SE3(3,N)*~S£G( <x  )  3397 

N  •  T  »  N  ”  1  ,393 

03  73  X  *  1 , N?  T  1394 

x».»\PT-X*1  1395 

'*•'**1  1336 

Sc5(1/S)*SEG(1,<X)  1397 

Si*.(7,N>®-SiGC,<K)  1393 

73  S<G(3,N)>SEG< I, x<)  1J99 

71  :®f.*3  1433 

CALI  P‘.  TPLT(S£G,I®£N, INDEX?,  V)  1431 

100  C  JNT  INuE  1432 

XcTUXS  U33 

END  1434 


SU6R0UTINE  ROTIA/L/TH/M) 

c  COMPUTES  POTATION  *AT»IX  A  FOR  ANGLE  TH 

C  ABOUT  X/Y  OR  Z  AXIS  AS  L  =  1/2/  OR  3. 

c 

c  ARGUMENTS: 

C  A:  3X3  ROTATION  vATRIx  TC  r*E  COMPUTED. 

C  L:  1/3  OR  3  TO  ROTATE  AiOUT  Y/Y  OR  I  AXIS. 

:  TH:  ANGLE  OF  ROTATION  III  RADIANS. 

C  .X:  1ST  DIMENSION  OF  A  IN  CALLING  PROS&A". 

c 

DIMENSION  A ( v / 3 ) 

C  =  CCS(TH) 

S  =  SIN(TH) 

IF  (L.EC.2)  S  =  -S 
DO  30  1=1/3 
IF(I. £5.3)30  TO  20 
DO  1C  J=I/2 

a<:,j*i)=o.o 
A<j+i,n=o.o 
IF(I+J*L.NE.5)G0  to  10 
A(I/J+1 )=3 
A ( J  * 1 /I)=-S 
10  CONTINUE 
20  A ( 1/ I ) =  C 

IF(I.E5.L)A(I,I)=1.0 
30  CONTINUE 
RETURN 
END 


1405 

REV  01  08/10/72  1 406 

1407 

1408 

1409 

1410 

1411 

1412 

nn 

1414 

1415 

1416 

1417 
141  3 
1419 
14  20 

1421 

1422 

1423 

1424 

1425 

1426 

1427 
1423 

1429 

1430 

1431 

1432 

1433 
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SUBROUTINE  $Ci.VA(R/AA11/AA22/AA12>  1434 

DIMENSION  R<2*3)  1435 

A11*»(1,1>**2  1436 

A12  =  2.0*?C,1  )*R(1/1)  1437 

A1 3*R(2#  1 ) **2  1 4  it 

A21*R  <1»2)**2  1439 

A?2s2.'J«St2/2)  •  3  <1  /  2)  1440 

A23=AC2,2>*»2  1441 

A31=3C1.,3)**2  1442 

A32*2.0*R(2/3> *3(1/3)  1443 

A33*P (2/3) **2  1444 

D£L=A11  *(A22*A3  5-A23»A32)-A12*(A21*A33-A23*A31  >♦  1445 

1A13*(A21*A32-A22*A31>  1446 

AA11=((A22-A12)*(A33-A23)-(A23-A13)*(A32-A22))/DEL  1447 

AA12=((A23-A13)*CA31-A21)-(A21-A11)*{A33-A23)>/8EL  1448 

AA22*((A21-A11  )*<A32-A22)-(A22-A12>*(A31-A21) >/DEl  14  49 

RETURN  1450 

END  1451 
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SUBROUTINE  SOLVR(A1yA2/A3/A4/A5*A6/A7/A3/P*RX/R2)  1452 

1453 

*****************************************  1454 

1455 

THIS  SUPRCUTINE  HILL  SOLVE  A  SET  OF  SIMULTANEOUS  EQUATIONS  1456 

TO  FIND  COMPONETS  Of  VECTOR  .R  THAT  SATISFY  THE  PROPERTIES  NEEDED  1457 

TO  DET'JRMINE  THE  E3UATICN  OF  THE  PROJECTED  ELLIPSE.  1458 

1459 

SEE  WRITEUP.  1 460 

1461 

. . 1462 

DIMENSION  P(3)  1455 

B  =  A1  *P(1)PA2*PC)*A3*P<3)  1464 

D=A4*P(1)+A5»P<2)»A6*P(3)  1465 

T1=A7»(D/d)**2*A6-2.0»A?*D/b  1456 

T2*2.0*A7O/(3)**2“2.0*A?/3  1467 

T3=A7«(1/fe)**2-1  1453 

RI=(-T2*S'JST(T2**2-4.  j* T1  «T3  ) )  / C 2 . 0*T1 )  1459 

RX*-D*PI/B-1 .0/d  1470 

RETURN  1471 

END  1472 


$U9  “CUT  I N£  TITLE  1473 

DIMENSION  ID<19/23)/ICCL0R(21 >  1474 

CONMON/OiiU3/ID£3LG<30>/NISG/OSVFL5/3NLINE/TERM/3D;»S/OFLINE  1475 

INTEGER  CNLINE/DEVFLG/TERM/ilDRS/OFLINE  1476 

NL1NE*2E  1477 

SII£».335  1478 

X*1  • 175  1479 

Y  * 1 0 . 3  1430 

1 4  i  1 

INITIALIZE  PLOTTING  PACKAGE  1432 

1433 

CALL  PLOT ( 0. 0/3. 0 / *3)  14c4 

RE  AD  <  5/ 1 )  NFRME  1435 

1  FORMAT (12)  1435 

IF (NFRMe. £9.0)  RETURN  1437 

DO  300  <=1/NFR3E  1483 

DO  50  I=1/NLINE  1  4£  9 

READC5/233)  (ID(J/I)/J-1/d)/IC0L0R(I)  1490 

WRIT3(6/2J0)  (ID(J/I)/J*1/5)/ICCL0i»<l>  14»1 

230  F0RMATC7A4/A2/I2)  149  2 

50  CONTINUE  1493 

DO  IOC  I =1 /NLlNE  1494 

Y=Y-.5  1495 

IFOEVFLG.ED.OFLINE.OR.DEVFLG.E9.ODRS)  CALL  NEWP EN  ( I  COLOP ( I  )  )  1496 

CALL  SYK30L(X/Y/SIZc/ID(1/I)/G./33)  1497 

100  CONTINUE  1498 

IFOEVFLG.EO. ONLINE)  CALL  PLOT ( 1 2. / 3. / - 3 )  1499 

IFOEVFL3.cS. TER  <)  CALL  PLO T ( 0. 3/0. 0 / -3 )  1500 

IFOEVFLG.E3.30RS)  CALL  NFRAMc  1501 

IFOEVFLG.ES. OFLINE)  CALL  “LOT  (14./0./-3)  1502 

330  CONTINUE  1503 

RETURN  1504 

END  1595 
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! 


fu 


I 


SUO ’OU  T I NE  T»OIST(PP2/:#:N)  1506 

3 1  MENS  I  OS  f*P2  <  3>/A<J>w?»1  (3)  1507 

CC.V3N/POLYGON/NPLAh£/IFLAG/N»»P<9?),PO(J,A,eO>/Pn»*»6J>/  1508 

*C0iVeC<i/4/9O>/»0$(2/VO)/SI5N(9O)  150? 

1510 

THIS  S.MIG'JTJV:  TESTS  A  POINT  OS  THE  REJECTION  PLANE  1511 

DEFINED  ii  PP2  AGAINST  A  PCLTSON  OS  THE  PROJECTION  PLANE  1512 

DEFINED  IV  I.  A  FIAO  'IN*  !3  RETURNED  TO  INDICATE  IF  THE  1515 

POINT  WAS  INSIDE  OP  OGTSIDE  THE  »CLY33S.  1514 

1515 

IN  *  1  POINT  WAS  INSIDE  POLYGON  1516 

1517 

IN  ■  2  ®3!NT  WAS  OUTSIDE  ROLY30N  1518 

151? 

15*1  1520 

NPTS1*N???(I)  1521 

00  20  J J=1/2  1 522 

20  P?1 <JJ)=®:5(J J/I)  1523 

DO  100  L*1/NPTS1  1524 

00  30  'n  =1/2  1525 

30  R(,\)*PP2(N)-PP1  (N)  1528 

SIGN2*C0NVEC(1 ,L*t)*R(2)-CONVEC(2rL/I>*R<1)  1527 

I F ( S 1 3N2*  SI 3N ( X )  .LT.  C.)  GO  TO  150  1528 

IF(ASS(iIGN2»SIGS(I)).LT.1.E-11)  GO  TO  150  152? 

30  100  %  * 1 / 2  1S30 

100  ?P1  (N)xoai (N)+CONV£C(N/L/I)  1531 

RETURN  1532 

150  IN*2  1533 

RETURN  1534 

END  1535 


SU3  *OUT I  NS  TRAN  SI <  R  /  P  )  1534 

COr"CN/Et.LIPS!/NST!*S<03)/t£LP/A<3/J/33>/SCGLP< J/«0>/V*<3>/  15  V 

0(3/3/9G)/9VP<3/ 3>/RA<3)/N5EG  ij  jj 

CJr  'CN/VISWP/VPDt  3>/O  V«*3<3/3)/lVP/V*2n)  1530 

DIMENSION  DO(3/3}/P<3>/*<3)/A2<;>/SEGL*2<3>  1540 

I  f  <  I  £L®  .GT.  30)  id  IC  10  1S41 

CALL  03TT<D'/P/0<1/1/IcL»>/D.>/i/J/3>  1542 

20  CONTINUE  1543 

CALL  NATID0/R/P2/3/3/1/3/:/3)  I544 

CALL  MATCDVP/SE SLR < 1 / I  EL® > / SEGL? 2/ 3/ ! / 1  /  3/ 3/ 3 >  1545 

DO  1  1*1/3  1544 

1  PtI)*S£uLP2(I) ♦R2(I)”VPc(I)  1547 

RETURN  1543 

10  00  11  1*1/3  1549 

DO  11  J*1/3  1550 

11  DPC/J)  *  DVPCI/J)  1551 

GO  TO  20  1552 

£N9  1553 
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HEAL  FUNCTION  X I NTC ?< X/V/XS AV/ V  SAV/VTEVH )  1554 

. . . . * . . .  1555 

1554 

THIS  FUNCTION  CALCULATES  THE  X  INTERCEPT  AT  1TC*F  1557 

15  55 

. . * . * . . . .  155* 

X1M-XS4V  15^0 

rl«v-vs4v  1541 

IF  (VI  .:<E.  ?.0>  9F  4C  TS* XI / V 1  1542 

IF  (VI. '3. 5.0)  ?FACT«*O.J  1 5  A  5 

ViiVTe^-YSAv  1544 

XtNTCf,*V2».*F4CT,*X$AV  1545 

HETJRN  1564 

•NO  1547 


subroutine  xyiinj. »/0/C/S,k,jflag»  isti 

DIMENSION  *J<3>»  15o9 

till  H.OJ  1 5  T  3 

1571 

. . . . . . . . . **....  15  7’ 

1  '■  7  J 

AL  0  H  A  AND  i  £  T  A  0  i  L  4  T  E  TM-  1  AND  :  :0‘®CNESTS  Of  r  10  1574 

TH?  X  COMPOSE.?  Of  N.  1575 

.  . . . . . . . . .  1  5  ’ 5 

1  5  ’  7 

ALPHA  *  K'JC2>  /  <J(1)  157.'’ 

3ETA  «  *■  J  {  3  >  /  rj(l)  15’9 

1550 

*»* . 15«1 

15-2 

S3LV  t  T  H  £  ESJATIOTi  FOP  ELLIPSE  * N '  TO  FIND  A  ®OINT  OS  15*3 

THE  ;LLI®S£  Th  4  T  .ILL  PET-1NINE  A  VECTOR  N.  15BA 

1  55  5 

. . 1505 

1  5?  7 

T1»A*P*ALPHA«ALPHA*C*SETA«0ETA  1503 

T7«2*A.S(1 )*2*P»ALPha.S(2)®2*C*J£TA«5(S)  1559 

TI«A.J(1 ).«2®3*j(2).*2*C»S(3)»*2-1  1590 

TEH®.  T2  •  T2  -  4  •  T1  •  T3  1591 

159’ 

. . 159J 

1594 

NC  SOLOMON  HEANS  ELLIPCOID  'N'  NOT  TOUCHED  3T  LIKE  OF  1595 

SIGHT  PAY.  1576 

1597 

. . 1593 

1599 

IFITE'-’.LT.C.O)  GO  TO  2  160C 

T; v®  =  SORT ( T  E v?)  1601 

1602 

............................................... ..................  1603 

1604 

F I N  3  TWO  POS5I3LE  M  VECTORS  OcCAJSe  A  SAY  ENTERING  A  1605 

SOLID  VUS  T  ALSO  LEAVE.  160o 

1607 

1*06 

1509 

h(1.1)«<T2*TSF'»)/(2*T1)  1610 

:<{;»i  )®alpha..h(1/1  >  i6ii 

H<3/1  )*  JcTA.*»(1.1  )  lol  ’ 

7:  (1»2)  =  (T?*TEH ’)/(’» T1)  1613 

H(2,:)=4L°HA*r(1,2)  1614 

r(  j/2)«esTA.x(l/2)  1615 

JFHS*0  1616 

.1 1 T  U  P  N  1617 

2  J CLAG*1  1 61 C 

RETURN  1619 

END  1420 
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A  £  A  L  FUNCTION  Y  I  N  T  C  *  (  * ,  Y»  X  5  A  V/ T  i  *  V/ A  Tl*!» )  1621 

* . . . . . . . . *  1622 

1623 

THIS  FdNCTIOS  CALCULATES  Th£  Y  JSlTEACEAT  AT  XTE"P  162* 

1625 

. .  162* 

Al»x-xSAV  1627 

Y1»Y-YSAV  1625 

if  (ai.ne.o.o)  ®fa:t»*yi//i 

If  CXI. ‘3. O.O)  Af  AC  T  X*0.0  16 JJ 

x:«XT£"»-r JAV  1  o  3  1 

y:stcp«x2»pfacta*ysav  isj’ 

« ZTUXN  16JJ 

EN3  16H 


168 


SU?»OUTiNC  TI<XU,A,3,C,  S/N,  JFLAi) 

JHCNSJO'.'  1U<3>,S<3),K<3,2> 

V*AL  PU/H 

C  •  • . . . . . . 

c 

:  AlPHA  »EL*T£J  THE  X  COHPONET  TO  Tm*  »  CO*PON£T 

C 


ALPHAMUtil/XJO 

C  . . . . 

:  solve  the  sjjamon  *oa  ellipse  vain  *»o  to 

£  UNO  A  *OlNT  0.«  T H t  ELLIPSE  THAT  HKL  DETERMINE  A  VECTOR  H. 


T1«3*C*ALPHA»ALRHA 

T2«2«3*S(2>»2»C*AL®H*»5<J> 

Tj»A»S<1)»»2*J*3(23*»i*C*S(3>»»J”1 

TE.V.?*T2*T2*4*T1»TS 


C  NO  SOLUTION  HE  AH  S  :LLI*SOIO  *N*  NOT  TOUCHED  i Y  LINE 

C  CF  SIGHT  PAT. 

C 


IFUEmp.LT.O.O)  GO  TO  2 
T£*B*SJ»T<T£«P> 

c 

C  FIND  two  possible  r  vectors  6ECAUSE  a  PAT  entering 

C  A  3CLI0  UST  ALSO  LEAVE. 


*<2,1 )*(T2*TENP)/(2«T1) 

H<:,1)«AL®HA»*.<2,1) 
1<2»2>«<T2-T£PP)/C2«T1> 
HC 1,2) *0.2 
*1CJ,2>«ALP>iA**<2,2> 

J  F LA  j *3 
RETURN 
2  JFLA3*1 
hET'JRN 
END 


HIS 

1636 

1632 

1633 
163* 
164C 

1641 

1642 

1643 

1644 
1 64  J 
1046 
1647 
1 64,1 
1  o4  ? 
1633 

1631 

1632 

1633 

1634 
1633 
1636 
1657 
1653 

1659 

1660 
1661 
1662 
1  ooS 

1664 

1665 

1666 
1667 
1663 
166? 
1670 
1  <  71 

1672 

1673 

1674 

1675 

1676 

1677 
1673 


\ 
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suioctinc  t<Hu/A/i*C/S*««jM.AG)  m?» 

::*(NiI3N  »j(S)/S(3)/<'(3/2>  1*40 

AMI  r,*u  1611 

:  . . . . . .  me? 

:  i«si 

:  SOI  Vi  TH!  EJJAT10N  n»  *LLI*li  *  V  •*£*(  K»0  AN5  163i 

:  to  t:  tin>  *  *oi\t  0-4  r*c  *llj/»s:  mat  his 

;  .ill  ;et:a*in£  a  vtctja  •  .  us* 

:  us? 

:  . . . .  mi 

T 1  ■  C  MS* 

t:o*c»su>  moo 

!?•»•*  Cl  ).0»J.SCJ1*00.5<  3»..?-1  1 6  ?1 

T!*»«T?.T»-A.T1.TJ  1692 

£  . . . . . . .  1691 

:  i  *94 

:  NO  SCLJTI3N  NEANS  ELLI’GOIO  '  N*  NOT  TOJCMED  *V  LINE  1*95 

:  Of  SIGHT  AAV.  1696 

C  16«9 

:  . .  1696 

IfCTEHP.LT.O.O)  GC  TO  2  1699 
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APPENDIX  A 

H1D0EN  LINE  PROBLEM  BETWEEN  TWO  ELLIPSOIDS 


Ellipsoids  are  plotted  as  a  set  of  contour  lines.  These  lines  consist  of  a 
series  of  short  vectors  that  are  sequentially  plotted  to  form  a  contour 
line.  The  hidden  line  problem  can  be  reduced  to  the  problem  of  finding 
what  vectors  are  hidden  from  the  viewpoint.  This  problem  can  be  broken 
down  further  If  It  Is  assumed  that  each  vector  Is  short  enough  that  only 
the  point  representing  the  head  of  the  vector  needs  to  be  checked;  thus  the 
problem  reduces  to  checking  points  to  see  If  they  are  hidden  from  the 
viewpoint.  Figure  A.l  shows  the  coordinate  systems  and  vectors  used  to 
solve  the  hidden  point  problem.  The  following  equations  are  used  to  solve 
the  hidden  point  problem  using  ellipsoids. 


P21  *  [Di]  P2 


PIT  *  COi]  pT 
r2 1  *  [Di]  [D2]T  r2 

5  -  F21  *  P2T  +  PIT  *  0 


£  -  p2i  +  ^21  -  P 1 1 


t 


In  ellipsoid  No.  1  system 


[Di]  (P^)  +  [D j]  [D2]T  F*2  -  CD  J  Pt 


VPi  *  [Di ]  VP 
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z 


Figure  A.l  Coordinate  Systems  and  Vectors  Used  to  Solve  the 
Hidden  Line  Problem 


^  -  vFT  ♦  pTT  -  o 
*  ■  vpT  -  PIT 

^1n  ellipsoid  No.  1  system  "  ^  “  ^i3  Pi 

?  ♦  fi  ■  $  -  $  ■  £  -  ft  ^  +  ^ 

using  the  ellipsoid  equation 


r  •  [A]  r  *  1 

?T  [A]  r  *  1 


Since  P  and  M  are  In  the  same  direction,  P  +  M  Is  also  in  in  the  same 
direction,  then  a  unit  vector  in  the  direction  of  M  can  be  obtained  from 

P  +  M  =  S  -  V 


therefore 


where 


51  -  Vi 

52  -  V2 

" . HAG' 

s3  -  V3 
■  "HAG' . 


mag  =  [(Si  -  Vi)2  +  (S2  -  V2)2  +  (S3  -  v3)2]1/2 


giving  u  in  the  direction  of  M  must  obey  the  following  relationship, 


Pl  Mi 
P2  ^2 


Pi  Ml 
P3  WT 


P2 

M2  =  Ml  —  =  Mi 
1  ui  iot 


P  3 

M3  =  Mi  — -  =  Mlg 


back  to 


{?  -  M)  [A]  (?  -  it)  -  1 
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S  -  M  = 


/Si  *  Mi 

52  *  <*Mi 

53  -  0Mij 


therefore 


($  -  ti)T  CA]  (?  -  tf)  =  1 


expands  to 


“  (Si  -  Mi)2  +  (S2  -  oMi)2  +  —  (S3  -  0Mi)2  =  1 

a  b2  c2 


A  (Si  -  Mi)2  +  B(S2  -  aM  1)2  +  C(S3  -  3M ! ) 2  =  1 


ASi 2  +  AM12  -  2AS  iM  1  +  BS22  +  Ba2Mi2  =  2BS2cMi  +  CS32  + 
C02Mi2  -  2C0S  3M 1  =  1 


(A  +  Ba2  +  C02)  Mi2  -  (2ASi+  2BSZ  +  2CgS3)  Mi  +  AS12  +  BS22  + 
CS32  -1=0 
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let 


Tl  *  A  +  Ba2  +  C02 

T2  =  -  (2AS i  +  2BS2a  +  2C6S3) 

T3  *  AS22  +  BS22  +  CS32  -1 
-T2  ± 

Mi  =  - 

only  T22  -  4TiT3  must  be  evaluated  to  test  for  Mi  on  the  ellipsoid. 

4 

If  Mi  is  real,  then  M  is  given  by 

/Mi 

ft  =  aMi 

BMi 


If  there  is  a  real  M,  then  it  must  be  determined  if  both  Ms  are  in  the 

4- 

opposite  direction  of  P. 

♦  ♦ 

Case  No.  1— if  both  Ms  are  in  the  opposite  direction  of  P,  then  the  point 

is  not  hidden. 

Case  No.  2— if  either  ft  is  in  the  same  direction  as  P,  then  the  point  is 
hidden. 
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Using  the  dot  product  to  determine  the  relative  directions  of  M  and  P  gives 

?  •  =  <  0  ,  the  point  Is  not  hidden. 

♦  t  ♦ 

P  =  S  -  V  -  M 

(Si  -  Vi  -  Ml  S2  -  V2  -  M2,  S 3  -  V3  -  M3) 

i$  •  =  S1M1  -  V1M1  -  M1M1 

+  s2m2  -  v2m2  -  m2m2 

+  S3M3  -  V3M3  -  M3M3 


If  |fi|  is  close  to  zero,  the  points  on  each  ellipsoid  are  very  close 
together,  and  it  is  not  necessary  to  hide  the  point  on  the  contour  vector. 

If  u1  equals  zero,  then  an  a  and  8  cannot  be  found  with  the  given 
expressions.  Therefore,  the  following  cases  must  be  examined. 

Case  No.  1 


Hi  =  0  U2  =  0 

Mi  =  0 
M2  =  0 

(§  -  tf)T  A($  -  fl)  -  1 


Si 


M 


\ 


S-fi.  s2 
\s3 


/ 
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AS!2  +  BS22  +  C  (S3  -  M3)2  -  1 
ASi2  +  BS22  +  CS32  +  CM32  -  2SC3M3=  1 

T1  =  C 
T2  =  -2CS3 

T3  =  ASi2  +  BS22  +  CS32  -  1 


Use  these  T  values  back  in  the  quadratic  formula  used  earlier. 
Case  No.  2 

ui  =  0  u2  *  0 

Mi  =  0 


$  -  fi  - 


/Sl  \ 

s2  -  m2 

[Sl  -  M3i 


M3  y3 
Wz  =  U2 

M  3  =  oM2 


a 


M3 

M2 


ASi2  +  BS22  +  BM22  =  2BS2M2  +  CS32  +  Ca2M22  -  2Cc*12S2  =  1 


T1  =  B  +  Ca2 


T2  =  2BS2  -  2CaS3 
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T3  =  AS i2  +  BS22  +  CS32  -  1 

I 

I 

1 

Use  these  T  values  In  the  quadratic  formula  and  proceed. 


m 


APPENDIX  B 

DISCUSSION  OF  EQUATIONS  USED  BY  PRJELR 

PRJELR  stands  for  project  ellipsoid  routine.  The  function  of  PRJELR  is  to 
circumscribe  a  projected  shadow  of  an  ellipsoid  with  a  rectangle.  The 
resulting  rectangle  Is  used  as  a  polygon  by  the  overlap  routines  to  deter¬ 
mine  what  objects  overlap  after  they  are  projected. 

GENERAL  APPROACH 


1.  Assume  that  all  ellipsoids  project  an  elliptical  shadow.  This 
is  a  good  assumption  when  the  viewpoint  is  far  away  from  all 
objects  and  the  viewpoint  Z  axis  is  nearly  directly  pointed  at 
all  ellipsoids. 


2.  Find  three  radial  vectors  of  the  ellipsoid  pointing  to  a  surface 
point  that  forms  the  contour  of  the  projected  shadow.  Add  three 
further  conditions— one  radial  vector  is  in  the  X-Z  plane  of  the 
viewpoint  coordinate  system,  another  radial  vector  is  in  the  Y-Z 
plane  of  the  viewpoint  coordinate  system,  and  one  radial  vector 
is  in  the  X  =  Y  and  Z  plane  of  the  viewpoint  coordinate 
system.  These  planes  are  defined  as  if  the  viewpoint  coordinate 
system  were  at  the  center  of  the  ellipse.  This  gives  three 
radial  vectors  as  indicated  in  Figure  B.l,  rj,  r2,  r3. 


3.  Project  all  three  vectors  onto  the  projection  plane  and  solve 
for  an  ellipse  matrix. 

4.  Circumscribe  the  resulting  ellipse  with  a  rectangle. 


Take  [A ( 3 , 3 ) ]  ellipsoid  matrix  and  transform  to  [A' (3,3)]  in  the 
viewpoint  coordinate  system. 


A'  =  [DVP]  [D]T  [A]  [D]  [DVP]T 
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VIEWPOINT 

COORDINATE 

SYSTEM 


VIEWPOINT 

COORDINATE 

SYSTEM 


Figure  B.l.  Three  Radial  Vectors 
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,9  r»,\ 


where 


1 


[OVP]  =  direction  cosine  matrix  that  transforms  from  the 
inertial  to  the  viewpoint  frame  of  reference. 


direction  cosine  matrix  that  transforms  from  the 
Inertial  to  the  ellipsoid  frame  of  reference. 

Find  the  following  three  vectors 


/rx\ 

\  .  1 

1°  \ 

\  .  i 

tr'\ 

-► 

r  = 

'  0 

II 

ft 

ry 

"54 

II 

1  rzj 

1  ' 

\rz) 

1  1 

lrz/ 

and 

CD]  = 


that  also  have  the  properties  of  being  radial  vectors  defined  by  and 
that  the  associated  vector  $  from  the  viewpoint  to  the  tip  of  r  is  normal 
to  the  normal  vector  for  the  point  defined  by  r  on  the  ellipsoid.  See 
Figure  9. 


> 

n 


i*  =  0 


rTA'r  =  1 


-► 

n 


=  A'r 


A'r  •  ?  =  ?  •  A'r  =  0 


{SS  +  r)TA'r  =  0 
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S$TA'r  +  rTA'r  *  0  rTA'r  =  1 

s!TA'f  -  -  1 

This  Is  the  equation  that  subroutine  SOLVR  uses  to  return  components  for  i* 
that  represent  rj,  7 2,  ^3* 

Let 

Case  No.  1  be  ri 
Case  No.  2  be  r^ 

Case  No.  3  be  F3 
in  all  three  cases 

/  A'u 

(SSX>  SSY,  ssz)  A'21 

\  A '  3 1 

expanding  the  left  hand  part  gives 

(SSxA'n  +  SSyA* 21  +  SSzA’31,  SSxA'12  +  SSyA'22  +  SSZA'32 
SSxA'13  +  SSyA  '23  +  SS^A 1 3 3 )  r  =  -1 
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With  each  case  of  r  the  expression  reduces  to 


(<*iSSx  +  o^SSy  +  a3 SSZ)  r^y  +  (o4SSx  +  o5SSy  +  a6SSz)  rz  =  -1 
a4  =  A' 13 

a5  =  A'23  true  for  all  cases 

a6  =  A' 33 

Case  No.  1 
ai  =  A'n 

“2  =  A 1 21 

a3  =  A 1 3 1 

Making  another  substitution,  let 

=  (o^SS^  +  a2SSy  +  a3SSz) 

e2  ■  (a4SSx  +  a5SSy  +•  a6SSz) 

Then  rX/Y  can  be  solved  for  by  the  following  equation 

e2  1 

rx/Y  =  '  FT  rz  ‘  eT 


Case 

No.  2 

Case 

No 

.  3 

«i  “ 

CM 

“1  * 

A' 

11  + 

A' 

12 

“2  = 

CSI 

CM 

< 

“2  = 

A' 

21  + 

A1 

22 

“3  * 

A'  32 

03  = 

A' 

31  + 

A' 

32 
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now  ri,  r2,  and  P3  Is  written 


Since 


4 

r\ 


I  @2 

37 


rz  " 


0 

r-, 


4:' 


4- 

T2  * 


/  ° 

02 


TT  r2  ‘  TT 


^7  rz  -  ^7 


4 

r3  = 


/  »  2  ! 
37  ^  01 

\  ^ 


rTA'r  =  1 


ri  can  be  found  for  any  one  of  the  three  cases,  if  the  abo.ve  expression  is 
expanded,  the  following  general  form  results. 


r‘rX/Y2  +  T0rX/YrZ  +  T0fZ2  *  1 
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Case  No.  1 


Y 1  s  A'  11  Y  2s  2A ' 13  Y  3a  A'33 


Case  No.  2 

Yl  s  A1 22  Y2  =  2A1 23  Y3  =  A1 33 


Case  No.  3 


Yl  =  A'u  +  2A‘12  +  A'22  Y2  =  2(A'13  +  A'23)  y3  =  A'33 


now  by  substituting 


■  X/Y 


02 

FT  rZ 


1 

Pi 


results  in  the  following  expression 


fBi  ,  2  &2 

Yl  "  *  FT  rz  '  FT  +  Yz  '  FT  rZ 


i_r 

Bi  z 


Y3r72  =  1 


expanding  and  combining  like  terms  yields 


/ e2  \  2 

Y1  FT  +Y3 


&2  B2 

Y2?_rz2  +  2T 


1  a 


i_V 


frz  +  Yi  — 1  -1  =  0 
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and  constant  term  be 


Let  the  coefficients  for  the  terms  r^, 
represented  by  Tl,  5T2,  and  T3,  respectively. 


T2  ±  T 22  -  (4T1T3)172 
rz  - - 2TT - 

$2  $2 

rZ/Y  =  ‘  87  rZ  ‘  37 

Values  for  rx,  rY>  and  r^  can  be  obtained  depending  upon  what  case  is  being 
sol ved. 

SOLVR  subroutine  will  return  vector  components  for  any  of  the  three 
cases. 

Definition  of  call  to  SOLVR  subroutine 


CALL 

SOLVR 

(Ai,  A2,  A3, 

A4  * 

A5» 

A6»  a7,  a8, 

ss. 

Rl, 

R3) 

Case 

No. 

.  1 

Case 

No. 

_2 

Case 

No. 

_2 

Ai  = 

A1 

11 

A 1  = 

A' 

12 

Ai  = 

A' 

n  + 

A'  12 

II 

CM 

C 

A1 

21 

II 

CM 

< 

A' 

22 

A  2  = 

A' 

21  + 

A '  22 

11 

CO 

«£ 

A' 

31 

A3  = 

A' 

32 

A3  = 

A' 

31  + 

A  32 

At*  = 

A' 

13 

A4  = 

A' 

13 

-P 

II 

A' 

13 

As  • 

A' 

23 

As  = 

A' 

23 

As 

=  A 

'  23 

A6  = 

A' 

33 

A6  = 

A' 

33 

o*» 

11 

A' 

33 

J 
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a7  =  A' u 

A7  =  A'22 

A7  =  A*  n  +  2A1 1?  +  A' 

CO 

II 

GO 

c 

^8  =  A ' 23 

^8  =  A1 13+  A'  23 

ss  =  ss 

SS  =  SS 

SS  =  SS 

R1  =  X  component 
of  ri 

RI  =  Y  component 
of  r2 

RI  =  X  and  Y  component 
of  r3 

R3  =  Z  component 
of  ri 

R3  =  Z  component 
of  r2 

R3  =  Z  component 
of  r  3 

After  SOLVR  constructs  all  three  cases  for  vector  these  three  vectors 
are  used  to  find  an  ellipse  matrix  [a]  on  the  projection  plane.  The  pro¬ 
perties  of  the  three  ?  vectors  are  such  that  they  satisfy  the  three- 
dimensional  ellipsoid  and  when  projected  satisfy  the  ellipse  on  the  projec 
tion  plane. 

Also  the  pr  , Arties  are  such  that  the  coefficients  of  the  ellipse  matrix 
can  be  found.  This  is  accomplished  by  SOLVA  subroutine. 


Projection  of  r  onto  the  Projection  Plane 


r'  =  p'  -  S'  p 


*  rX  **Y  *  rY 
Sx  +  i"z  *  +  rZ 


S'  = 


r'  = 


'SX  +  rx 


S^  Sy  +  Sy 


^TT7  +  r^ 
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r  ,Tar '  *  1 


Ca] 


/a 1 1  a12 


a2i  a22 


Since  there  are  three  P  vectors  and  only  three  of  the  four  components  of 
[a]  are  Independent,  the  components  of  [a]  are  obtained  by  solving  three 
equations  simultaneously  In  subroutine  SOLVA.  SOLVA  returns  the  components 
of  [a]. 


To  circumscribe  the  ellipse  with  a  rectangle,  the  major  and  minor  axis 
vectors  must  be  found.  These  vectors  are  found  by  solving  for  the  eigen¬ 
vectors  of  [a], 

aP  =  aP 


This  condition  Is  true  only  for  the  vectors  that  represent  the  major  and 
minor  axis  of  the  ellipse. 


/ all  a12  ' 

1 

'rx’ 

\ a12  a22  t 

'  ' 

A l 

/ ai irx  +  ai2rY' 

/*rx 

^  a12rX  +  a22ry^ 

\  rY 

A 


rx\ 


Ar^  -  a 1 1 r^  -  a^Ty  =  0 


Ary  _ai2r^  -  a22ry  a  0 


(X  -  an)  rx  -  ai2ry  ■  0 


[1] 


(X  -  a22)  ry  -  ai2rx  •  0 


[2] 


The  only  w«y  No.  1  and  No.  2  can  be  zero  Is  If 


/  a12 


\Xi  -  an 


r2  = 


fx2  -  a22' 


an 


/r* 


\rY 


These  are  the  major  and  minor  axes  vectors. 


These  vectors  must  also  satisfy  the  ellipse  equation 


ar  =  Xr 


+T  + 
r'ar 


=  1 


rTXr 


=  1 
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since  X  is  a  scalar 


Xr  r 


2 


x 


x  + 


1 

X 


1 

X 


Both  eigenvalues  can  be  found  by  solving 


x  -  an 

-ai2 

-ai2 

X  - 

X  = 


ail  +  a22  ±  [(an  +  a22)2  -  Mana22*  a^2)]1^2 


Using  these  two  eigenvalues,  ri  and  r2  are  determined  and  must  be 
normalized  by  the  relation 


These  equations  are  used  to  circumscribe  the  ellipse  with  a  rectangle. 


Pi  =  ?t  +  r2 


P2  =  -r i  +  r2 
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P3  -  -ri  -  r2 


P4  -  ri  -  r2 

CONVEC  (1 ,1,K)  -  Fj-Pt*rJ-rt-rt-rJ»  -2  Ft 

CONVEC  (1 ,2,K)  -  Pj-Pj-rJ-rT-Fj+rT  -  -2rJ 

CONVEC  (I ,3,K)  *  K  -  P3  s  rt  -  rj  +  rj  +  Ft-  2pJ 

CONVEC  (I ,4,K)  =  pt-Pj=rt+rt-Ft+Ft  =  2rt 

I  =  1,2  for  vectors  on  the  projection  plane 
K  1,2, 3, 4,  ....  90 

K  is  an  ellipsoid  number  or  polygon  number. 
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APPENDIX  C 

INTERSECTION  OF  A  THREE  SPACE  VECTOR  AND  PLANE 


Figure  C.l  shows  the  typical  vector  problem  that  represents  the 
Intersection  of  a  three-space  vector  and  a  plane.  The  problem  Is  to  deter¬ 
mine  If  the  Intersection  point  lies  along  vector  ?  or  beyond  the  tip  of 
K  The  figure  shows  that  If  Tau  Is  less  than  one,  the  Intersection  point 
lies  between  the  two  points  VP  and  PT. 

A 

Tau  Is  found  with  the  following  equations  (where  N  Is  the  normal  unit 
vector  to  the  plane). 


N  •  (f  *  0 


e  =  n  -  ? 


N  •  (ri  -  ?)  «  0 

N  •  rx  -  N  •  P  -  0 

x  -  j  '  * 
n  •  i* 

t  >  1  point  Is  not  blocked  by  the  Intersection  point. 


APPENDIX  D 

INTERSECTION  OF  LINE  SEGMENTS  IN  A  PLANE 
Given  two  line  segments, 

P  i f*  2  sod  P  3P4 

where 

Pi  *  («i .  yi ) 

Consider  all  parallel  lines  to  be  nonintersecting. 

The  Regular  Configuration  --  Neither  of  the  line  segments  Is  vertical. 
The  line  which  contains  P1P2  has  the  equation 

y2*yi  y-yi 

x2  -  X!  *  X  -  Ki 
which  simplifies  to 


y  -  (x  -  xi)  mi  ♦  yi 


where 


mi 


yi  -  yi 

X2  -  Xi 


Likewise  the  line  containing  P3P4  has  the  equation 
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y  -  (x  -  x3 )  m2  +  y3 


where 

y**  -  y3 

"2 

Since  at  the  point  of  Intersection 

Po  *  (xo.  yo) 

yo  ■  (x0  -  xx)  m3  +  yi 

yo  =  (x0  -  x3)  m2  +  y3 

then  equating  and  solving  for  xo  yields 

y3  -  yi  +  m ix i  -  m2x3 

x0  = - 

mi  -  m2 

To  determine  If  the  point  of  Intersection  of  the  two  lines  Is  on  each  line 
segment,  note 

Wj  <■>'.•  P1  ♦  *(pj  -  V  for0<t<1 

f 
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Then  let 


t 


xo  -  xi 

X2  -xi 


xo  -  x3 

s  ■  — - 

X4  -  X3 

then  if  0  <  t  <  1  and  0  <  s  <  1  ,  there  Is  Intersection. 

But,  since  we  will  say  that  the  two  segments  do  not  Intersect  If  the  point 
of  Intersection  Is  one  of  the  endpoints,  then  1f0<t<land0<s<l, 
there  Is  Intersection. 

If  one  of  the  lines  Is  vertical,  the  regular  procedure  will  not  work  since 
there  will  be  a  zero  denominator  In  one  of  the  n^'s.  Therefore,  provided 
the  nonvertical  segment  Is  not  horizontal,  make  the  substitution 


*  yv  x1  for  each  Pi  -  (x^  y^ 


Using  the  P^ '  endpoints,  the  regular  procedure  will  then  determine  If  there 
Is  an  Intersection. 

The  only  case  not  covered,  so  far.  Is  the  case  where  one  segment  Is 
vertical  and  the  other  Is  horizontal.  Without  loss  of  generality, 
assume  F jfPJ  Is  vertical  and  Fp^  Is  horizontal. 

In  this  case 

t 

Po  -  (xi.  y3) 
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so 


(XI.  y3 )  ■  (xi.  yi)  +  t  C(xi,  y2)  -  (xi,  y i)3  -> 


t 


y3  -  yi 
y2  -  yi 


likewise 


C(X1.  y3)  ■  ( x 3 »  y3)  +  s  (x4,  y3)  -  (x3»  y3)]  ■> 


xi  -  X3 
s  *  x4  -  X3 

then  If  0  <  t  <  1  and  0  <  a  <  1  ,  there  Is  an  Intersection. 
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